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ABSTRACT 


The  branched  crack  problem  for  both  an  isotropic  and  anisotropic  material  is  solved 
using  the  method  of  dislocations,  and  the  stress  intensity  factors  and  T-stress  in  front  of  the 
branched  crack  are  evaluated  numerically .  The  T-stress  based  fracture  criteria,  developed 
by  Cotterell  and  Rice  for  a  flat  crack,  is  modified  by  incorporating  an  experimentally  deter¬ 
mined  critical  T-stress  value  (Tent).  Based  on  this  Tcrit,  a  modified  T-stress  based  fracture 
criteria  is  proposed.  This  criteria  is  applied  to  the  branched  crack  and  the  direction  of  growth 
of  the  branched  crack  is  discussed  by  comparing  the  theoretically  evaluated  T-stress  values 
with  some  available  experimentally  determined  Tcrit  values. 

The  solution  for  the  branched  crack  problem  is  obtained  in  terms  of  a  singular  integral 
equation  which  is  solved  using  three  different  numerical  schemes,  the  merits  of  which  are 
discussed.  The  nature  of  the  stress  singularity  at  the  reentrant  wedge  comer  of  the  branched 
crack  is  analyzed  and  is  verified  for  the  isotropic  case.  It  is  shown  that  the  T-stress  and  the 
stress  intensity  factors  are  insensitive  to  the  order  of  the  singularity  assumed  at  the  reentrant 
wedge  comer  of  the  branched  crack  in  either  an  isotropic  or  anisotropic  material.  The  T- 
stress  for  the  isotropic  case  is  obtained  in  terms  of  applied  load,  kink  length  and  kink  angle. 
For  the  anisotropic  case  the  T-stress  also  depends  on  the  relative  stiffness  properties  of  the 
fibre  and  matrix.  For  a  uniaxial  loading  case,  by  applying  the  modified  T-stress  based  crite¬ 
ria  to  the  branched  crack,  it  is  demonstrated  that  for  a  short  kink  length  the  kink  will  turn  from 
its  initial  direction  and  realign  with  the  main  crack.  If  the  loading  is  biaxial  then  the  kink 
growth  direction  depends  strongly  on  applied  transverse  stress  Oxx°°-  For  a  longer  initial 
kink  length  the  kink  growth  direction  depends  on  the  kink  angle,  loading  and  relative  stiff¬ 


ness  of  the  fibre  and  matrix. 


TABLE  OF  CONTENTS 


Page 

TITLE  PAGE .  i 

ABSTRACT .  ii 

LIST  OF  TABLES .  vi 

LIST  OF  FIGURES .  vii 

CHAPTER 

I.  INTRODUCTION .  1 

Background  .  1 

Literature  Review  On  the  Branched  Crack 

Problem-Isotropic  Materials  .  3 

Conformal  Mapping  Techniques  .  4 

Dislocation  Technique  .  5 

Mellin  Transform .  6 

Perturbation  Technique .  7 

Literature  Review  On  the  Branched  Crack 

Problem-Anisotropic  Material .  8 

Review  of  Salient  Fracture  Criteria-Isotropic  Material  ....  10 

Review  of  Fracture  Criteria- Anisotropic  Materials .  11 

How  is  the  Present  Study  Different  from  Other 

T-stress  Based  Studies  ? .  14 

An  Overview  of  this  Study .  14 


II.  ISOTROPIC  MATERIAL .  16 

Analytical  Formulation  .  16 

Description  of  Dislocation  Function  .  16 

Method  of  Superposition  Explained .  19 

Formulation  of  the  Problem .  19 

Verification  of  Existence  of  Williams’  Singularity 

at  the  Internal  Wedge  Comer .  24 

Extension  of  Theorems  for  Singular  Integral 

to  Bounded  Integrals  .  24 

Proof  For  Verification  of  Williams  Singularity 

.  25 


IV 


Table  of  Contents  (Continued) 

Page 

III.  ANISOTROPIC  MATERIAL  .  32 

Analytical  Formulation  .  32 

Development  of  Dislocation  Function .  32 

Formulation  of  Stress  Functions  for  the  Main  Crack 
Which  Opens  Up  Due  to  the  Presence 

of  a  Dislocation  .  34 

Formulation  of  Stress  Functions  for  the  Main 
Crack  Which  Opens  Up  Due  to 

Loads  at  Infinity .  34 

Formulation  of  the  Singular  Integral  Equation  .  35 

IV.  EVALUATION  OF  STRESS  INTENSITY  FACTORS 

AND  T-STRESS  .  37 

Isotropic  Material .  37 

Evaluation  of  the  T-S tress  for  Isotropic  Material  .  39 

Anisotropic  Material  .  44 

Evaluating  the  T-stress  for  Anisotropic  Material  .  45 

V.  RESULTS  AND  DISCUSSION-I .  37 

Isotropic  Material .  48 

Development  of  Fracture  Model  .  55 

Kink  Closing  Angle .  66 

VI.  RESULTS  AND  DISCUSSION-H .  68 

Anisotropic  Case  .  68 

VII.  CONCLUSIONS .  83 

APPENDICES .  85 

A.  ISOTROPIC  MATERIAL .  86 

Analytical  Formulation  for  a  Dislocation  Function .  86 

Formulation  of  Stress  Functions  for  the  Main  Crack 
Which  Opens  Up  Due  to  the  Presence  of  a 

Dislocation  .  88 

Formulation  of  the  Solution  for  a  Infinite  Plane 

with  a  Crack .  99 

Formulation  of  the  Solution  for  the  Kinked  Crack .  100 


B.  EVALUATION  OF  SINGULAR  INTEGRALS 


102 


Table  of  Contents  (Continued) 


Page 


Formula  .  102 

Some  Typical  Singular  Integrals  .  103 

C.  ANISOTROPIC  MATERIALS  .  106 

Formulation  of  Governing  Equations  for  an 

Anisotropic  Elastic  Body  .  106 

Analytical  Formulation  For  Dislocation  Function .  109 

Formulation  of  Stress  Functions  for  the  Main  Crack  Which 

Opens  Up  Due  to  the  Presence  of  a  Dislocation  .  112 

Formulation  of  Stress  Functions  for  the  Main  Crack  Which 

Opens  Up  Due  to  Loads  at  Infinity  .  121 

Formulation  of  the  Singular  Integral  Equation  .  123 

D.  EVALUATION  OF  SOME  CONTOUR  INTEGRALS .  126 

E.  NUMERICAL  TECHNIQUES  .  134 

Method  I:  Solution  by  Gauss-Chebychev  Integration 

Formula .  135 

Method  11:  Solution  by  Jacobi  Polynomials .  135 

Method  nr.  Solution  by  Double  Series  Jacobi 

Polynomials .  136 

LITERATURE  CITED  .  138 


LIST  OF  TABLES 

Table  Page 

1.  Load  Case  1  (No.  of  Quadrature  Points  =  121,  Method  I) .  50 

2.  Load  Case  1  (No.  of  Terms  in  Polynomial,  n  =  40,  Method  II) .  50 

3.  Load  Case  1  (No.  of  Terms  in  Polynomial,  ni=  6,  n2  =4, 

[see  E.  14]  Method  m) .  50 

4.  Two  Roots  of  the  Williams’  Characteristic  Equations  (36)  and  (37) 

Used  in  Numerical  Analysis  .  51 

5.  Load  Case  2  (No.  of  Quadrature  Points  =  121)  .  51 

6.  Load  Case  3  (No.  of  Quadrature  Points  =  121)  .  51 

7.  Comparison  of  Results  with  Mellin  [34]  and  Isida  [86] 

for  2c// =  oo  (Method  I)  .  54 

8.  Comparison  of  Results  with  Mellin  [34]  and  Isida  [86] 

for  2c//  =  oo  (Method  IH)  .  54 

9.  Comparison  of  T-stress  Obtained  by  the  Three  Methods  for 

Unit  Load  at  Infinity  (//2c=0.005)  .  54 

10.  Comparison  of  Eq.  (76)  with  Numerical  T-stress 

Eq.  (67)  (2c//  =  oo) .  55 

11.  Roots  of  Bogy’s  [87]  Characteristic  equation  (99) .  71 

12.  Ki,  Kn  and  T-S tress  Calculated  Using  the  Three 

Methods  (//2c  =  0.005,  E1/E2  =1.1)  .  71 


LIST  OF  FIGURES 


Figure  Page 

1 .  Geometry  of  the  Branched  Crack  Problem .  2 

2.  Edge  Dislocation .  17 

3.  Screw  Dislocation .  17 

4.  In  Plane  Shear .  18 

5.  Normal  Stress  .  18 

6.  Anti  Plane  Shear .  18 

7.  Superposition  Steps  in  Solving  the  Kinked  Crack  Problem .  20 

8.  Tangential  Jump  .  21 

9.  RadialJump  .  21 

10.  Geometry  of  the  Wedge  as  Defined  in  Timoshenko  et.al.  [84] .  30 

11.  Geometry  Of  Anisotropic  Branched  Crack  Problem .  33 

12.  Configuration  Describing  Near-Kink-Tip  Stresses  for  Evaluating 

Stress  Intensity  Factors .  37 

13.  Near  Crack-Tip  Stresses  for  the  Branched  Crack .  40 

14.  Components  of  Remote  Loads  Parallel  to  and  Perpendicular 

to  the  Kink .  44 

15.  Loading  Condition  1  .  52 

16.  Loading  Condition  2  .  52 

17.  Loading  Condition  3  .  53 

18.  Comparison  of  an  Arbitrary  cracked  Coupon  with  an 

Infinite  Center-Cracked  Panel  (Richardson  [94]) .  57 

19.  Definition  of  Class  I  and  Class  n  fracture .  59 

20.  Variation  of  Normalized  Kj  with  Normalized  Kink  Length  /  (k=0)  . .  62 


vm 


List  of  Figures  (Continued) 

Page 


21.  Variation  of  Normalized  Ki  with  Normalized 

Kink  Length  /  (k=5) .  62 

22.  Variation  of  Normalized  Kn  with  Normalized  Kink 

Length  I  (k=0) .  63 

23.  Variation  of  Normalized  Knwith  Normalized  Kink 

Length  /  (k=5) .  63 

24.  Variation  of  Normalized  T-stress  with  Normalized 

Kink  Length /(k  =  0) .  64 

25.  Variation  of  Normalized  T-stress  with  Normalized  Kink 

Length  /  (k  =  5) .  64 

26.  Geometry  of  PMMA  Test  Coupons .  65 

27.  Kink  Tip  Closing  Angle  (0c)  as  a  Function  of 

Normalized  Kink  Length .  67 

28.  Geometry  Showing  Parameters  for  Anisotropic  Branched 

Crack  Problem .  69 

29.  Variation  of  Normalized  Kj  with  Normalized  Kink  Length  / 

(k  =  0,Ei/E2=  10,  9^  =  90) .  74 

30.  Variation  of  Normalized  Ki  with  Normalized  Kink  Length  I 

(k  =  0,Ei/E2=  1.08,  0^  =  0) .  74 

31.  Variation  of  Normalized  Kj  with  Normalized  Kink  Length  I 

(k  =  O,Ei/E2=lO.O,0m  =  O) .  75 

32.  Variation  of  Normalized  Kn  with  Normalized  Kink  Length  / 

(k  =  O,Ei/E2=lO,0m  =  9O) .  75 

33.  Variation  of  Normalized  Kn  with  Normalized  Kink  Length  I 

(k  =  O,Ei/E2=l.O8,0in=O) .  76 

34.  Variation  of  Normalized  Kn  with  Normalized  Kink  Length  I 

(k  =  O,Ei/E2=lO,0na  =  O) .  76 


IX 


List  of  Figures  (Continued) 

Page 

35.  Variation  of  Normalized  T-stress  with  Normalized  Kink  Length  / 

(k  =  0,Ei/E2=  10,0^  =  90) .  78 

36.  Variation  of  Normalized  T-stress  with  Normalized  Kink  Length  / 

(k  =  0,  E1/E2  =  1 .08,  0m  =  0) .  78 

37.  Variation  of  Normalized  T-stress  with  Normalized  Kink  Length  I 

(k  =  O,Ei/E2=lO,0m  =  O) .  79 

38.  Variation  of  Normalized  Kj  with  Material  Angle  0m 

(k  =  0,Ei/E2=  10) .  79 

39.  Variation  of  Normalized  Kn  with  Material  Angle  0m 

(k  =  0,  E1/E2  =10) .  80 

40.  Variation  of  Normalized  T-stress  with  Material  Angle  0m 

(k  =  0,Ei/E2=  10) .  80 

A.l.  Geometry  of  Dislocation  Function .  86 

A.2.  Geometry  of  Pressurized  Main  Crack  .  89 

A.3.  Definition  of  the  Function  X(z)  .  91 

A.4.  Contour  G .  91 

A.5.  Contour  G  (L1+L2) .  94 

C.l.  Geometrj' of  Dislocation  Function  .  110 

C.2.  Geometry  of  Pressurized  Main  Crack .  113 

C. 3.  Geometry  for  the  Function  X(z)  .  116 

D. l.  Contour  Assumed  for  Evaluating  Integral  D.l .  126 

D.2.  Contour  Assumed  for  Evaluating  Integral  D.20 .  130 

D.3.  Contour  Assumed  for  Evaluating  Integral  D.26,  D.33 .  131 


CHAPTER  I 


INTRODUCTION 

Background 

Linear  Elastic  Fracture  Mechanics  (LEFM)  proposes  that  for  brittle  materials  a  single 
parameter,  the  stress  intensity  factor  (Kjc),  is  sufficient  for  predicting  fracture.  It  is  also  im¬ 
plied  that  the  Kic  is  a  property  of  the  material  similar  to  Young’s  modulus  and  the  ultimate 
strength.  Over  the  years  LEFM  has  been  applied  to  a  wide  range  of  materials  with  reasonably 
good  success.  But  it  has  been  observed  by  a  number  of  researchers  [1-4]  that  for  some  mate¬ 
rials  both  Kic  snd  crack  growth  direction  can  depend  on  the  coupon  geometry.  Albritton 
and  Goree  [5]  demonstrated,  through  a  number  of  experiments,  a  significant  influence  of  ge¬ 
ometry  on  Kic  for  silicon  carbide  reinforced  aluminum  coupons.  Three  types  of  coupons 
were  tested,  namely  compact-tension  (CT),  center-crack-tension  (CCT)  and  single-edge- 
notch  tension  (SENT).  It  was  observed  that  crack  turning  occurred  for  CT  coupons  while 
it  did  not  occur  for  SENT  and  CCT  coupons  with  the  same  configuration.  These  results  led 
to  further  investigation  of  the  problem  by  Richardson  and  Goree  [6] .  They  proposed  an  em¬ 
pirical  two-parameter  (the  parameters  are  Kic  ^nd  T-stress)  fracture  criterion  which  ac¬ 
counted  for  the  geometry  of  the  coupon.  The  validity  of  the  model  was  demonstrated  for 
Polymethyl  methacrylate  (PMMA)  [6].  This  model  however  did  not  account  for  turned 
cracks. 

In  the  present  study  the  crack  turning  behavior  in  isotropic  and  anisotropic  materials 
is  analyzed  and  the  T-stress  evaluated  in  front  of  a  branched  crack  (or  kink).  The  geometry 
of  the  branched  crack  problem  with  the  relevant  labels  is  shown  in  Figure  1 .  An  infinite  body 
(isotropic  or  anisotropic)  with  a  main  crack  of  length  2c  and  a  kink  (or  branched  crack)  of 
length  I  is  shown.  The  branched  crack  makes  an  angle  0  with  the  X  axis.  Bi-axial  loads  a 
and  ka  are  applied  remotely  along  the  Y  and  X  axis  respectively.  The  T-stress  is  defined 
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as  the  constant  term  in  the  asymptotic  stress  expansion  in  front  of  the  kink  or  the  stress  com¬ 
ponent  which  acts  parallel  to  and  along  the  kink.  Cotterell  and  Rice  [7]  evaluated  the  T— 
stress  in  front  of  an  unbranched  or  flat  crack  and  proposed  a  criterion  which  determines  the 
direction  of  crack  growth.  This  T-stress  based  criterion  is  modified,  and  the  modified  criteri¬ 
on  is  extended  in  the  present  study  to  the  case  of  the  branched  crack.  A  detailed  description 
of  the  development  of  this  T-stress  based  model  is  given  in  Chapter  V. 


In  order  to  study  the  influence  of  the  T-stress  on  the  branched  crack  it  is  necessary  to 
solve  the  branched  crack  problem  in  an  infinite  panel  for  both  isotropic  and  anisotropic  mate¬ 
rials.  In  the  remainder  of  this  chapter  a  review  of  the  literature  on  the  branched  crack  problem 
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is  presented.  A  brief  survey  of  related  fracture  models  for  predicting  crack  growth  and  stabil¬ 
ity  which  are  available  in  the  open  literature  is  also  given. 

Literature  Review  On  the  Branched  Crack  Problem-Isotropic  Materials 

A  large  volume  of  literature  is  available  on  the  subject  of  the  branched  crack  in  an  iso¬ 
tropic  material.  It  is  interesting  to  note  that  crack  branching  was  first  studied  for  the  dynamic 
crack  growth  problem  before  the  focus  shifted  to  crack  branching  during  quasi-static  crack 
growth.  The  term  quasi-static  (or  static)  is  used  to  describe  the  condition  when  the  crack-tip 
velocity  is  small,  and  there  is  no  sudden  change  in  velocity  and  the  dynamic  effect  that  it  pro¬ 
duces  can  be  neglected. 

Yoffee  [8]  was  one  of  the  earliest  to  study  crack  branching  during  dynamic  crack 
growth  in  a  brittle  material  like  glass.  She  was  the  first  to  propose  that  fracture  will  take  place 
in  a  direction  normal  to  the  maximum  tensile  stress  in  a  dynamic  crack  growth  situation. 
Based  on  this  assumption,  a  critical  velocity  of  crack  growth  was  determined,  above  which 
crack  branching  occurred.  Schardin  [9],  Clark  and  Irwin  [10]  and  Congleton  and  Fetch  [11] 
studied  similar  problems  both  analytically  and  experimentally.  Mostovy,  et  al.,  [4]  and  Carter 
[12]  studied  crack  branching  due  to  stress  corrosion  in  high  strength  steel. 

The  analytical  study  of  static  crack  branching  can  generally  be  grouped,  based  on  the 
methodology  used,  into  (1)  conformal  mapping  techniques,  (2)  dislocation  techniques,  (3) 
Mellin  transforms  and  (4)  perturbation  techniques. 

Using  conformal  mapping,  the  geometry  of  the  branched  crack  is  first  mapped  on  to 
an  equivalent  circle  and  then  the  mapped  problem  is  solved.  This  method  is  used  to  calculate 
the  stress  intensity  factors  and  the  energy  release  rates  at  the  branched  or  kinked  crack  tip. 
The  disadvantage  of  this  method  is  that  its  accuracy  decreases  as  the  kink  length  becomes 
infinitesimally  small. 

In  the  dislocation  technique  the  kink  is  modeled  as  an  edge  dislocation.  With  this 
method  it  is  possible  to  calculate  the  fracture  parameters  for  infinitesimally  small  branched 


cracks  also.  The  dislocation  technique  is  used  in  the  formulation  of  the  problem  in  the  pres¬ 
ent  study  since  a  range  of  kink  lengths,  from  infinitesimally  small  to  very  large  can  easily 
be  modelled. 

Mellin  transform  and  perturbation  techniques  can  also  be  used  to  solve  infinitesimally 
small  kink  length  problems.  Perturbation  techniques  can  be  used  to  obtain  solutions  for  finite 
kink  lengths  if  more  relevant  terms  are  included  in  the  series  solution.  Sometimes  a  com¬ 
bination  of  the  above  methods  are  used  to  solve  the  branched  crack  problem. 

In  the  next  few  paragraphs  a  review  of  the  literature  on  branched  crack  is  given,  which 
traces  the  development  of  studying  branched  crack  problems  in  an  isotropic  and  anisotropic 
material. 

Conformal  Mapping  Techniques 

Andersson  [13,14]  was  the  first  to  calculate  the  stress  intensity  factors  at  the  tips  of  a 
star  shaped  contour  in  an  infinite  tensile  sheet,  with  remote  loading,  using  conformal  map¬ 
ping  and  Muskhelishvili  [15]  complex  variable  methods.  Dudukalenko  and  Romalis  [16] 
used  conformal  mapping  techniques  to  map  a  finite  crack  with  a  branch  onto  a  circle.  They 
used  maximum  energy  release  criteria  to  determine  the  crack  turning  direction.  Palaniswa- 
my  and  Knauss  [17]  used  conformal  mapping  techniques  to  study  the  branched  crack  prob¬ 
lem.  They  compared  the  stress  and  energy  based  fracture  criteria  for  mixed-mode  fracture 
problems  and  concluded  that  the  energy  based  criteria  was  conservative  when  compared  to 
the  maximum  stress  criteria.  They  also  performed  experiments  relating  to  mixed-mode  frac¬ 
ture  and  discussed  the  results  in  light  of  their  analytical  solution.  Chatterjee  [18]  solved  the 
branched  crack  problem  by  using  conformal  mapping  method.  He  calculated  the  stress  in¬ 
tensity  factors  in  front  of  the  kink  and  the  main  crack  for  small  kink  length  of  the  order  of 
1/1000  th  the  main  crack  length.  For  smaller  kink  angles  it  was  reported  that  the  numerical 
scheme  became  unstable.  Chatterjee  also  was  the  first  to  show  that  the  stress  singularity  at 
the  kink  comer  was  indeed  the  Williams  [19]  singularity.  It  should  be  kept  in  mind  that  the 
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methodology  used  by  Chatterjee  for  solving  the  problem  was  conformal  mapping.  A  similar 
proof  for  verifying  that  the  singularity  at  the  kink  comer  is  that  given  by  the  Williams  [19] 
singularity  for  the  branched  crack  problem  (solved  by  dislocation  techniques)  will  be  dis¬ 
cussed  in  this  study.  Hussain,  et  al.,  [20]  used  the  conformal  mapping  technique  to  study  the 
branched  crack  problem  in  an  infinite  plane  and  calculated  the  crack  growth  path  based  on 
Griffith’s  energy  release  criteria  when  the  branched  crack  was  very  small  compared  to  the 
main  crack  length.  Kitagawa,  et  al.,  [21]  and  Kitagawa  and  Yuuki  [22]  also  analyzed  the 
problem  of  the  branched  crack  using  conformal  mapping.  They  examined  bent,  curved, 
branched  and  zig-zag  crack  geometries.  The  influence  of  biaxial  loading  on  the  stability  of 
crack  growth  was  also  addressed  by  them.  They  used  the  Erdogan  and  Sih  [23]  maximum 
hoop  stress  criteria  to  determine  the  angle  of  growth  of  a  kinked  crack.  With  their  model  they 
were  able  to  analyze  only  finite  kink  geometry. 

Dislocation  Technique 

Lo  [24]  was  the  first  to  solve  the  branched  crack  problem  using  the  dislocation  method 
and  he  obtained  the  stress  intensity  factors  for  the  entire  range  of  kink  length  varying  from 
infinitesimally  small  to  large  kink  length.  The  results  obtained  by  Lo  [24]  agreed  with  the 
results  of  Chatterjee  [18],  Kitagawa,  et  al.,  [21]  for  finite  kink  lengths  but  disagreed  with 
those  of  Dudukalenko,  et  al.  [16],  and  Hussain,  et  al.,  [20]  for  infinitesimally  small  kink 
lengths.  Karihaloo,  et  al.,  [25]  used  the  dislocation  method  to  reexamine  the  problem  ad¬ 
dressed  by  Cotterell  and  Rice  [7]  (who  had  obtained  first  order  solutions)  and  constructed 
a  second-order  solution  for  the  branched  crack.  For  an  infinitesimally  small  kink  they  were 
able  to  express  the  stress  intensity  factors  in  front  of  the  kink  in  terms  of  the  stress  intensity 
factors  in  front  of  the  main  crack  in  a  manner  similar  to  Cotterell  and  Rice  [7].  They  were 
also  able  to  obtain  the  correction  terms  for  large  kink  angles  by  means  of  which  the  results 
obtained  by  Cotterell  and  Rice  [7]  could  be  extended  for  kink  angles  more  than  40  degrees. 
It  should  be  noted  that  Table.  1  in  their  paper  (Karihaloo,  et  al.  [25])  contains  some  numerical 
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errors  but  can  be  corrected  if  the  formula  they  give  earlier  in  the  paper  is  used  to  calculate 
the  quantities  in  the  table.  Hayashi  and  Nemat-Nasser  [26],  using  dislocations,  expressed 
the  stress  intensities  and  the  energy  release  rate  of  the  kink  in  terms  of  the  stress  intensities 
existing  in  front  of  the  main  crack  and  some  coefficients  whose  numerical  values  were  given 
in  a  tabular  form.  Datsyshin  and  Savruk  [  27]  analyzed  the  problem  of  multiple  arbitrarily 
oriented  cracks  using  a  method  similar  to  the  dislocation  method.  They  obtained  the  solution 
in  terms  of  singular  integral  equations  where  the  unknown  was  the  slope  of  the  crack  opening 
displacement  just  as  in  the  dislocation  method.  They,  however,  did  not  study  the  case  when 
the  cracks  touch  each  other.  Theocaris  and  loakimidis  [28]  analyzed  the  symmetrically 
branched  crack  in  an  infinite  body  using  a  method  similar  to  that  used  by  Datsyshin  and  Sa¬ 
vruk  [27].  They  also  found  that  their  results  for  stress  intensity  factors  compared  well  with 
the  experimental  results  of  symmetrically  branched  cracks  which  was  determined  using  the 
method  of  caustics  [29].  Later  Theocaris  [30]  also  analyzed  the  problem  of  asymmetrically 
branched  crack  and  found  good  agreement  with  the  stress  intensity  factors  obtained  exper¬ 
imentally  using  the  method  of  caustics  [29].  The  disadvantage  of  using  the  method  used  by 
Theocaris  to  solve  branched  crack  problem  is  that  instead  of  obtaining  one  complex  singular 
integral  equation  as  in  the  case  of  dislocation  method,  two  complex  singular  integral  equa¬ 
tions  are  obtained. 


Mellin  Transform 

Khrapkov  [31],  using  Mellin  transforms,  obtained  the  solution  for  a  branched  crack 
emanating  from  a  semi-infinite  notch.  Bilby  and  Cardew  [32]  obtained  the  stress  intensity 
factors  in  front  of  the  kink  tip  in  terms  of  the  stress  intensity  factors  of  the  main  crack  for 
infinitesimally  short  kink  lengths.  The  stress  intensity  factors  for  the  in-plane  and  anti¬ 
plane  problem  were  calculated  using  the  methodology  developed  by  Khrapkov  [31].  They 
included  the  constant  stress  term  (T-stress)  parallel  to  the  main  crack  in  the  analysis  and 
made  use  of  the  maximum  energy  release  criterion  to  predict  the  initial  angle  of  crack  kink- 
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ing.  Later  Bilby,  et  al.,  [33]  solved  the  problem  of  forked  cracks  using  the  same  methodology 
as  discussed  in  [3 1  ] .  They  compared  their  results  for  a  single  branched  crack  with  Chatterjee 
[18]  and  Kitagawa,  etal.,  [21].  Their  results  for  a  branched  crack  agreed  with  Chatterjee  [18] 
but  disagreed  with  Kitagawa,  et  al.  [21].  Melin  [34]  used  the  Mellin  transform  technique  to 
obtain  the  stress  intensity  factors  accurately  in  front  of  an  infinitesimally  small  branched 
crack. 

Perturbation  Technique 

Banichuk  [35],  using  perturbation  techniques  and  Muskhelishvili  complex  stress 
functions  [15],  obtained  the  solution  for  a  weakly  curved,  infinitesimally  small  crack,  in  a 
brittle  material  using  the  local  symmetry  of  the  stress  state  (i.e.,  Kii=0)  and  the  Griffith  ener¬ 
gy  criterion.  Gol’dstein  and  Salganik  [36]  used  perturbation  methods  to  calculate  the  stress 
intensity  factors  in  front  of  an  infinitesimally  short  kink.  They  used  the  criterion  of  local 
symmetry  (i.e.,  Kii=0)  to  determine  the  growth  direction  of  the  kink.  Cotterell  and  Rice  [7] 
used  methods  similar  to  the  ones  used  by  Banichuk  [35]  and  Gol’dstein,  et  al.  [36],  to  solve 
the  infinitesimally  small  branched  crack  with  a  small  deviation  from  the  main  crack  direc¬ 
tion.  They  also  incorporated  the  constant  stress  acting  parallel  to  the  main  crack  (T-stress) 
in  their  formulation  and  determined  a  stability  criterion  for  crack  branching  using  the  T- 
stress.  For  an  infinitesimally  small  kink  Cotterell  and  Rice  [7]  compared  their  calculated 
stress  intensity  factors  with  those  obtained  by  Bilby,  et  al.,  [33]  and  found  good  agreement 
even  for  larger  (up  to  40  degrees)  angles  of  kink.  For  finite  kinks  their  results  matched  well 
with  those  of  Kitagawa,  et  al.,  [21].  They  also  expressed  the  stress  intensity  factors  in  front 
of  the  kink  in  terms  of  the  stress  intensity  factor  of  the  main  crack  and  coefficients  which  are 
functions  of  the  the  kink  angle.  These  coefficients  are  the  terms  containing  the  angular  varia¬ 
tion  in  the  Irwin-Williams  asymptotic  stress  expansion  in  front  of  the  crack  tip.  The  stress 
intensity  factors  obtained  by  Cotterell  and  Rice  [7]  were  of  first  order  approximation.  Sumi, 
et  al.,  [37]  using  perturbation  methods  and  an  alternating  scheme  developed  a  solution  to  de- 
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termine  crack  growth  path  in  a  finite  geometry.  They  assumed  a  shape  for  the  crack  extension 
and  by  using  the  principle  of  local  symmetry  determined  the  crack  growth  path  in  a  finite 
geometry.  They  also  noted  that  for  a  first  order  approximation,  the  far  field  stresses  in  a  finite 
body  were  not  affected  by  the  shape  of  the  crack  extension.  Karihaloo,  et  ah,  [38]  made  use 
of  the  perturbation  technique  to  obtain  the  curvature  of  the  infinitesimally  small  growing 
kink.  They  reported  that  the  curvature  depended  not  only  on  the  T-stress  term  but  also  on 
the  derivatives  of  the  stress  intensities  of  the  main  crack  with  respect  to  the  main  crack  length. 
Leblond  [39]  and  Amestoy  and  Leblond  [40]  obtained  a  more  exact  expression  for  the  stress 
intensity  factors  in  front  of  the  kink  by  using  more  terms  in  their  asymptotic  expansion.  By 
making  use  of  the  principle  of  local  symmetry  and  maximum  energy  release  rate  criterion 
they  proposed  a  method  to  trace  the  path  of  a  growing  mixed-mode  crack. 

Literature  Review  On  the  Branched  Crack  Problem-Anisotropic  Material 
Branched  cracks  in  an  anisotropic  body  have  been  studied  considerably  less  than  the 
isotropic  case.  Eshelby  [4 1  ]  was  the  first  to  give  expressions  for  the  elastic  displacement  and 
energy  due  to  a  dislocation  in  an  anisotropic  body.  This  work  was  followed  by  Frank  [42] 
who  laid  down  the  foundations  of  crystal  dislocation  theory.  Eshelby,  et  al.,  [43]  derived  the 
general  solution  of  the  elasticity  equations  for  an  arbitrary  homogeneous  solid  where  the 
elastic  state  is  independent  of  one  of  the  three  Cartesian  coordinates.  Following  this  earlier 
work  Stroh  [44]  developed  the  expressions  for  the  stresses  in  an  general  anisotropic  body 
(where  the  field  is  independent  of  one  of  the  three  Cartesian  coordinates)  due  to  a  dislocation, 
a  pile  of  parallel  dislocations  and  a  crack.  Atkinson  [45]  investigated  the  behavior  of  a  mov¬ 
ing  semi-infinite  plane  crack  in  an  anisotropic  material  where  the  direction  of  propagation 
lies  in  the  plane  of  material  symmetry.  Atkinson  [46]  followed  this  work  with  a  study  of  the 
interaction  between  a  crack  and  a  dislocation  in  an  infinite  anisotropic  body  where  again  the 


field  is  independent  of  one  of  the  coordinate  axis. 

A  general  anisotropic  elasticity  theory  was  also  developed  independently  by  Lekhnits- 


kii  [47],  Green  and  Zema  [48]  and  Savin  [49].  The  foundations  of  fracture  mechanics  as  ap¬ 
plied  to  cracks  in  anisotropic  bodies  was  first  introduced  by  Sih,  et  al.,  [50]  and  Wu  [51]. 
Steen  Krenk  [52]  developed  the  general  solution  for  a  finite  number  of  collinear  cracks  in 
an  anisotropic  infinite  body  using  the  Muskhelishvili  complex  variable  formulation.  From 
the  analysis  he  concluded  that  if  the  loads  on  the  crack  surface  were  self  equilibrating  then 
the  normal  and  shear  stresses  on  the  lines  of  the  crack  were  independent  of  the  elastic 
constants.  He  also  developed  expressions  for  stress  intensity  factors  in  front  of  the  crack  tip. 
Tan  and  Bigelow  [53]  used  the  boundary  force  method  to  analyze  an  edge  crack  and  an  in¬ 
clined  single  edge  crack  in  a  finite  anisotropic  body. 

Obata,  et  al.,  [54]  solved  the  branched  crack  problem  using  dislocation  techniques. 
Their  method  was  similar  to  that  developed  by  Lo  [32]  for  an  isotropic  case.  Their  analysis 
was  valid  only  for  infinitesimally  small  branched  crack.  They  calculated  the  stress  intensity 
factors  and  energy  release  rate  in  front  of  the  kink  and  from  the  analysis  concluded  that  the 
behavior  of  the  stress  intensity  factor  was  similar  to  that  of  the  isotropic  case  whereas  the 
behavior  of  the  energy  release  rate  was  different.  Their  conclusion  was  based  on  analyzing 
the  problem  in  which  the  stiffer  material  axis  was  always  perpendicular  to  or  inclined  at  a 
small  angle  to  the  main  crack.  Obata,  et  al.,  [54]  did  not  consider  anisotropy  in  the  fracture 
toughness  of  the  material. 

Chiang  [55]  gave  an  approximate  solution  for  determining  the  stress  intensity  factors 
at  the  kink  tip  in  terms  of  the  stress  intensity  factors  of  the  main  crack  for  the  case  where  the 
kink  is  infinitesimally  small.  The  asymptotic  stress  field  around  the  kink  tip  was  derived  for 
the  degenerate  case  when  the  anisotropy  becomes  isotropic.  Chiang  [55]  also  notes  that  the 
maximum  stress  intensity  necessarily  does  not  occur  when  the  kink  angle  is  zero. 

Gao  and  Chiu  [56]  used  Stroh’s  [44]  formulation  to  solve  the  infinitesimally  small 
branched  crack  problem  in  an  anisotropic  body.  With  their  model  they  were  able  to  analyze 
wavy,  jogged  and  kinked  cracks.  They  provided  first  and  second  order  perturbation  solutions 
for  stress  intensity  factors  at  the  tip  of  the  kink  in  terms  of  the  stress  intensity  factors  of  the 
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main  or  macro  crack.  They  also  included  the  T-stress  acting  parallel  to  the  macro  crack  in 
their  formulation.  It  was  observed  that  the  stress  intensity  based  fracture  criteria  leads  to 
physically  unreasonable  branching  angles  whereas  the  strain  energy  based  fracture  criteria 
gave  physically  reasonable  results  for  crack  branching  angles.  Based  on  the  energy  criterion 
they  were  able  to  show  that  the  T-stress  can  be  used  in  the  same  way  as  in  the  isotropic  case 
to  determine  the  stability  of  crack  propagation.  They  also  did  not  include  a  variation  in 
strength  in  their  formulation. 

Xu  and  Keer  [57]  used  a  perturbation  technique  in  a  manner  similar  to  Cotterell  and 
Rice  [33]  and  Sumi,  et  al.,  [35]  to  express  the  stress  intensity  factors  in  front  of  the  kink  in 
terms  of  the  stress  intensity  factors  of  the  main  crack.  They  incorporated  anisotropy  in 
strength  in  their  model.  By  using  the  normal  stress  ratio  criterion  developed  by  Buczek  and 
Herakovich  [58]  and  Gregory  and  Herakovich  [59]  they  were  able  to  show  that  the  T-stress 
may  not  behave  the  same  way  as  in  the  isotropic  case.  For  example,  they  showed  that  for 
some  stiffness  ratio  and  material  orientation  a  negative  T-stress  caused  unstable  crack 
growth  which  is  contrary  to  what  is  observed  in  isotropic  materials. 

Since  the  present  study  involves  development  of  a  fracture  criteria  to  determine  the 
behavior  of  a  kinked  crack,  a  brief  survey  of  the  development  of  several  fracture  criteria  is 
given  below. 

Review  of  Salient  Fracture  Criteria-Isotropic  Material 

In  the  next  few  paragraphs  some  of  the  important  fracture  criteria  are  outlined  for  an 
isotropic  material.  The  fracture  criteria  can  be  generally  classified  as  either  stress  based  or 
energy  based.  All  the  criteria  discussed  below  are  only  for  plane  problems. 

Griffith  [60,61]  was  the  first  to  develop  the  energy  based  fracture  criteria.  He  postu¬ 
lated  that  the  crack  in  a  material  grew  if  the  rate  of  change  of  elastic  energy  of  a  material  with 
respect  to  the  crack  length  exceeded  a  critical  energy  rate.  He  determined  that  every  material 
had  a  specific  critical  energy  rate,  called  the  critical  energy  release  rate,  which  if  exceeded 


caused  the  crack  to  grow. 

Irwin  [62]  used  the  stress  intensity  factors  to  determine  when  the  crack  grew  and 
showed  the  equivalence  between  the  energy  release  rate  criteria  and  the  stress  intensity  fac¬ 
tors.  All  the  work  for  some  time  was  related  to  mode  one  fracture. 

Erdogan  and  Sih  [29]  and  Panasyuk,  et  al.,  [63]  were  the  first  to  study  the  mixed  mode 
static  crack  problem.  Erdogan  and  Sih  [29]  argued  that  the  crack  travels  in  a  direction  per¬ 
pendicular  to  the  maximum  hoop  stress  which  also  corresponds  to  zero  shear  stress.  Based 
on  experiments  they  performed,  they  observed  that  their  model  predicted  experimental  re¬ 
sults  reasonably  well  even  though  some  discrepancy  existed.  Williams  and  Ewing  [64]  stu¬ 
died  this  problem  further  and  showed  that  by  including  the  non-singular  stress  term,  namely 
the  T-stress,  the  experimental  results  obtained  by  Erdogan  and  Sih  [29]  could  be  matched 
much  better  with  their  analytical  results.  Since  the  T-stress  was  included,  the  hoop  stress 
had  to  be  calculated  at  a  certain  distance  from  the  crack  tip.  This  distance  was  considered 
to  be  a  material  parameter  which  varied  with  materials  but  was  fixed  for  a  particular  material. 
Finnie  and  Saith  [65]  pointed  out  an  error  in  the  Williams  and  Ewing  [64]  calculation  which 
substantially  improved  the  correlation  between  experiment  and  theory.  Sih  [66-67] 
introduced  a  fracture  criterion  based  on  strain  energy  density  called  the  S-criterion.  He  pro¬ 
posed  that  fracture  occurs  along  the  direction  in  which  the  strain  energy  density  attains  a 
minimum  value  at  a  specific  distance  from  the  crack  tip.  This  distance  is  a  material  parame¬ 
ter.  A  sample  of  the  literature  dealing  with  the  energy  or  the  stress  based  fracture  criteria  is 
listed  in  reference  [68-72  ]. 


Review  of  Fracture  Criteria-Anisotropic  Materials 
The  fracture  criteria  for  anisotropic  bodies  are  much  less  developed  than  those  for  iso¬ 
tropic  bodies.  In  fact,  it  has  been  only  recently  that  a  number  of  task  groups  have  been  estab¬ 
lished  by  ASTM  to  determine  if  isotropic  fracture  criterion  could  be  extended  to  orthotropic 
bodies.  Albritton  and  Goree  [5]  carried  out  experiments  on  silicon  carbide  reinforced  alumi- 
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num  and  determined  that  the  ASTM  standards  used  for  isotropic  material  cannot  be  applied 
to  an  orthotropic  material  like  silicon  carbide  reinforced  aluminum. 

Sih,  et  al.,  [50]  developed  the  equations  for  the  near-crack-tip  asymptotic  stress  fields 
for  both  in-plane  and  anti-plane  problems.  They  also  developed  methods  for  determining 
stress  intensity  factors  at  the  crack  tip  in  an  anisotropic  body.  The  stress  intensities  were  de¬ 
fined  in  such  a  manner  that  the  form  of  the  stress  intensity  factor  was  similar  to  the  isotropic 
case  but  the  actual  magnitudes  of  the  stresses  were  different  since  they  depended  on  the  mate¬ 
rial  parameters.  Based  on  this,  Sih,  et  al.,  [50]  concluded  that  same  fracture  criteria  devel¬ 
oped  for  isotropic  materials  could  be  extended  to  the  anisotropic  materials.  They  did  not  con¬ 
sider  the  variation  in  strength  for  the  anisotropic  material. 

Maiti  [73]  used  the  stress  intensities  for  an  infinitesimal  small  kink  to  derive  the  energy 
release  rate  for  an  infinitesimal  small  kink  and  studied  the  variation  of  the  energy  release  rate 
with  the  kink  angle. 

Saouma,  et  al.,  [74]  studied  the  inclined  crack  problem  (mixed-mode)  in  an  anisotrop¬ 
ic  body.  They  maximized  the  ratio  of  the  circumferential  stress  to  the  critical  tensile  strength 
of  the  material  and  observed  that  the  material  anisostropy  significantly  affected  crack  exten¬ 
sion  when  compared  to  the  predictions  made  by  a  similar  isotropic  model. 

Zhang  [75]  studied  the  problem  of  relating  the  fracture  toughness,  parallel  and  perpen¬ 
dicular  to  the  fiber  of  an  orthotropic  material.  He  proposed  including  the  shear  stress  term 
for  deriving  the  stress  state  in  front  of  a  crack  perpendicular  to  a  fiber.  By  including  the  shear 
stress  term,  the  theoretical  prediction  matched  closely  to  the  experimental  results. 

Sarkar  and  Maiti  [76]  provided  a  formula  for  relating  the  critical  fracture  toughness 
of  an  orthotropic  material  due  to  a  crack  parallel  to  the  fibers,  to  the  critical  fracture  tough¬ 
ness  when  the  crack  is  inclined  at  an  angle  to  the  fibers. 

Zhang,  et  al.,  [77]  proposed  a  new  energy  based  fracture  criterion  to  predict  mixed 
mode  crack  initiation  and  propagation  called  the  Z-criterion.  This  criterion  suggests  that 
mode  I  crack  initiation  is  controlled  by  dilational  strain  energy,  mode  II  crack  initiation  is 
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controlled  by  distortional  strain  energy  and  mixed  mode  crack  initiation  and  propagation  by 
the  damage  energy  density  factor  called  the  Z-criterion.  Zhang,  et  al.,  [78]  extended  Tsai- 
Hill  and  Norris  criteria  to  predict  crack  turning  in  orthotropic  materials. 

Obata,  et  al.,  [54]  solved  the  problem  of  an  infinitesimally  small  kinked  crack  in  an 
orthotropic  material.  They  base  their  formulation  on  Lekhnitskii’s  [47]  approach.  Based  on 
their  results,  they  concluded  that  the  behavior  of  the  stress  intensity  factors  for  an  orthotropic 
material  is  similar  in  nature  to  isotropic  materials.  But  they  note  that  the  variation  in  energy 
is  different  from  the  isotropic  case. 

Chiang  [55]  solved  the  problem  of  a  semi-infinite  crack  with  a  kink  in  an  orthotropic 
material.  He  expressed  the  stress  intensity  factors  in  front  of  the  kink  in  terms  of  the  the  stress 
intensity  of  the  main  crack  and  also  derived  the  asymptotic  stress  field  around  the  crack  tip 
when  the  material  becomes  isotropic. 

Ellyin  and  El  Kadi  [79]  proposed  a  critical  strain  energy  density  criterion  that  predicts 
crack  growth  direction  closer  to  experimental  results.  They  also  compare  their  results  with 
the  normal  stress  ratio  criterion  and  minimum  strain  energy  density  and  find  that  their  model 
predicts  crack  turning  angle  better  than  the  above  mentioned  models. 

Gao  and  Chiu  [56]  use  Stroh’s  formulation  to  solve  the  kinked  crack  problem  in  an  ani¬ 
sotropic  body.  They  point  out  to  some  erroneous  conclusions  by  Obata,  et  al  [54]  which  are 
described  below.  Gao  and  Chiu  [56]  observe  that  when  the  main  crack  is  parallel  to  the  fiber 
direction  the  local  maximum  for  Ki  does  not  occur  when  the  kink  angle  is  zero  which  is  con¬ 
trary  to  what  was  reported  by  Obata,  et  al.,  [54].  They  show  that  using  stress  intensity  factors 
to  predict  crack  turning  angles  leads  to  some  physically  impossible  conclusions  which  are 
contrary  to  the  isotropic  case.  The  energy  based  criteria  is  shown  to  be  a  much  better  way 
of  calculating  the  kink  angles.  This  point  was  missed  by  Obata,  et  al.,  [54]  since  the  formula 
they  used  to  calculate  the  energy  release  rate  was  incorrect.  Gao,  et  al.,  [56]  also  show  that, 
based  on  energy  criteria,  the  Cotterell  and  Rice  [7]  T-stress  criteria  to  predict  crack  growth 
direction  in  isotropic  material  can  be  also  extended  to  orthotropic  materials.  This  conclusion 
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should  be  treated  with  caution  since  they  did  not  include  a  variation  in  fracture  toughness 
in  their  model. 

Xu  and  Keer  [57]  used  a  perturbation  technique  to  obtain  the  solution  of  an  infinitesi¬ 
mally  small  kinked  crack  in  an  orthotropic  material.  They  use  the  normal  stress  ratio  criteri¬ 
on  to  predict  crack  growth  direction  and  show  that,  for  some  stiffness  ratios,  negative  T- 
stress  leads  to  unstable  crack  growth  which  is  contrary  to  the  isotropic  case. 

This  brief  survey  of  existing  literature  suggests  that  much  works  needs  to  be  done  in 
the  area  of  fracture  mechanics  as  related  to  orthotropic  materials. 

How  is  the  Present  Study  Different  from  Other  T-stress  Based  Studies? 

In  the  review  of  literature  so  far  discussed,  for  isotropic  and  anisotropic  materials,  it 
should  be  noted  that  the  T-stress  calculated  in  front  of  the  main  crack  tip  has  been  used  only 
to  determine  the  stability  of  the  crack  growth.  In  other  words  the  constant  term  in  the  asymp¬ 
totic  stress  expansion  in  front  of  a  flat  crack  has  been  used  to  determine  the  stability  of  the 
crack  once  it  branches  or  kinks.  What  this  study  attempts  to  do  differently  is  to  first  calculate 
the  local  T-stress  in  front  of  the  branched  crack  or  kink  rather  than  calculate  the  T-stress 
when  the  kink  is  flat  (or  coincides  with  the  main  crack).  By  doing  this,  the  stability  of  the 
branched  or  kinked  crack  itself  can  be  discussed  in  terms  of  the  local  T-stress,  which  is  a 
function  of  the  length  and  angle  of  the  kink  and  the  remote  loads.  Also  Cotterell  and  Rice’s 
[7]  T-stress  criteria  is  modified  to  include  an  experimentally  determined  Tcrit  value  and  the 
modified  criteria  is  applied  to  the  branched  crack,  and  the  behavior  of  the  branched  crack  is 
discussed  in  light  of  the  modified  T-stress  criteria. 

An  Overview  of  this  Study 

This  thesis  is  divided  into  two  major  subdivisions;  (1)  The  Isotropic  Case  and  (2)  The 
Anisotropic  Case. 


The  Isotropic  Case:  First  a  description  of  the  problem  is  given  followed  by  the  for- 
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mulation  of  the  analytical  solution,  which  is  similar  to  the  one  derived  by  Lo  [24] .  The  differ¬ 
ence  between  this  formulation  and  Lo  [24]  is  that  a  remote  biaxial  loading  is  included  in  the 
model.  The  analytical  solution  is  obtained  in  terms  of  a  singular  integral  equation  in  which 
the  slope  of  the  kink-opening-displacement  is  the  unknown.  This  is  followed  by  a  proof 
showing  that  the  singularity  at  the  internal  comer  is  given  by  the  roots  of  the  Williams’  char¬ 
acteristic  equation  [19].  Then  the  singular  integral  equation  is  solved  for  the  unknown  slopes 
of  the  crack  opening  displacement  using  three  different  numerical  schemes.  The  advantages 
and  disadvantages  of  each  numerical  scheme  is  discussed  in  detail  based  on  computational 
time,  numerical  stability  and  rate  of  convergence.  Next,  the  method  used  to  calculate  the 
stress  intensity  factors  and  the  local  T-stress  in  front  of  the  kink  tip  is  outlined.  This  is  fol¬ 
lowed  by  the  description  of  the  T-stress  based  fracture  criteria  to  predict  how  an  already 
kinked  crack  behaves.  The  influence  of  biaxial  load  on  the  local  T-stress  is  also  studied. 

The  Anisotropic  Case:  The  details  in  this  subdivision  are  very  similar  to  that  outlined 
in  the  isotropic  Case.  A  brief  description  of  the  problem  is  followed  by  a  general  analytical 
solution  to  a  branched  crack  problem  in  an  ortho  tropic  body.  This  is  an  extension  of  the 
methodology  established  by  Obata,  et  al.,  [54].  The  Obata,  et  al.,  [54]  solution  can  be  used 
only  for  infinitesimal  small  kink  lengths  where  as  the  present  formulation  is  capable  of  han¬ 
dling  any  kink  length.  The  solution  is  obtained  in  terms  of  a  singular  integral  equation  where 
the  slope  of  the  crack  opening  displacement  is  the  unknown.  The  singular  integral  equation 
is  solved  using  three  different  numerical  schemes.  Again  the  effectiveness  of  each  numerical 
scheme  is  discussed  in  the  light  of  computational  time,  numerical  stability  and  rate  of  conver¬ 
gence.  The  methods  for  calculating  the  stress  intensity  factors  and  the  T-stress  are  estab¬ 
lished.  The  influence  of  anisotropy  on  the  T-stress  and  stress  intensity  factors  are  discussed 
and  a  T-stress  based  fracture  criteria  is  proposed  which  helps  in  the  better  understanding  of 
crack  turning  in  an  anisotropic  material. 


CHAPTER  II 


ISOTROPIC  MATERIAL 
Analytical  Formulation 

The  formulation  of  the  branched  crack  problem  in  an  isotropic  material  follows  the 
procedure  first  developed  by  Lo  [24].  Lo’s  formulation  is  extended  to  include  bi-axial  loads. 
Given  below  are  the  necessary  equations  in  the  final  form.  For  more  details  on  the  derivation 
of  the  equations,  the  reader  is  referred  to  Appendix  A. 

The  geometry  of  the  problem  with  the  relevant  labels  is  shown  in  Figure  1  (page  2). 
Throughout  the  development  of  the  analytical  solution  for  the  above  branched  crack  prob¬ 
lem,  Muskhelishvili’s  [15]  complex  variable  approach  is  used.  Dislocations  are  used  to  mod¬ 
el  the  kinked  crack.  A  brief  description  of  the  dislocation  function  is  given  next. 

Description  of  Dislocation  Function 

Physically,  on  an  atomic  level,  dislocations  occurs  because  of  some  vacant  spaces  in 
an  atomic  structure.  In  other  words  the  continuity  in  the  matter  is  lost.  A  line  dislocation 
can  be  defined  as  a  boundary  between  a  perfect  (unslipped)  and  imperfect  (slipped)  region 

in  a  crystal.  The  direction  and  magnitude  of  the  slip  is  given  by  the  Burgers  vector  b.  This 
line  dislocation  is  classified  into  two  types,  namely  the  edge  dislocation  and  the  screw  dis¬ 
location,  depending  on  the  direction  of  the  Burgers  vector.  If  the  Burgers  vector  is  perpendic¬ 
ular  to  the  dislocation  line  (AD)  it  is  known  as  an  edge  dislocation  (Figure  2)  and  if  the  Burg¬ 
ers  vector  is  parallel  to  the  dislocation  it  is  known  as  a  screw  dislocation  (Figure  3). 

Mathematically  an  edge  dislocation  is  used  to  represent  inplane  shear  (Figure  4)  and 
normal  stress  (Figure  5)  acting  on  a  crack  surface,  while  the  screw  dislocation  is  used  to  rep¬ 
resent  antiplane  shear  (Figure  6). 
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From  a  macroscopic  view  point  a  crack  can  be  assumed  to  be  made  up  of  a  series  of 
these  dislocations.  Therefore,  the  stress  fields  due  to  a  single  dislocation  can  be  viewed 
mathematically  as  a  Green’s  function  for  a  point  discontinuity;  and  to  obtain  the  stress  field 
for  a  crack  of  finite  length  I,  the  Green’s  function  is  integrated  over  the  length  /. 

Method  of  Superposition  Explained 

Superposition  is  used  to  solve  the  branched  crack  problem  (Figure  7).  First  the  stress 
functions  (OD(z),\itD(z))  due  to  a  single  edge  dislocation  in  an  infinite  isotropic  body  is  for¬ 
mulated,  this  is  step  1  in  Figure  7.  The  stresses  due  to  this  line  dislocation  on  the  main  crack 
which  is  of  length  2c  is  obtained.  The  negative  of  this  stress  is  applied  on  the  main  crack 
which  is  step  2  in  Figure  7.  The  stress  potentials  (OM(z),’V|rM(z))  due  to  this  loading  condi¬ 
tions  are  next  developed.  By  summing  up  step  1  and  2,  the  stress  free  main  crack  in  the  pres¬ 
ence  of  a  single  dislocation  is  obtained.  To  obtain  the  stress  functions  for  a  kink  the  sum  of 
the  stress  functions  obtained  in  step  1  and  2  are  integrated  over  the  length  of  the  kink.  Next 
the  stress  fields  (<I>L(z),tltL(z))  due  to  the  main  crack  alone  with  loads  at  infinity  are  obtained 
as  step  3  in  Figure  7.  The  stress  potentials  for  the  kinked  crack  problem  is  given  by  the  addi¬ 
tion  of  the  above  three  solutions.  Finally  a  singular  integral  equation,  in  terms  of  the  disloca¬ 
tion  density  function  ©(to),  is  obtained  by  satisfying  the  stress  free  boundary  condition  on 
the  kink.  The  unknown  ©(to)  is  solved  for  by  using  a  suitable  numerical  scheme. 

Formulation  of  the  Problem 
Step  1 

The  stress  functions  due  to  a  single  dislocation  in  an  infinite  plane  can  be  written  as 
(see  equations  A.1-A.16). 

(1)d(z)  =  0(to)ln(z  -  Zo)  (1) 

tl)D(z)  =  0(to)ln(z  -  Zo)  -  0(to)zo/(z  -  Zo).  (2) 
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where  z  defines  the  position  of  a  point  in  the  infinite  plane  and  Zq  defines  the  position  of  a 
dislocation.  0(to)  is  the  unknown  function  which  is  given  by  the  jumps  in  displacements 
across  the  crack  surfaces. 

0(to)  =  Eb(to)e‘®/<ra(>t  +1)}  •  (3) 

b(to)  =  [Ur"^  -  Ur“]  +  i[u^  -  U^].  (4) 

where,  z=te‘®^,  Zo=toe^®,  E=Youngs  Modulus,  b(to)  is  the  Burgers’  vector,  0  is  the  kink  angle 
and  K  is  a  constant  whose  value  depends  on  the  plane  stress  or  plane  strain  assumption.  The 

components  of  the  Burgers’  vector  b(to)  is  shown  in  Figure  8  and  Figure  9. 


Figure  8.  Tangential  Jump 


The  stresses  are  related  to  the  differentials  of  the  above  functions  ((t)^,  ■V|/'d)  and  the  dif¬ 
ferentials  can  be  written  as. 


(j)'D(z)  =  ^>d(z)  =  ©(to)/(z  -  Zo)- 


(5) 
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Step  2 

The  stresses  acting  on  the  main  crack  due  to  this  dislocation  is  next  evaluated  and  the 
negative  of  these  stresses  are  applied  on  the  main  crack  as  shown  in  Figure  7.  The  corre¬ 
sponding  stress  functions  are  given  as  (see  A.20-A.72), 

<1)m(z)  =  -  [a(to)F(z,  Zo)  +  a(to)F(z,  Zo)  +  a(to)(zo  -  Zo)G(z,  Zo)]>  (7) 

Wm(z)  =  $m(z)  - 

where, 


F(z,Zo)  =  I 


1  - 


X(Zo) 

X(z) 


/(z  -  Zo)» 


G(z,  Zo)  =  Zo)' 


Step  3 

The  stress  functions  for  step  3,  i.e.  the  main  crack  with  remote  infinite  loading  in  both 
the  X  and  Y  directions  can  be  written  down  as  (see  A.73-A.85), 


o 


Ol(z)  = 


z  -I-  c  ^  ^ 


1 


X(z) 


a 


W.Cz)  =  ^ 


yy 


zc^ 

(X(z))- 


+  (1  -  k) 


(9) 


(10) 


where, 


X(z)  =  yz(z  +  2c),  k  =  ^  (11) 

Oyy 

In  the  above  equations,  c  is  the  half  length  of  the  main  crack,  and  k  is  the  biaxial  stress 
ratio.  By  the  principle  of  superposition  the  stress  functions  developed  in  steps  1,2  and  3  can 
be  added  to  obtain  the  complete  solution  of  the  problem,  as  given  below. 

(12) 


<I)(z)  =  Od(z)  +  <I)m(z)  +  Ol(z)- 
W(z)  =  Wd(z)  +  'Fm(z)  +  'Fl(z). 


(13) 
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Once  the  stress  functions  are  known,  the  stresses  at  any  point  in  the  infinite  plane  can 
be  written  in  terms  of  the  stress  functions  in  polar  co-ordinates  as  follows. 

Oqq(z)  +  iOje(z)  =  <I>(z)  +  <I>(z)  +  e^‘®[z<l>'(z)  +  'F(z)]*  (14) 

an-(z)  -  iOj^(z)  =  <h(z)  +  <I)(z)  -  e^*®[z<l>'(z)  +  W(z)]-  (15) 

On  the  kink  surface  both  the  normal  stress  (099)  and  the  shear  stress  (ar9)  vanish.  Im¬ 
posing  this  condition  and  using  Plemelj  formula  [  1 5]  the  following  singular  integral  equation 
is  obtained  (see  A.86-A.90). 

i  1 

OQQit)  +  ia^(t)  =  +  j  M(t,to;0(to))dto  +  ae”e(z)  +  ia^(z)  =  0, 

0  0 

(16) 

where, 

M(t,  to,  a(to))  =  <>m(z)  +  ^m(z)  +  e2'®[(z  -  z)<I>m'(z)  -  ^^(z)  +  ^(z)]-  (17) 

The  above  singular  integral  equation  is  in  terms  of  the  unknown  0(to)  whose  behavior 
at  the  end  points  0  and  I  has  to  be  determined.  From  the  geometry  of  the  crack  problem 
(Figure  1 ,  page  2)  it  can  be  observed  that  the  tip  of  the  kink  (A)  has  the  well  known  square 
root  singularity.  At  the  other  end  of  the  kink  (i.e.  at  the  origin),  the  upper  side  of  the  kink, 
indicated  by  the  letter  U,  is  an  external  comer  of  a  wedge  and  according  to  Williams  [19], 
does  not  give  rise  to  any  stress  singularity.  The  lower  side  of  the  kink,  represented  by  C,  is 
an  internal  comer  of  a  wedge  and  (following  Williams  [19])  one  would  expect  it  to  give  rise 
to  singular  stresses.  It  has  been  shown  by  Chatterjee  [1 8],  by  using  conformal  mapping  tech¬ 
niques,  that  the  lower  side  of  the  kink  represented  by  C  indeed  gives  rise  to  singular  stresses. 
A  similar  exercise  using  the  dislocation  technique  is  lacking  in  the  literature.  It  has  been  cus¬ 
tomary  for  researchers  working  the  branched  crack  problem  with  the  dislocation  technique 
to  recognize  that  the  Williams’  singularity  exists  at  the  wedge  internal  comer  without  proof. 
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In  the  section  below,  a  complete  proof  verifying  the  existence  of  the  Williams’  singularity 
at  the  internal  comer  is  given. 


Verification  of  Existence  of  Williams’  Singularity  at  the  Internal  Wedge  Comer 
Extension  of  Theorems  for  Singular  Integrals  to  Bounded  Integrals 
For  developing  the  proofs  to  verify  the  Williams’  singularity  at  the  lower  side  of  the 
kink,  it  is  necessary  to  extend  theorems  developed  for  singular  integrals  to  bounded  integrals. 
The  bounded  integral  near  the  integration  limits  tends  to  a  bounded  (finite  non-zero)  part 
and  a  part  which  vanishes.  It  is  essential  to  know  how  the  vanishing  term  of  the  bounded 
integral  goes  to  zero. 

Following  Muskhelishvili  [83] ,  consider  the  behavior  of  a  Cauchy  integral  of  the  form, 

i 


F(z)  = 


«(to) 
to  -  Z 


dto’ 


(18) 


0 

near  the  limits  of  integration.  The  function  a(to)  is  bounded  every  where  in  the  region  of 
integration  (0  - 1)  (i.e.  it  satisfies  Holders  condition  in  the  range  0  -  /)  except  near  the  end 
points  0  and  1.  Near  the  end  points  a(to)  can  be  assumed  to  be  of  the  form 


a(to)  = 


?^(to) 

t?(/  -  to)P 


(19) 


where  X,(to)  is  bounded  in  the  entire  range  of  integration,  a  and  P  are  assumed  to  be  always 
between  0  and  1  in  order  to  have  bounded  integrals.  If  the  above  conditions  are  satisfied, 
then,  near  the  end  of  integration  limit  (say  0),  F(z)  can  be  written  as  follows  (Muskhelishvili 
[83],  see  eq.  29.5,  pg  74). 
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limF(z)  =  lim 

z->0  z^O 


X(to) 

t“(l  -  to)^ 


dtp 

to  -  z 


sin(jra)  z-*o 


+  bf> 


(20) 


where  bf  is  a  bounded  function.  A  similar  relation  ean  be  written  for  the  behavior  of  the  func¬ 
tion  at  the  other  end  1.  If  the  restriction  that  a  >  0  is  removed  for  the  lower  end  (i.e.  a  <1 
only)  then  a  can  have  negative  values  and  hence  in  the  limit  the  integral  has  a  part  that  va¬ 
nishes  and  a  bounded  part.  It  should  be  noted  that  the  same  proof  used  to  get  the  above  rela¬ 
tionship  (20)  can  be  used  even  when  a  is  less  than  0.  Hence  any  asymptotic  value  for  a 

Cauchy  integral  can  be  extended  for  the  range  -  “  =  a  <  i.  n  should  be  noted  that  when 
a  =  0  logarithmic  singularities  occur  and  equation  (20)  is  no  longer  valid. 


Proof  For  Verification  of  Williams  Singularity 
In  order  to  verify  the  singularity  at  the  wedge  comer  it  is  sufficient  to  examine  the  sum 
of  the  hoop  and  radial  stresses.  This  simplifies  the  expression  used  for  the  asymptotic  analy¬ 
sis.  Using  the  stress  functions  the  sum  of  the  radial  and  hoop  stresses  can  be  written  as, 

arr(z)  +  09e(z)  =  4Re[<I>(z)]-  (21) 

Since  all  the  stresses  vanish  at  the  external  comer  of  the  wedge  the  sum  of  the  radial  and  hoop 
stresses  must  also  vanish  at  the  external  comer.  Therefore  using  (5),  (7),  (9)  and  (12)  in  (21) 
the  following  expression  is  obtained. 

Re{0(z)}  =  Re[l,(z)  +  l^(z)  +  ^z)  +  Uz)  +  l5(z)  -f  l^(z)  +  l7(z)  +  Cj]  =  0  (22) 
where, 


i 

=  i  f  Q(to) 

2  J  z  -  Zp^^o’ 


Il(z)  =  T 


T  _  II  ®(to)  t. 

-  2  z-Zp^^°’ 


I 
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1,(7)  = 


1  ©(to)X(Zo) 


2  X(z)(z  -  Zo) 


I4(Z)  =  ^ 


1  0(to)X(Zo) 


2  X(z)(z  -  Zo) 


l5(^)  =  -  t 


1  0(to)(Zo  -  Zo) 


(z  -  Zo)2 


I6(Z)  = 


1  I  Q(to)(Zo  -  Zo)X(Zo)^ 

2  X(z)  (z  -  Zo)^ 


1  0(to)(Zo  -  Zo)(Zo  +  c) 


2  X(z)X(zo)  (z  -  Zo)' 


r  =  -  i  + 

'  4  2X(z)  ’ 

where,  X(z)  is  given  by  (1 1). 

Using  the  limiting  value  of  each  integral  near  the  end  point,  i.e.  z  =  0,  the  asymptotic 
values  of  each  can  be  written  down.  Integrals  with  X(z)  in  their  denominator  require  special 
attention.  It  is  these  integrals  which  have  to  be  evaluated  using  the  behavior  for  bounded 
integrals  as  discussed  earlier.  For  example,  consider  integral  I3  which  can  be  written  as. 


1  0(to)X(Zo) 


^  2  X(z)(z  -  Zo) 


1  0(to)X(Zo) 

2X(z)  (z-zo)  ' 


Now  (see  B.5,  B.6), 
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0(to)X(zo) 
(z  -  Zo) 


dto  = 


jr  P(0)ie  “  “  ™)]e  -  i[0(i  -  «)] 

cos(3ta)  , 


+  bfj^ 


therefore, 


i;™  T  _  i;„  1  ©(to)X(Zo) 
lim  13(2)  lim  ^  ry  \ 

t^o  t~*o2  X(z)(z  “  Zo) 


3rP(O)ie-‘f“(0a-’^)]e-‘Wi-«)l  1  bfj 

= - ;;; - 7 r -  +  — p-  (z4; 

2cos(jia)  ft 

where  bfs  is  a  bounded  function.  If  a  is  assumed  to  be  less  than  1/2  then  the  first  term  in 
equation  (23)  vanishes.  If  the  vanishing  term  is  neglected  this  would  leave  only  the  bounded 
part  which  from  equation  (24)  results  in  only  the  bfa/t^^^  jf  j^g  vanishing  term  is 

retained  and  it  is  multiplied  by  X(z)  then  the  first  term  in  equation  (24)  results,  plus  the  term 
containing  bfs.  The  other  integrals  can  be  evaluated  in  a  similar  manner.  Refer  to  Appendix 
B  for  more  details. 


r  T/N  1-  1  I  ©(to)  j  jtP(O)e-‘f«(9a-Jt)]e-'[0(i-a)]  j 

limlj(z)  =  lim:^  v  -  /  dtp  = - - - r-r — r -  hm;^  +  bf,-  (25) 

t-o  ^  t-*o2  z  Zo  o  2sm(jta)  ^ 


T  /  X  V  1  ©(to)  j  jtp(O)e->[a(0a-Ji)]e-i[e(a-i)]  1 

hml2(z)  =  hm:^  - ^dto  = - r-p — r -  +  bfo-  (2^) 

t^o  ^  t^o2  z  -  Zo  2sm(jta)  t-^ot'*  ^ 


i:„  T  _  i:„  1  ©(to)X(Zo) 
linil4(z)  lini—  .  irrcito 

t^o  ^  t-o2  X(z)(z  -  Zo) 


n:p(O)ie“iW®a-ii;)]e-i[0(a-i)]  j  ^  bf4 
2cos(jta)  !™t^  ~J‘ 


limljfz)  =  -  limi  ^dto 

t—O  t-*o2  (z  -  Zo)^ 


0 
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+  {Aj  +  iA2}y  =  0  (32) 

where,  Bf={bfi+bf2+bf5-l/4},  Ai  =  Re{bf3+bf4+bf6+bf7+bfc}  and  A2=Im{bf3+bf4+ 
bf6+bf7+bfc}. 

As  t  tends  to  zero  the  coefficient  of  the  1/t^^^  term  and  the  coefficient  of  l/t“  term 
should  vanish.  Considering  the  coefficient  of  the  l/t“  term  it  can  be  written  as, 

Re  P(O)e“''^“*^®'“’^^~?^sin[0(a  —  1)  +  Jca] 


+  P(0)sin(9)e”‘f“^®"  2)(l  -  a)e  2)+jia]  _  q.  (33) 

Since  if  the  real  part  of  a  complex  function  is  zero  the  real  part  of  the  complex  conju¬ 
gate  should  also  be  zero.  Therefore, 

Rej^|3(O)e‘f“*'®''“’’^^“5^sin[0(a  —  1)  -I-  Jta] 

-1-  .p(O)sin(0)e't“(®»-’'>-t\l  -  =  q.  (34) 


The  above  equation  (34)  can  be  written  in  matrix  form  as  follows, 


rsin[0(a  -  1)  +  Jta]e 


-i(0(a-2)) 


sin(0)(l  -  a)e  ^’^“e 


-ijta_i[a(6a-Ji)-f)'! 


ji(0(a-2))gjn(0)(i  _  cc)e4a(0a-Jt)-f)ei^a  sin[0(a  -  1)  + 


If  p(0)  and  P(0j  do  not  vanish,  then  the  determinant  of  the  above  equation  should  van¬ 
ish  (i.e.  at  least  one  non-trivial  solution  exists).  Therefore, 
sin[0(a  -  1)  +  Jta]  +  (1  -  a)sin0  =  0 

sin[0(a  -  1)  +  Jta]  -  (1  -  a)sin0  =  0 

If  equations  (36)  or  (37)  or  both  are  satisfied  by  choosing  a  suitably,  then  the  coeffi¬ 
cient  of  l/t“  term  vanishes.  Considering  the  1/t^''^  terms,  an  equation  of  the  form. 


Ai  =  0> 


is  obtained  as  t  tends  to  zero. 


Now  let  z  approach  the  internal  comer  (C  in  Figure  1)  of  the  wedge.  X(z)  is  the  only 
function  which  changes  sign  as  the  internal  comer  is  approached.  Performing  an  asymptotic 
analysis  as  before,  the  following  relationship  is  obtained. 


Oniz)  +  afin(z)  =  4Re 


P(O)e“4a(9a-ii:)-f]sin[0(l  —  a)  +  Jta] 
sin(2jta) 


3(0)e~i[“(®»“^)“f>e“^(®(“  2))sin(0)(i  -  a)e  1 
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-  {A,  +  iAj}^  (39) 

The  coefficient  of  the  l/t“  terms  in  the  above  equation  is  not  satisfied  if  equations  (36) 
or  (37)  is  satisfied.  From  equation  (38)  the  coefficient  of  the  t^^^  term  vanishes.  Therefore 
as  the  reentrant  comer  is  approached  the  stresses  are  l/t“  singular. 

Consider  equations  (36)  and  (37)  which  are  the  same  expression  obtained  for  an  iso¬ 
tropic  wedge  by  Timoshenko[84].  According  to  Timoshenko  [84]  the  Williams  characteris¬ 
tic  equations  for  a  wedge  can  be  written  as, 

k  sin  2a  +  sin2ka  =  0>  (^9) 

k  sin  2a  -  sin2ka  =  0-  (^1) 

The  geometry  of  the  wedge  is  given  below  in  Figure  10. 


Figure  10.  Geometry  of  the  Wedge  as  Defined  in  Timoshenko,  et  al.,  [84] 


Using  the  following  change  of  variables, 


k  =  1  —  a » 


(42) 
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a=IL^,  (43) 

and  substituting  them  (equations  (42)  and  (43))  in  equations  (40)  and  (41),  equations  (36) 
and  (37)  can  be  recovered. 

It  has  been  conclusively  shown  that,  in  order  for  the  stresses  to  vanish  at  the  external 
comer  of  the  wedge,  the  Williams’  characteristic  equations  for  the  interior  isotropic  wedge 
has  to  be  satisfied.  Since  the  order  of  the  stress  singularity  at  the  kink  comer  (C  and  U  in 
Figure  1,  page  2)  is  determined,  the  singular  integral  equation  (16)  can  now  be  solved  by 
using  a  suitable  numerical  scheme. 


CHAPTER  m 


ANISOTROPIC  MATERIAL 
Analytical  Formulation 

The  formulation  of  the  branched  crack  problem  in  an  anisotropic  material  follows  the 
procedure  developed  by  Krenk[52]  and  Lo  [24].  This  problem  was  also  solved  by  Obata, 
et  al.,  [54]  using  a  different  procedure,  but  the  solution  they  developed  was  restricted  to  the 
infinitesimally  small  kink.  In  contrast,  the  solution  developed  below  can  be  used  for  any 
kink  length.  The  necessary  equations  in  the  final  form  are  given  below.  For  more  details 
on  the  derivation  of  the  equations  the  reader  is  referred  to  Appendix  C. 

The  geometry  of  the  problem  with  the  relevant  labels  is  shown  in  Figure  1 1 .  An  infi¬ 
nite  anisotropic  body  with  a  main  crack  of  length  2c  and  a  kink  (or  branched  crack)  of  length 
I  is  shown.  The  branched  crack  makes  an  angle  0  with  the  X  axis.  Bi-axial  loads  a  and  ko 
are  applied  remotely  along  the  Y  and  X  axis  respectively. 

Throughout  the  development  of  the  analytical  solution  for  the  above  branched  crack 
problem,  Muskhelishvili’s  [14]  complex  variable  approach  is  used.  The  problem  is  solved 
using  the  method  of  superposition  in  the  same  manner  as  for  the  isotropic  case  (Chapter  II). 


Development  of  Dislocation  Function 
The  dislocation  stress  functions  are  given  as  (see  C.40-C.55), 

^  SjBx  By  2 

d(  l)  _  S2)(Sj  —  S3)(Sj  —  S4)Cjj  (Zj  —  Zqj) 

826^  By  I 

2)  ^  d(  2)  ^g^  _  Sj)(S2  ~  S3)(S2  —  S4)Cjj  (Z2  ~  ZQ2) 

where  Bx  and  By  are  the  components  of  the  Burgers’  vector. 


(44) 

(45) 
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Z;  =  X  +  Sjy,  Zoi  =  Xo  +  SiYo-  (46) 

where,  i  =  1,  4.  Xo  and  jo  indicates  the  position  of  the  dislocation  in  the  infinite  plane. 

Formulation  of  Stress  Functions  for  the  Main  Crack  Which  Opens  Up  Due  to  the 

Presence  of  a  Dislocation 


The  stress  fields  on  the  main  crack  due  to  a  single  dislocation  are  next  evaluated  by 


using  the  dislocation  stress  functions.  Then  the  negative  of  these  stresses  are  applied  on  the 
main  crack  and  the  stress  functions  are  developed,  following  the  formulation  of  BCrenk  [52]. 


After  satisfying  the  condition  for  single  valuedness  of  displacement,  the  following  expres¬ 
sions  are  obtained  (see  C.56-C.117). 


^m(z)  = 


1  1  1 
Cii  4jri(s2  -  Sj) 


SjBx  By 
(Sj  -  SjXSj  -  S4) 


X(Zoi)\ 

X(z)  ) 


_ S3BX  By _ 1  /  ^  _  X(Zq3)\ 

(S3  -  S4)(S3  -  Si)  Z-Zo3\'  X(z)  ) 

S4BX  -  By _ 1  (  .  X(Zq4)\ 

(S4  -  S3)(S4  -  Si)  z-Z04\  X(z)  )  ■■ 

vu  ('7'!  =  J _ 1  1  ^zBx  -  By _ 1  / .  _  X(Zq2)\ 

^  Cl  1  4ji;i  (Sj  -  S2)  1  (S2  -  S3)(S2  -  S4)  Z  -  Z02  X(z)  j 


S3BX  By _ 1  (  ^  X(Zq3)\ 

(S3  -  S4)(S3  -  S2)  Z-Z03(^'  X(Z)  ) 

S4BX  -  By _ 1  /  .  _  X(Zo4)\  1 

(S4  -  S3)(S4  -  S2)  Z-Z04|^'  X(Z)  jj' 


Formulation  of  Stress  Functions  for  the  Main  Crack  Which  Opens  Up  Due  to  Loads  at 

Infinity 

The  stress  functions  for  the  problem  of  the  main  crack  with  remote  loadings  can  be 


obtained  using  the  same  general  approach  as  discussed  for  the  isotropic  case  and  are  given 
below  (seeC.118-C.133). 
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Ol(z) 

Wl(z)  =  2Re 


(  ^2<^yy  z  ,  .Ox”  +  (aia2  -  PlP2)Oy“ 

2  Re  - r^^TTT  +  i" 


2(S2-Si)X(z)  2(S2  -  SiXPi  +  P2) 


SlOyy 


z  .  Ox”  +  (0102  -  PlP2)Oy” 


2(Si  -S2)X(z)  ■  ^  2(Si  -  S2)(Pi  +  P2) 


(47) 


(48) 


Formulation  of  the  Singular  Integral  Equation 
By  satisfying  the  stress-free  condition  on  the  branched  crack,  the  following  expres¬ 
sions  are  obtained. 


Oee  - 


■f 


t  -  tr 


dto 


i 

-I-  I  2Re 
0 


2  2 
(sjSinB  -t-  COS0)  0]^(zj)  -f  (s2sin0  -I-  cos0)  W}^(z2) 


dto 


+  2  Re 


2  2 

(sjsin0  -I-  COS0)  Ol(zi)  +  (s2sin0  +  cos0)  'Rl(z2) 


=  0 


(49) 


T  M21®x(to)  M22By(to) 
=-h - ^ ^ - dtr 

0 


t  -  to 


I 

+  1 2Re 
0 


^(1  —  Sjjsm0cos0  —  Sj  cos20jd>j^(zj) 


-I-  ((1  —  S2)sin0cos0  -  S2Cos20jWjy[(z2) 


dto 


+  2  Re 


|l  -  Sjjsin0cos0  —  Sj cos20|<I>l(zj) 


-f  |(l  —  s|)sin0cos0  -  S2Cos20jWL(z2)  =  0 


(50) 


where, 


Z;  =  a  -I-  t(cos  0  -f  Sj  sin  0), 
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Zqj  =  a  +  to(cos0  +  SiSinS), 

and  i  =  1,  2.  My,  M12,  M21  and  M22  are  given  in  the  Appendix  C  (C.136-C.139). 


CHAPTER  IV 

EVALUATION  OF  STRESS  INTENSITY  FACTORS  AND  T-STRESS 

Isotropic  Material 

The  mode  I  and  mode  11  stress  intensity  factors  are  evaluated  by  extracting  the  singular 
parts  of  the  stresses  at  the  tip  of  the  kink.  Only  the  singular  parts  of  the  dislocation  function 
contribute  to  the  stress  singularity  at  the  kink  tip.  So  only  the  stress  functions  associated  with 
the  dislocation  functions  are  examined.  The  stress  intensity  factors  at  the  kink  tip  are  defined 
in  terms  of  hoop  and  shear  stresses  in  front  of  the  kink  as  follows, 

Kj  +  iKy  =  y^(r  -  /)5[a0e  +  ia^^]  (51) 

where  r  and  0  are  defined  in  Figure  12.  The  above  stresses  have  already  been  normalized 
with  respect  to  the  applied  uniaxial  load  Oyy°°. 


Figure  1 2.  Configuration  Describing  Near-Kink-Tip  Stresses  for  Evaluating  Stress  Intensity 
Factors 
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Using  the  stress  functions,  the  stresses  are  written  as, 

Oge  +  io^g  =  Od(z)  +  Od(z)  +  e^*®[z<I)i3'(z)  +  W^Cz)]. 

where. 


(52) 


z  =  re‘®,  Zo  =  toe'®  (0  <  to  <  0-  (5^) 

Substituting  for  the  stress  functions  using  equations  (5)  and  (6)  and  simplifying,  the  follow¬ 
ing  equation  is  obtained. 

/ 


J  r  -  to 


,  where 


P(to) 


and 


(54) 


(55) 


a  and  P  are  the  stress  singularities  at  the  kink  comer  and  the  kink  tip  respectively.  The  above 
equation  is  non-dimensionalized  using  the  following  transformation  equations. 


to  =  4(^0  +  1),  r  —  jrih,  +  1),  dto  -  xdlo- 


(56) 


After  substituting  the  transformation  equations  in  equation  (54)  and  simplifying,  the  follow¬ 
ing  equation  is  obtained, 

1 

(57) 


O00  +  la^e  - 


■  »Rr*’  ( 


P(?o)e'® 


d^o 


(1  +  |o)“(l  -  io)P?  - 


The  limiting  value  of  the  above  integral  as  ^  approaches  1  is  evaluated  using  techniques  in 
Muskhelishvili  [83]  for  singular  integrals. 


lim 


1 

I 


P(io)e 


i0 


d|o  ^  i 


?-*!  (1  +  U“(l  -  -  ^0  l-i2«sinjTP(|  -  1)P 


(58) 


-1 


Using  the  above  equation  in  equation  (57)  and  then  substituting  into  equation  (51)  and  sim¬ 
plifying,  the  following  is  obtained  where  P  =  1/2. 

(59) 


Ki  -I-  iKji  - 


43T^/^P(l)e‘® 


jra  + 1/2 
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In  a  similar  manner,  using  the  method  developed  by  Sih,  et  al.,  [50]  the  stress  intensity 
factors  for  the  main  crack  can  also  be  evaluated.  Only  the  stress  functions  due  to  the  main 
crack  with  loads  at  infinity  and  the  stress  functions  due  to  the  interaction  between  the  main 
crack  and  a  dislocation  contribute  to  the  singular  stresses.  Therefore 

^Mi  “  =  -2ij2K  lim(z  -  +  Ol(z)]  (60) 

Z-»Zi 

where  Kmi  and  Kmii  are  mode  I  and  mode  II  stress  intensity  factors  for  the  main  crack  respec¬ 
tively.  In  this  case  zi  =  -2c  and  c+0.  Substituting  for  the  stress  functions  using  equations 
(7),  and  (9)  and  simplifying,  the  following  equation  is  obtained. 


_  r,..  ^  \  \  \  0(to)X(Zo)^,  ,  1  I 

Km,  -  .Km,,  -  -  2.y2„  I  ^  I 


/ 

l2c\ 


0(to)X(Zo) 


i  0(to)(Zo  -  Zo)X(Zo) 

2  J2c  J  (2c  +  Zo)^ 


i 

A.  ,  i  f  0(to)(Zo  -  Zo)(Zo  +  c)  ,  iv/cj 
X(Zo)(2c  +  Zo)  2J2 


All  the  integrals  in  the  above  equation  are  bounded  and  can  be  evaluated  numerically  once 
the  unknown  ©(to)  is  calculated.  If  the  kink  vanishes  then  0(to)  =  0,  therefore 

^Ml  “  i^Mll  “ 

which  is  the  stress  intensity  of  the  main  crack  alone  under  uniaxial  tensile  load.  When  the 
kink  becomes  flat  i.e.  0  =  0  then  the  mode  I  stress  intensity  factor  Kim  the  tip  of  the  main 
crack  becomes  equal  to  the  mode  I  stress  intensity  factor  at  the  kink  tip  (Ki),  i.e. 


Kmi  ~ 


Evaluation  of  the  T-Stress  for  Isotropic  Material 
The  T-stress  is  defined  as  the  constant,  ‘p  ’  independent  term  in  the  asymptotic  expan¬ 
sion  of  the  stress  field  around  the  tip  of  the  crack  (see  equations  63  -  65).  In  the  case  of  the 


branched  crack  problem  two  types  of  T-stress  can  be  defined.  One  is  the  T-stress  associated 
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with  the  main  crack  alone,  i.e.  without  the  branched  crack  and  the  other  is  the  local  T-stress 
defined  at  the  tip  of  the  branched  crack.  For  a  coordinate  system  as  shown  in  Figure  13,  the 
asymptotic  stress  expansion  is  given  up  to  the  constant  term  in  the  following  equations. 


Figure  13.  Near  Crack-Tip  Stresses  for  the  Branched  Crack. 


Oqq 

_  Ki 

Jlnq 

[|cos| 

1  (b 

-  ;|cOS3^]  + 

Kn  [. 

5  .  4)  3  4)  o 

•^sin^  +  ^sm3y]  +  Tcos^4>’ 

(63) 

_  Ki 

JlnQ 

[|cos| 

3  .  4)  3  4)  T 

-  ^sm^  -  •|sm3^]  +  Tsin^4^» 

(64) 

_  Ki 
JllXQ 

[ism| 

1  4^ 

+  ^sin3  j]  + 

Kii  .1 
y2jT0  4 

4)  3  4) 

cos-j  +  ^cos3-2]  —  Tsin4)cos4>- 

(65) 

lf(|)  =  0 

and  p  =0,  then 

Oqq  -  T- 

(66) 

i.e.  the  T-stress  is  given  by  the  difference  in  the  hoop  stress  and  radial  stress  evaluated  at  the 
kink  tip. 


From  Muskhelishvili  [15], 

-  age  =  2Re([z<D'(z)  +  W(z)]e2ie),  (67) 

where. 
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O  =  0)0  +  ^>M  + 

^  +  'Pm  + 

Consider  the  contribution  from  the  dislocation  function  only,  which  can  be  written  as 


2Re{[zOD'(z)  +  WD(z)]e2‘®)  =  2Re< 

^  1 

[  eie0(to)  -  e-'00(to),. 

r  -  to  ' 

J 

lo  J 

(68) 


where  z  =  re'®.  The  integrand  is  a  pure  imaginary  number  and  therefore  there  is  no  contribu¬ 
tion  to  the  T-stress  from  the  dislocation  function.  Next,  consider  the  contribution  from  0)m 
and 

2Re([zO>M'(z)  +  =  2Re(^(z)  -  O^Cz)  +  (z  -  z)(Dm'(z)1-  (69) 

After  substituting  for  the  above  stress  functions,  by  observation,  the  following  singular  inte¬ 
grals  are  isolated  as  ‘r’  approaches  7’.  The  remaining  integrals  in  the  above  equation  (69) 
are  bounded  and  converge  to  a  finite  limit  which  contributes  to  the  T-stress. 


i 


‘Pm(^)  “  2 


Q(to)  ,  i 
z  —  Zq  2 


0 


=  I 


-_®(i2L_dt  -  i 

(Z  -  Zo)"  2 


©(to)X(Zo)  dtp 

X(z)  z  -  Zc 


0(to)X(Zo)(z  +  c)  dtp 

X(z)3  2 


(70) 


-if  ®(^o' 

2  X 


,)X(Zo)  dto 


^(^)  (z  -  ZoY 


(71) 


i 


^m(2)  “  ^  I  z-  °Zo  2 


0(to)X(zo)  dtp 

X(z)  z  -  Zc 


i  I 

1  f  0(to)(Zo  -  Zo)  ,  1  f  0(to)(Zo  -  Zo)X(Zo: 

2}  °  2j  XW 


)  dtp 

(z  -  Zof 
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1  I  0(to)(zo  -  Zo)(z  +  c)  dtp 

2  X(z)X(Zo)  z  -  z 


Therefore,  substituting  the  above  stress  functions  in  equation  (69)  and  simplifying,  the  fol¬ 
lowing  expression  is  obtained. 

2Re{^(z)  -  <1>m(z)  +  (z  -  z)<1)m'(z)}  = 

i  i _ 

if  ©(to)  ,  1  f  0(to)X(zo)  dto 

-  2  +  2  X(z)  z 


1  f  ©(to)  1  ©(to)X(Zo)  dto 

2  z  -  Zo^^°  2  X(z)  z  -  Zc 


1  j  0(to)(Zo  -  Zo)  ,  1  0(to)(Zo  -  Zo)X(Zo)  dto _ 

2j  (z-zo)^  2j  X(z) 


J_  J  0(to)(zo  -  Zo)(z  +  c)  dtp  _  1_  0(to)X(Zo)(z  +  c)  dtp 
2  X(z)X(Zo)  z  -  Zo  2  X(z)3  z  -  Zc 


1  ©(to) 


©(to)X(Zo) 

dto 

X(z) 

1 

N 

o 

V 

0  0  J 

In  the  limit  as  z  tends  to  I  the  singular  part  of  each  pair  of  integrals  given  above  cancel  each 
other.  Therefore,  the  contribution  to  the  T-stress  from  Om  and  Wm  is  only  a  constant.  Hence, 
if  the  real  part  of  equation  (69)  is  calculated  at  the  crack  tip  it  will  be  a  constant  which  contrib¬ 
utes  to  the  T-stress.  Consider  the  contribution  from  and  to  the  T-stress  which  is  giv¬ 
en  by  equation  (74)  for  a  biaxial  loading  condition. 


2Re([za)L'(z)  +  WL(z)]eZi0}  =  2Re 


c2(z  -  z)  (k-l)V2i0.„- 
2X(z)3  i 
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where  k  is  the  biaxial  stress  ratio.  The  above  equation  can  also  be  evaluated  at  the  kink  tip 
(i.e.  z=/ei®).  Therefore,  by  adding  equations  (74)  and  (69),  the  local  T-stress  can  be  eva¬ 
luated  at  the  kink  tip.  Some  simplifications  to  the  above  expression  for  the  T-stress  can  be 
made  if  an  infinitesimally  small  kink  length  is  assumed.  Under  this  assumption  the  contribu¬ 
tion  to  the  T-stress  from  Om  and  becomes  small.  Consider  the  contribution  from  Ol 
and  Wl  only  which,  for  infinitesimally  small  kink  length  can  be  expressed  as, 

r  /xsinOsin,  „  „ 

2Rej[z<I)L'(z)  -f  WL(z)]e2>e}  =  ^  -  -  (k  -  1)  cos  20  Voyy  (75) 

Therefore, 

T  =  .  -  (k  -  l)cos2e  -Oy^y  +  6  (Km.'J’m)' 

where  5(<1>m*^m)  represents  a  small  but  finite  contribution  to  the  T-stress  from  Om  and  'Rm- 
Using  further  manipulations  the  local  T-stress  in  front  of  the  incipient  kink  can  be  expressed 
in  terms  of  the  main  crack  mode  I  stress  intensity  factor  Kmi  (i.e.  Kmi  is  calculated  in  front 
of  the  main  crack  as  if  the  kink  did  not  exist)  as  follows. 

T  =  cosf^)  —  cos(^)  —  cos29a^(k  —  1)  4-  6(0]^,Wjy[).  (77) 

V2;  \2/j 

From  the  above  equation,  for  a  uniaxial  load  (i.e.  k  =  0),  it  can  be  seen  that  when  0  = 
0,  from  equation  (69),  6(Om»^m)  =  0  and  T  =  — ayy°°.  Which  means  that  when  the  kink  is 
flat  i.e.  when  then  the  axis  of  the  kink  coincides  with  the  axis  of  the  main  crack  the  T-stress 
is  the  negative  of  the  applied  remote  stress.  Now  as  0  is  varied  the  T— stress  becomes  positive 
rapidly  because  I  is  very  small.  Note  that  the  above  expression  for  the  T-stress  can  be  used 
only  for  an  infinitesimally  small  kink  length.  Equation  77  can  also  be  expressed  as  follows. 


-  (aH  -  a“)  -h  8(a>M,^M)- 


(78) 
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where  0||°°  and  OnT|°°  are  the  components  of  the  remote  loads  parallel  to  and  perpendicular 
to  the  kink  respectively  as  shown  in  Figure  14. 


Figure  14.  Components  of  Remote  Loads  Parallel  to  and  Perpendicular  to  the  Kink. 


0 
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M2|Bx(to)  "I"  M22By(to) 

r  -  to 


dto' 


(82) 


Let, 


Bx(to)  ~ 
By(to) 


Px(to) 

tg(/  -  tof 

Py(to) 


(83) 


(84) 


t?(/  -  to)P 

Consider  the  integral  associated  with  099.  Using  the  following  transformations,  the  integral 
equation  can  be  normalized  in  the  range  from  -1  to  1 . 


to  =  |(lo  +  1),  r  =  |(i  +  1),  dto  =  |dio- 

r  M,iMo)  +  Mi2py(lo) 


-1 


(1  +  ^o)“(l  -  lo)P(?  -  lo) 


d^o- 


(85) 

(86) 


As  ^  tends  to  1  the  singular  part  is  given  as, 

M„px(l)  +  Mi2Py(l) 


^06  “ 


(87) 


/a  +  p  (I  _  i)p 

Substituting  the  above  equation  in  (79)  and  simplifying,  after  letting  P  =1/2,  the  following 
is  obtained. 

^  2j^3/2  MitPx(l)  +  M;2Py(l)  j-gg^ 

In  a  similar  manner. 


X.  _  ._3/2  M2iPx(1)  +  M22py(l) 
-  2Jt 


(89) 


Evaluating  the  T-stress  for  Anisotropic  Material 
Although  the  method  of  calculating  the  stress  intensity  factors  are  the  same  for  the  iso¬ 
tropic  and  anisotropic  cases,  the  method  of  calculating  the  T-stress  is  different  for  anisotrop¬ 
ic  material  since  the  equations  do  not  lend  themselves  to  simplifications  as  in  the  isotropic 
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case.  The  T-stress  is  calculated  by  evaluating  the  stresses  on  the  flanks  of  the  kink  or 
branched  crack.  The  only  stress  which  is  non-zero  on  the  flank  of  the  branched  crack  is  the 
radial  stress.  The  radial  stress  on  the  flank,  near  the  kink  tip,  contains  singular  terms  and  a 
constant  term  which  is  the  T-stress.  By  isolating  the  singular  term  the  constant  term  can  be 
calculated  more  accurately.  The  radial  stress  for  any  kink  orientation  9  is  given  as. 


On-  =  2  Re  <J>(Zi)(si  COS0  —  sin0)^  +  'P(z2)(s2  cos  0  —  sin0)^  » 


(90) 


where. 


4)(z,)  =  <I.d(z,)  -f  4>m(Zi)  +  «’l(zi)' 

'{'(z,)  =  -PcCz,)  -f  WmCz,)  >Pl(Z2).  (92) 

The  subscript  D  indicates  stress  functions  associated  with  the  dislocation  only,  the  subscript 
M  indicates  stress  function  associated  with  the  interaction  due  to  the  main  crack  and  disloca¬ 
tion  and  subscript  L  is  associated  with  the  main  crack  under  remote  loading.  The  radial  stress 
On-  can  thus  be  split  up  into  three  components, 

Orr  =  og  +  ojf  +  O^,  (93) 

where  the  superscripts  have  their  usual  meaning.  Only  the  dislocation  component  of  the 
stress  (o^it)  gives  rise  to  singular  stresses  as  the  crack  tip  is  approached  along  the  kink.  The 
remaining  components  are  all  bounded  on  the  flank  of  the  kink.  The  radial  stress  due  to  the 
dislocation  is  given  as, 

=  2Re[<I>j5(Zj)(SjCOS0  —  sin0)^  4-  Wj)(z2)(s2  cos  0  —  sin0)^J-  (94) 
Substituting  for  the  stress  functions  and  simplifying. 


= 

Urr  — 


M^3lBx(to)  T  M32By('-o) 


Zl  -  Zol 


dto» 


(95) 


M31 


2  Re 


1  1 
2jtiCll 


(Sicos0  —  sin0)Zsj 

(Sj  -  S2)(Si  -  S3)(Si  -  S4) 


where, 


47 


+ 


(S2COS0  -  sin0)2s2 

(S2  -  Sj)(S2  -  S3)(S2  -  S4) 


M32  =  -  2Re 


J__l_ 

23TiCii 


(Sicos0  -  sin  0)2 

(Si  -  S2)(Si  -  S3)(Si  -  S4) 


+ 


(S2COS0  —  sin  0)2 


(S2  -  Si)(S2  -  S3)(S2  -  S4) 


and  zi  is  any  point  on  the  infinite  plane  except  on  the  kink  or  branched  crack. 

For  calculating  stresses  on  the  flank,  zi  is  allowed  to  be  on  the  line  of  discontinuity 
(i.e.  the  branched  crack).  Therefore,  using  the  Plemelj  formulae, 

i 


=  ±  Jii|M3jBx(t)  +  M32: 


By(t)]  +-j-- 


M3|Bx(to)  +  M32By(to) 


dto' 


(96) 


0 

The  above  relation  is  valid  at  all  points  on  the  line  of  discontinuity  (0  -  1)  except  at  the  end 
points.  But  the  stresses  we  are  interested  in  has  to  be  evaluated  at  the  end  point.  This  is  not 
a  problem  since  accurate  T-stresses  can  be  obtained  by  removing  the  singular  component 
and  then  evaluating  the  stresses  very  close  to  the  end  point  but  not  at  the  end  point  itself.  The 
singular  component  comes  from  the  first  terms  in  equation  (96)  and  the  second  term  is 
bounded  for  all  values  of  t  and  to  except  at  the  end  points.  Therefore,  by  neglecting  the  first 
term  in  equation  (96)  and  evaluating  the  sum  of  the  terms  in  equation  (93)  accurate  T- 
stresses  can  be  calculated  on  the  flank  of  the  kink  since  the  numerical  evaluation  is  no  longer 
obscured  by  terms  which  are  very  large. 


CHAPTER  V 


RESULTS  AND  DISCUSSION-I 
Isotropic  Material 

The  singular  integral  equation  (16)  is  solved  using  the  three  methods  discussed  in  Ap¬ 
pendix  E.  As  discussed  in  Appendix  E,  in  the  first  method  (MI)  a  stronger  singularity  is  as¬ 
sumed  at  the  intersection  of  the  kink  comer  and  the  main  crack,  i.e.  a  square  root  singularity 
is  assumed  at  the  kink  comer  while  in  reality  the  singularity  there  is  less  than  square  root. 
In  the  other  two  methods  (Mil  and  Mill)  the  correct  singularity,  which  is  given  by  the  roots 
of  the  equations  (36)  and  (37)  is  assumed  at  the  kink  comer.  The  numerical  results  for  stress 
intensity  factors  (mode  I  and  mode  11)  given  by  the  three  methods  are  compared  with  the 
solution  obtained  by  Isida,  et  al.,  [86]  for  a  range  of  kink  to  main  crack  length  ratio  (//2c) 
and  kink  angles  (0).  The  boundary  force  method  combined  with  perturbation  techniques  was 
used  by  Isida,  et  al.,  [86]  to  obtain  the  solution  for  the  branched  crack  problem  with  various 
loading  conditions.  The  loading  configurations  for  which  the  comparisons  are  made  are 
shown  in  Figure  15,  Figure  16  and  Figure  17.  The  reason  for  choosing  Isida’s  [86]  among 
the  other  available  solutions  in  the  literature  is  because  numerical  values  can  be  obtained  ac¬ 
curately  from  tables  as  opposed  to  obtaining  numbers  from  a  graph  as  in  other  cases.  The 
numerical  values  compared  are  the  mode  I  and  II  stress  intensity  factors  (namely  Ki  and  Kn) 
which  are  normalized  with  respect  to  the  stress  intensity  factor  of  the  main  crack.  One  typical 
comparison  using  method  I  (MI)  subjected  to  load  case  1  for  a  kink  to  main  crack  length  ratio 
of  0.005  is  shown  in  Table  1 .  Numerical  results  using  methods  11  (Mil)  and  HI  (Mill)  are 
compared  in  Table  2  and  3.  From  the  numerical  values  given  in  these  tables  the  following 
observations  can  be  made.  Method  I,  although  it  assumes  a  stronger  (square  root)  singularity 
at  the  kink  comer,  still  gives  accurate  results.  Method  II  uses  the  correct  singularity  at  the 
kink  comer  but  it  requires  a  large  number  of  terms  for  convergence.  This  is  due  mainly  to 


the  presence  of  other  less  dominant  roots  of  the  Williams  characteristic  equations  (36)  and 
(37)  which  results  in  the  unboundedness  of  the  slope  of  the  bounded  function.  Once  these 
roots  are  taken  care  of,  as  in  Method  III,  the  convergence  rate  is  very  fast  (almost  as  fast  as 
method  I).  The  convergence  is  fast  because  not  only  the  bounded  function  but  also  the  singu¬ 
lar  nature  of  the  slope  of  the  bounded  function  is  taken  care  of  in  the  series  representation 
of  the  unknown  as  discussed  in  the  Appendix  E.  Table  4  shows  the  two  roots  of  the  Williams’ 
characteristic  equation  used  in  the  numerical  method.  The  fourth  column  in  the  table  gives 
the  slope  of  the  bounded  function  and  the  negative  value  is  indicative  of  the  fact  that  the  de¬ 
rivatives  are  unbounded  at  the  kink  comer.  It  is  due  to  this  fact  that  the  method  II  has  a  very 
slow  convergence.  Numerical  results  for  load  cases  2  and  3  shown  in  Figures  16  and  17  are 
compared  with  Isida,  et  al.,  [86]  in  Tables  5  and  6  for  a  kink  to  main  crack  length  ratio  of 
0.005  using  method  I. 

Tables  7  and  8  compares  the  mode  I  and  mode  II  stress  intensity  factors  obtained  by 
using  method  I  and  III  respectively  with  the  solutions  obtained  by  Isida  et  al.,  [86]  and  Mellin 
[34]  when  the  kink  becomes  infinitesimally  small.  All  the  results  given  are  for  a  unit  loading 
at  infinity,  normal  to  the  main  crack.  The  difference  in  the  stress  intensity  factors  between 
the  three  methods  is  always  less  than  1%. 

All  the  three  methods  give  accurate  values  for  mode  I  and  11  stress  intensity  factors. 
Method  II  (not  shown)  and  HI,  where  the  correct  singularity  is  assumed  compare  more  favor¬ 
ably  with  the  result  obtained  by  Mellin  [34]  and  Isida  [86]  for  larger  values  of  kink  angles. 

From  the  above  comparison  of  mode  I  and  II  stress  intensity  factors  it  can  be  concluded 
that  the  nature  of  the  stress  singularity  at  the  kink  comer  does  not  influence  the  stress  intensi¬ 
ty  factors  at  the  kink  tip.  A  similar  conclusion  about  the  T-stress  can  be  made  from  Table  9 
where  it  can  be  observed  that  the  difference  in  the  T-stress  values  calculated  by  the  three 
methods  are  negligible  (<  .5%).  Though  the  particular  case  shown  is  for  a  kink  to  main  crack 
length  ratio  of  0.005  the  results  hold  tme  for  all  range  of  kink  to  main  crack  length  ratios. 
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Table  1.  Load  Case  1  (No.  of  Quadrature  Points  =  121,  Method  I) 


//2c 

0 

Ki[MI] 

Kn[MI] 

KiPsida] 

KnPsi- 

da] 

Diff% 

Ki 

Diff% 

Kn 

0.005 

0 

1.0025 

0.0 

1.0025 

0.0 

0 

0 

0.005 

10 

0.9911 

0.0862 

0.9912 

0.0862 

0.01 

0.0 

0.005 

20 

0.9579 

0.1676 

0.9580 

0.1676 

0.01 

0.0 

0.005 

40 

0.8344 

0.2989 

0.8344 

0.2989 

0.0 

0.0 

0.005 

60 

0.6594 

0.3688 

0.6592 

0.3690 

0.03 

0.05 

0.005 

80 

0.4677 

0.3705 

0.4672 

0.3708 

0.11 

0.08 

0.005 

90 

0.3759 

0.3484 

0.3752 

0.3489 

0.18 

0.14 

Table  2.  Load  Case  1  (No.  of  Terms  in  Polynomial,  n  = 

40,  Method  11) 

//2c 

0 

Ki[Mn] 

Kn[MII] 

KiPsida] 

Kn[Isi- 

da] 

Diff% 

Ki 

Diff% 

Kii 

0.005 

0 

1.0027 

0.0 

1.0025 

0.0 

0.02 

0.0 

0.005 

10 

0.9912 

0.0862 

0.9912 

0.0862 

0.0 

0.0 

0.005 

20 

0.9579 

0.1676 

0.9580 

0.1676 

0.01 

0.0 

0.005 

40 

0.8345 

0.2990 

0.8344 

0.2989 

0.01 

0.03 

0.005 

60 

0.6593 

0.3691 

0.6592 

0.3690 

0.02 

0.03 

0.005 

80 

0.4673 

0.3710 

0.4672 

0.3708 

0.06 

0.03 

0.005 

90 

0.3753 

0.3490 

0.3752 

0.3489 

0.03 

0.03 

Table  3.  Load  Case  1  (No.  of  Terms  in  Polynomial,  ni=  6,  n2  =4,  [see  E.14]  Method  III) 


//2c 

0 

Ki[Mni] 

KnLMIII] 

KiPsida] 

KnPsi- 

da] 

Diff% 

Ki 

Diff% 

Kn 

0.005 

10 

0.9912 

0.0862 

0.9912 

0.0862 

0.0 

0.0 

0.005 

20 

0.9579 

0.1676 

0.9580 

0.1676 

0.01 

0.0 

0.005 

40 

0.8344 

0.2990 

0.8344 

0.2989 

0.0 

0.03 

0.005 

60 

0.6592 

0.3691 

0.6592 

0.3690 

0.0 

0.03 

0.005 

80 

0.4673 

0.3710 

0.4672 

0.3708 

0.06 

0.03 

0.005 

90 

0.3752 

0.3490 

0.3752 

0.3489 

0.0 

0.03 

51 

Table  4.  Two  Roots  of  the  Williams’  Characteristic  Equations  (36)  and  (37)  Used  in  Numeri¬ 
cal  Analysis 


e 

cti 

a2 

Ctl-0t2-l 

10 

0.0999 

0.0 

-0.9001 

20 

0.1813 

0.0 

-0.8187 

40 

0.3028 

0.0 

-0.6972 

60 

0.3842 

0.0 

-0.6158 

80 

0.4371 

0.0195 

-0.5824 

90 

0.4555 

0.0914 

-0.6359 

Table  5.  Load  Case  2  (No.  of  Quadrature  Points  =  121) 


//2c 

0 

Ki[MI] 

Kn[MI] 

Ki[Isida] 

Kn[Isi- 

da] 

Diff% 

Ki 

Diff% 

Kn 

0.005 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.005 

10 

-0.0027 

0.0154 

-0.0027 

0.0154 

0. 

0.0 

0.005 

20 

-0.0105 

0.0291 

-0.0105 

0.0291 

0.0 

0.0 

0.005 

40 

-0.0377 

0.0460 

-0.0376 

0.0459 

0.27 

0.22 

0.005 

60 

-0.0698 

0.0438 

-0.0698 

0.0438 

0.0 

0.0 

0.005 

80 

-0.0937 

0.0246 

-0.0936 

0.0246 

0.11 

0.0 

0.005 

90 

-0.0993 

0.0113 

-0.0993 

0.0114 

0.0 

0.88 

Table  6.  Load  Case  3  (No.  of  Quadrature  Points  =121) 


//2c 

0 

Ki[MI] 

Kn[MI] 

Ki[Isida] 

Kn[Isi- 

da] 

Diff% 

Ki 

Diff% 

Kii 

0.005 

0 

0.0 

1.0025 

0.0 

1.0025 

0.0 

0.0 

0.005 

10 

-0.2607 

0.9785 

-0.2607 

0.9785 

0.0 

0.0 

0.005 

20 

-0.5086 

0.9085 

-0.5086 

0.9083 

0.0 

0.02 

0.005 

40 

-0.9201 

0.6524 

-0.9199 

0.6520 

0.02 

0.06 

0.005 

60 

-1.1659 

0.3035 

-1.1655 

0.3033 

0.03 

0.06 

0.005 

80 

-1.2223 

-0.0495 

-1.2216 

-0.0490 

0.06 

1.02 

0.005 

90 

-1.1848 

-0.2013 

-1.1840 

-0.2004 

0.07 

0.45 
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Table  7.  Comparison  of  Results  with  Mellin  [34]  and  Isida  [86]  for  2c//  —  oo  (Method  I) 


0 

Ki[MI] 

Kn[MI] 

Ki[Isida] 

Kn[Isida] 

Ki[Mellin] 

Kn[Mel- 

lin] 

0 

1.0000 

0.0 

1.0000 

0.0 

1.0000 

0.0 

10 

0.9886 

0.0865 

0.9887 

0.0864 

0.9886 

0.0864 

20 

0.9552 

0.1680 

0.9553 

0.168 

0.9552 

0.1680 

40 

0.8313 

0.2994 

0.8313 

0.2994 

0.8313 

0.2994 

60 

0.6558 

0.3690 

0.6557 

0.3691 

0.6557 

0.3693 

80 

0.4640 

0.3700  n 

0.4637 

0.3703^ 

0.4636 

0.3705 

90 

0.3724 

0.3476 

0.3719 

0.3480 

0.3719 

0.3481 

Table  8.  Comparison  of  Results  with  Mellin  [34]  and  Isida  [86]  for  2c//  =  co  (Method  III) 


e 

Ki[MIII] 

KnlMEI] 

Ki[Isida] 

Kn[Isida] 

Ki[Mellin] 

Kn[Mel- 

lin] 

0 

1.0000 

0.0 

1.0000 

0.0 

1.0000 

0.0 

10 

0.9886 

0.0865 

0.9887 

0.0864 

0.9886 

0.0864 

20 

0.9553 

0.1680 

0.9553 

0.1680 

0.9552 

0.1680 

40 

0.8314 

0.2993 

0.8313 

0.2994 

0.8313 

0.2994 

60 

0.6554 

0.3692 

0.6557 

0.3691 

0.6557 

0.3693 

80 

0.4636 

0.3705 

0.4637 

0.3703 

0.4636 

0.3705 

90 

0.3717 

0.3481 

0.3719 

0.3480 

0.3719 

0.3481 

Table  9.  Comparison  of  T-stress  Obtained  by  the  Three  Methods  for  Unit  Load  at  Infinity 
(//2c=0.005) 


e 

MI 

Mil 

Mm 

0 

-1.0 

-1.0 

-1.0 

10 

0.1069 

0.1069 

0.1072 

20 

0.4185 

0.4185 

0.4197 

40 

1.5329 

1.5327 

1.5369 

60 

2.9709 

2.9692 

2.9769 

80 

4.2550 

4.2497 

4.2600 

90 

4.6975 

4.6900 

4.6992 
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Table  10  shows  the  comparison  of  T-stress  values  evaluated  directly  using  equation 
(76),  after  neglecting  5(Om,Wm),  with  the  values  obtained  by  solving  the  singular  integral 
equation  (16)  and  then  using  equation  (67).  The  difference  between  these  two  values  in¬ 
creases  with  the  kink  angle.  This  is  due  to  the  increase  in  6(<1>m,^m)  iii  equation  (76)  as  the 
kink  angle  increases.  5(<I)m,^m).  which  can  also  be  called  a  correction  term,  can  be  calcu¬ 
lated  for  each  kink  angle  by  subtracting  column  3  from  column  2  in  Table  10.  The  above 
method  of  expressing  the  T-stress  for  vanishingly  small  kink  length  is  similar  to  the  manner 
in  which  Kj  and  Kn  are  expressed  in  terms  of  sums  of  the  main  crack  stress  intensity  factors 
and  a  correction  term  (for  eg.  Karihaloo,  et  al,  [25])  which  is  given  in  a  tabular  form  as  a 
function  of  kink  angles.  From  Table  1 0  it  can  be  inferred  that  equation  (76)  is  valid  or  infini¬ 
tesimally  small  kink  length  and  small  angles.  For  larger  angles  the  correction  term  has  to 
be  used. 


Table  10.  Comparison  of  Eq.  (76)  with  Numerical  T-stress  Eq.  (67)  (2c/Z  =  oo). 


0 

T  (Numerical) 

T  (Eq.(76)) 

%Diff 

0 

-1.00 

-1.00 

0.0 

10 

76.53 

74.73 

2.4 

20 

301.06 

296.12 

1.66 

40 

1101.57 

1099.06 

0.23 

60 

2129.84 

2165.56 

1.65 

80 

3032.29 

3166.05 

4.22 

90 

3346.92 

3536.53 

5.36 

Development  of  Fracture  Model 


A  fracture  model  based  on  the  T-stress  which  will  predict  qualitatively  the  behavior 
of  a  branched  crack  is  discussed  next.  Recall  from  Chapter  IV  that  the  T-stress  is  a  ‘p  ’  inde 
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pendent  constant  term  which  occurs  in  the  asymptotic  expansion  of  stresses  near  the  kink  tip. 
For  a  flat  crack  in  an  infinite  panel  subjected  to  an  uniaxial  normal  load  cyyy°°  the  T— stress 
is  given  as  -ayy°°.  This  T-stress  acts  in  a  direction  parallel  to  the  flat  crack.  Similarly  for 
a  branched  crack  the  T-stress  at  the  kink  tip  can  also  be  described  as  the  stress  which  acts 
parallel  to  the  kink  direction.  For  a  branched  crack,  the  T-stress  is  a  function  of  kink  length, 
angle  and  the  loading  at  infinity  (Oyy°°,  axx°°)  •  This  functional  relationship  can  be  seen  from 
equation  (77)  which  is  for  an  infinitesimally  small  kink.  The  T-stress  can  either  be  tensile 
(positive  T-stress)  or  compressive  (negative  T-stress). 

At  this  stage  it  will  be  useful  to  state  the  Cotterell  and  Rice  [7]  T-stress  based  fracture 
stability  criterion  for  a  flat  crack  (unbranched  crack)  which  will  then  be  used  to  develop  the 
fracture  model  for  a  branched  crack.  Cotterell  and  Rice  [7]  studied  the  influence  of  the  sign 
of  the  T-stress  evaluated  in  front  of  a  flat  crack  on  the  behavior  of  the  flat  crack  when  it  grows 
in  its  own  plane.  They  showed  that  if  the  T-stress  in  front  of  the  kink  was  negative  (or  com¬ 
pressive)  then  the  crack  tends  to  grow  in  the  same  direction  it  started  growing,  i.e.  along  the 
plane  of  the  crack.  This  they  called  class  I  fracture.  If  the  T-stress  was  positive  then  the  crack 
deviates  from  its  own  plane,  this  fracture  they  termed  as  class  II  fracture. 

The  classification  of  the  type  of  fracture  based  on  the  sign  of  the  T-stress  is  based  on 
an  earlier  work  by  Cotterell  [93]  in  which  some  experimental  verification  was  done  regard¬ 
ing  the  influence  of  the  sign  of  the  T-stress  on  crack  growth  direction.  Using  photoelastic 
methods  Cotterell  [93]  showed  that  in  a  plastic  cleavage  type  specimen  a  flat  crack  ran 
straight  when  the  T-stress  in  front  of  the  crack  was  negative  and  the  crack  turned  when  the 
T-stress  was  positive. 

Richardson  [94]  used  the  T-stress  and  the  critical  mode  I  stress  intensity  factor  calcu¬ 
lated  in  front  of  an  edge  crack  in  a  finite  geometry  test  coupon  (isotropic,  PMMA)  to  develop 
a  relationship  between  the  geometry  of  the  coupon  and  critical  mode  I  stress  intensity  factor 
Kic-  This  resulted  in  an  empirical  curve  called  the  Equivalent  Ratio  Biaxial  Stress  (ERBS) 
curve.  The  ERBS  ratio  Be  is  evaluated  by  first  finding  an  equivalent  biaxial  remote  load  for 
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an  infinite  coupon  geometry  with  a  crack  which  would  give  the  same  crack  tip  stress  state 
as  that  of  a  finite  coupon  geometry  with  a  crack  and  then  finding  the  ratio  of  the  biaxial  re¬ 
mote  loads.  This  is  shown  in  Figure  18. 


Figure  1 8.  Comparison  of  an  Arbitrary  Cracked  Coupon  with  an  Infinite  Center-Cracked 
Panel  (Richardson  [94]) 


isotropic  body. 


Be  =  1  + 


TV^ 

Kj 


(98) 
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The  above  equation  (98)  relates  the  T-stress  and  the  mode  I  stress  intensity  factor,  Kj  eva¬ 
luated  in  front  of  the  edge  notched  finite  geometry  coupon  to  the  biaxial  stress  ratio  Be  related 
to  the  infinite  panel.  While  performing  experiments  on  PMMA  test  coupons  of  varying  di¬ 
mensions  Richardson  [94]  observed  that  for  large  positive  values  of  the  ERBS  ratio  Bg  crack 
turning  occurred  and  the  fracture  mode  was  no  longer  purely  mode  I.  Large  values  of  Bg 
results  when  the  T-stress  is  positive  (see  equation  98).  Richardson  [94]  also  observed  that 
crack  turning  did  not  occur  immediately  after  the  T-stress  became  positive  as  stated  by  Cot- 
terell  and  Rice  [7],  but  at  a  considerably  higher  positive  values  of  T-stress.  For  PMMA  a 
critical  T-stress  (Tgnt)  above  which  crack  turning  occurs  can  be  determined  by  using  the 
ERBS  curve  given  in  Richardson  [94]  and  is  found  to  be  0.4489  ksi.  For  materials  other  than 
PMMA  the  critical  T-stress  (Tgnt)  must  be  determined  experimentally.  Therefore  Richard¬ 
son’s  [94]  experiments  shows  a  clear  dependence  of  crack  turning  behavior  on  the  T-stress 
evaluated  in  front  of  the  crack  in  a  finite  geometry  coupon. 

The  T-stress  based  fracture  criterion  developed  for  a  flat  crack  can  be  extended  to  a 
branched  crack.  Consider  the  situation  in  which  a  kink  is  formed  at  an  angle  to  the  main  crack 
as  shown  in  Figure  1 .  If  the  local  T-stress  calculated  in  front  of  the  kink  tip  is  less  than  Tgnt 
then  the  kink  continues  to  run  in  the  same  direction  as  it  started  to  kink.  This  type  of  crack 
growth  is  called  class  I  fracture.  On  the  other  hand  if  the  T-stress  is  greater  than  Tgnt  in  front 
of  the  kink  then  the  kink  will  not  continue  to  run  in  the  direction  of  initial  kinking.  This  type 
of  fracture  in  which  the  kink  turns  from  its  initial  kinking  direction  is  known  as  class  II  frac¬ 
ture.  Hence,  based  on  the  calculated  T-stress  in  front  of  the  kink,  the  direction  of  kink 
growth  can  be  qualitatively  predicted.  The  types  of  fracture  are  illustrated  in  Figure  19. 

In  the  present  model,  the  local  T-stress  in  front  of  the  turned  crack  in  an  infinite  coupon 
can  be  used  to  predict  the  most  likely  direction  the  turned  crack  or  kink  is  going  to  grow. 
Based  on  this  understanding  of  the  behavior  of  the  branched  crack  in  an  infinite  plane  the 
behavior  of  the  turned  crack  in  a  finite  coupon  geometry  can  be  predicted. 
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The  behavior  of  mode  I  and  11  stress  intensity  factors  and  the  T-stress  as  a  function 
of  kink  angle  and  length  for  both  uniaxial  and  biaxial  loading  is  discussed  in  the  next  few 
paragraphs. 

Figures  20  and  21  shows  the  variation  of  the  normalized  mode  I  stress  intensity  factor 
with  the  normalized  kink  length  for  uniaxial  and  biaxial  loading  respectively.  The  biaxial 
stress  ratio  a°°xx/ct°°yy  is  represented  by  k.  Similarly  Figures  22  and  23  represent  the  varia¬ 
tion  in  Kn  and  Figure  24  and  Figure  25  the  variation  in  T-stress. 

Consider  the  uniaxial  loading  case  first  (i.e.  k  =  0,  Figures  20,  22  and  24).  From 
Figure  20  it  can  be  seen  that  the  mode  I  stress  intensity  factor  is  maximum  and  increases  with 
the  kink  length  when  the  kink  angle  is  zero  i.e.  0  =  0.  For  all  other  kink  angles  the  stress 
intensity  factor  is  lower  than  when  0  =  0,  and  initially  drops  with  an  increase  in  kink  length 
before  increasing  again.  Therefore  it  can  be  concluded  that  once  the  crack  turns  it  requires 
more  load  to  extend  the  crack  again  in  the  new  direction.  The  variation  of  mode  II  stress 
intensity  factor  with  kink  length  is  shown  in  Figure  22.  For  all  kink  lengths,  Kn  is  maximum 
when  the  kink  angle  is  between  45  and  60  degrees.  Also  Kn  is  never  zero  for  kink  angles 
other  than  zero.  The  T-stress  (Figure  24)  behaves  in  an  interesting  manner.  The  Tent  value 
for  a  PMMA  material  is  shown  on  the  same  figure.  For  very  short  kink  lengths  (less  than 
0.5%  of  the  main  crack  length)  and  kink  angles  above  zero  the  T-stress  is  greater  than  Tent 
and  as  the  kink  grows  the  T-stress  decreases.  For  kink  angles  less  than  60  degrees  and  kink 
lengths  more  than  40%  of  the  main  crack  length  the  T-stress  is  less  than  Tcrit-  Therefore, 
whenever  a  crack  deviates  from  its  path  to  form  a  kink  such  that  the  kink  length  is  very  small, 
since  the  T-stress  in  front  of  the  kink  is  more  than  Tcrit,  it  is  likely  that  the  kink  will  deviate 
again.  If  the  kink  length  is  long  enough  and  the  angle  of  deviation  is  also  shallow,  since  the 
T-stress  is  less  than  Tent,  the  kink  can  continue  to  run  in  the  same  direction  as  it  first  started 
to  grow.  It  should  be  kept  in  mind  that  the  natural  direction  for  the  crack  to  grow  when  an 
uniaxial  load  is  applied,  and  when  there  is  no  kink,  is  along  its  own  plane  (i.e.  along  Kn  = 


0).  The  crack  might  kink  from  its  plane  of  growth  due  to  material  inhomogeneities  or  due 
to  material  defects. 

The  biaxial  stress  (k  =  5)  state  is  also  examined,  and  the  results  for  mode  I  and  II  stress 
intensity  factors  and  T-stress  are  plotted  in  Figures  21, 23  and  25  respectively.  These  plots 
present  an  interesting  deviation  from  the  case  when  k  =  0.  The  mode  I  stress  intensity  factor 
(Figure  21)  increases  with  an  increasing  kink  length  for  all  kink  angles.  The  mode  II  stress 
intensity  factor  (Figure  23)  shows  a  local  minimum  (i.e.  Kn  =  0)  for  kink  angles  other  than 
zero  degrees.  Therefore  the  crack  can  deviate  and  run  along  directions  other  than  0  =  0.  But 
crack  growth  in  all  directions  is  not  stable  as  seen  in  Figure  25.  Only  for  kink  angles  greater 
than  45  degrees  does  the  T-stress  become  less  than  Tent  and  hence  the  crack  growth  is  stable. 

The  results  discussed  for  the  infinite  coupon  geometry  can  be  qualitatively  extended 
to  the  case  of  finite  coupon  geometry.  Consider  two  tension  specimen  of  varying  dimension 
as  shown  in  Figure  26.  From  Richardson  [94]  it  can  be  noted  that,  on  application  of  load 
and  before  the  crack  starts  growing,  a  compressive  T-stress  exists  in  front  of  the  edge  crack 
in  the  3  in.  tension  specimen  and  a  tensile  T-stress  exist  in  front  of  the  edge  crack  in  the  2 
in.  compact  tension  specimen.  The  behavior  of  the  cracks  in  the  test  coupons  with  compres¬ 
sive  and  tensile  T-stresses  can  be  understood  by  studying  the  branched  crack  in  an  infinite 
coupon  under  uniaxial  and  biaxial  load  respectively. 

First  consider  the  case  of  the  3  in.  tension  specimen  with  T/P  =  -0.6,  where  P  is  the 
applied  load.  Since  only  qualitative  predictions  will  be  made  about  crack  turning  angle  this 
coupon  geometry  can  be  compared  to  a  flat  crack  in  an  infinite  plane  where  T/a  =  -1 ,  where 
o  is  the  applied  load.  It  can  be  seen  from  Figure  24  that  for  6  =  0,  T/a  is  -1.  If  then  due 
to  the  presence  of  some  defect  or  inhomogeneity  in  the  material  the  crack  kinks  out  of  its 
plane,  the  T-stress  in  front  of  the  kink  becomes  greater  than  Tent  for  very  short  kink  lengths 
(less  than  0.5%  of  the  main  crack  length)  and  hence  the  kink  turns  back  towards  the  main 


Figure  26.  Geometry  of  PMMA  Test  Coupons 
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crack  direction  and  continues  to  grow.  The  kink  comes  back  to  the  main  crack  direction  be¬ 
cause  along  that  direction  Kn  =  0  (Figure  22)  and  the  T-stress  is  less  than  Tent  (Figure  24). 
If  the  kink  length  is  long  enough  (more  than  40%  of  the  main  crack  length),  then  according 
to  Figure  24  the  T-stress  becomes  less  than  Tcrit,  and  hence  the  kink  can  grow  at  angles  less 
than  60  degrees.  This  does  not  happen  for  the  3  in.  PMMA  coupons  because  the  kink,  initial¬ 
ly,  does  not  make  long  enough  excursions  from  it  plane  of  growth  to  really  experience  a  de¬ 
crease  in  T-stress  values.  Therefore,  in  the  absence  of  large  inhomogeneities  or  discontinui¬ 
ties  in  the  tension  specimen,  which  might  cause  the  kink  to  have  long  excursion  from  the 
main  crack,  the  crack  in  the  3  in.  tension  coupon  grows  straight  which  is  also  observed  exper¬ 
imentally. 

In  a  similar  manner  the  2  in.  tension  specimen  with  T/P  =  4  can  be  compared  to  an  infi¬ 
nite  coupon  geometry  with  k  =  5.  When  the  kink  angle  is  zero  the  T-stress  is  greater  than 
Tcrit  as  in  the  test  coupon.  Now  if  the  crack  deviates  from  its  plane  of  growth,  for  short  kink 
lengths  (i.e.  for  kink  lengths  less  than  0.5%  of  main  crack),  the  kink  can  either  turn  back  and 
run  along  the  main  crack  direction  or  it  can  deviate  further  away  from  the  main  crack.  It  is 
more  likely  that  the  kink  will  turn  away  from  the  main  crack  direction  and  run  in  a  direction 
along  which  the  T-stress  is  less  than  Tcrit  since  only  along  this  direction  stable  crack  growth 
can  occur.  From  Figure  23  it  can  be  observed  that  the  kink  angles  for  which  Kjj  =  0  (plane 
of  local  symmetry),  is  smaller  for  smaller  kink  lengths.  Therefore,  short  excursions  from 
the  main  crack  direction  results  in  shallow  turning  angles  as  seen  in  PMMA  test  coupons, 
and  longer  excursions  of  the  kink  results  in  sharp  turning  of  kink.  Ultimately  for  the  infinite 
coupon  the  kink  is  likely  to  turn  by  more  than  45  degrees  since  only  then  does  the  T-stress 
become  less  than  Tcrit- 

Kink  Closing  Angle 

Kink  closing  angle  (0c)  can  be  defined  as  that  kink  angle  at  which  the  mode  I  stress 


intensity  factor  (Ki)  vanishes  (becomes  zero).  Figure  27  shows  the  kink  closing  angle  as  a 
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function  of  kink  length  for  uniaxial  loading.  From  this  plot  it  can  be  inferred  that  kink  clo¬ 
sure  occurs  only  for  kink  angles  above  78  degrees  and  kink  lengths  less  than  0.1  times  the 
length  of  the  main  crack 


Figure  27.  Kink  Tip  Closing  Angle  (0c)  as  a  Function  of  Normalized  Kink  Length 


CHAPTER  VI 


RESULTS  AND  DISCUSSION-II 
Anisotropic  Case 

Before  discussing  the  numerical  results  for  the  anisotropic  case  it  will  be  useful  to  de¬ 
fine  some  of  the  parameters  which  are  used  in  the  presentation  of  the  results.  These  parame¬ 
ters  are  shown  in  Figure  28.  X-Y  are  the  geometric  axes  and  1-2  are  the  material  axes  and 
the  angle  between  them,  denoted  by  Gm,  is  the  material  angle.  The  exterior  angle  formed 
between  the  kink  and  the  main  crack  is  denoted  by  623  and  the  angle  between  the  external 
bisector  to  the  wedge  angle  ABC  and  1-axis  is  denoted  by  9b.  6  is  the  angle  between  the 
positive  X  axis  and  the  kink.  The  1 -direction  is  assumed  to  be  the  direction  of  higher  stiff¬ 
ness,  (typically  the  fibre  direction). 

The  singular  integral  equations  (49  and  50)  developed  for  the  anisotropic  case  are 
solved  using  the  three  methods  discussed  in  Appendix  E.  As  discussed  in  Appendix  E,  in 
the  first  method,  a  stronger  singularity  is  assumed  at  the  kink  comer  where  the  kink  and  the 
main  crack  meet.  In  the  other  two  methods  the  correct  singularity,  which  is  given  by  equation 
(99),  is  assumed  at  the  kink  comer.  The  development  of  this  equation  is  given  in  Bogy  [87] 
and  is  repeated  below  for  completeness. 

D  =  -  t5.[2{t4  -  tj}  +  {t2  +  tjtg]  =  0  (99) 

t,  =  [Ri(0a)-^Rl(-  ea)R2(0a)R2(-  63)-']^^+'), 

t2  =  [Ri(0a)Ri(-  0a)-'R2(0a)''R2(-  0a)]^^^'\ 

ts  =  (1  -  YiY2)^cos{(ai  -b  l)[4)i(0a) 

—  4^i(~  0a)  ~  <t^2(0a)  +  ^2^“  0a)]K 

t4  =  (Yi  -  Y2)^cos{(ai  +  l)[([)i(0a) 

-  (j)i(-  0a)  +  (l)2(0a)  -  ^2(-  ^a)]}, 
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t5  =  [Ri(0a)Ri(-  ea)R2(0a)R2(- 
te  =  (1  -  Y?)(l  -  Y2)’ 

where  7i  (i=l,2)  are  the  roots  of  the  following  equation, 

ShY"^  —  4SigY^  +  2(Sgg  +  Sjj)y^  —  4SjgY  +  Sj2  =  0> 

where, 


Sii  =  -^(Cji  +  C22  +  C55  -  2C12)’ 

^12  ”  4(^11  ^22  “  ^66  ~  2Cj2)’ 

^16  “4(^11  “  ^22)’ 

^66  “  4(^11  C22  +  2Ci2)> 

and  Cjj’s  are  the  same  parameters  as  used  in  equation  (C.26).  These  Cij’s  can  be  expressed 
in  terms  of  elasticity  constants  as  shown  below. 

c 

Ej 


"22 


_  1 


_  1^12 


^12  ~ 


^66  ~  ^12' 


Ri  =  (1  +  Yf  +  2YiCOS[2(0b  -  0)])l/2. 

sin0  +  YiSin(20jj  —  0*) 


(j)j(0)  =  atan 


0*  _  Tt  +  0 


COS0  +  YiCos(20j,  —  0*) 


,i  =  1,2- 


In  the  above  equations  Ei,  E2  are  the  stiffness  along  the  1  and  2  material  direction,  G12 


is  the  shear  modulus  and  \i\2  is  the  Poison’s  ratio. 
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Table  1 1  shows  the  roots  of  equation  (99)  for  different  kink  angles  and  material  angles 
for  a  stiffness  ratio  (E1/E2)  of  1.1.  The  numerical  results  for  stress  intensity  factors  (mode 
I  and  mode  II)  and  the  T-stress  are  compared  using  the  three  methods  and  the  results  are 
shown  in  Table  12.  In  method  1, 121  quadrature  points  are  used;  in  method  11, 40  terms  in 
the  Jacobi  polynomial  series;  and  in  Method  III  the  first  Jacobi  polynomial  series  has  6  terms 
and  the  second  Jacobi  polynomial  series  has  4  terms.  From  the  results  shown  in  Table  12 
it  can  be  observed  that  the  error  in  Kj,  Kn  and  T-stress  calculated  using  the  three  methods 
is  always  less  than  1  %  (as  seen  in  the  isotropic  case).  It  can  also  be  observed  that  the  singular¬ 
ity  assumed  at  the  kink  comer  does  not  affect  the  near-kink-tip  stress  state. 


Table  11.  Roots  of  Bogy’s  [87]  Characteristic  Equation  (99) 


0 

02a 

0b 

ai 

Ct2 

10 

95 

85 

0.0196 

0.0 

20 

100 

80 

0.1855 

0.0 

40 

110 

70 

0.3058 

0.0 

60 

120 

60 

0.3850 

0.0 

80 

130 

50 

0.4368 

0.0207 

90 

135 

45 

0.4551 

0.0923 

Table  12.  Kj,  Kn  and  T-Stress  Calculated  Using  the  Three  Methods  (//2c  =  0.005,  E1/E2  = 
1.1) 


MI  (121) 

ME  (n=40) 

Mill  (ni-i-n2  = 

6+4) 

0 

Ki 

Kn 

T 

Ki 

Kn 

T 

Ki 

Kn 

T 

0 

1.0025 

0.0 

-0.038 

1.0024 

0.0 

-0.038 

1.0025 

0.0 

-0.038 

10 

0.9915 

0.0836 

0.0826 

0.9915 

0.0836 

0.0827 

0.9915 

0.0836 

0.0830 

20 

0.9591 

0.1631 

0.4240 

0.9591 

0.1631 

0.4240 

0.9591 

0.1631 

0.4253 

40 

0.8382 

0.2937 

1.5674 

0.8382 

0.2938 

1.5672 

0.8384 

0.2938 

1.5718 

60 

0.6651 

0.3664 

2.9842 

0.6649 

0.3667 

2.9824 

0.6652 

0.3667 

2.9906 

80 

0.4736 

0.3714 

4.2815 

0.4731 

0.3719 

4.2762 

0.4732 

0.3718 

4.2868 

90 

0.3813 

0.3508 

4.7337 

0.3807 

0.3513 

4.7261 

0.3807 

0.3513 

4.7364 

Figures  29, 30  and  31  show  the  variation  of  normalized  mode  I  stress  intensity  factor 
(Kj)  as  a  function  of  kink  length  (/)  and  kink  angle  (0)  for  different  stiffness  ratio.  From  these 
plots  the  following  can  be  observed  regarding  kink  tip  closing  angle  (which  is  defined  as  that 
kink  angle  at  which  the  mode  I  stress  intensity  factor  becomes  zero).  The  kink  tip  closes  for 
only  a  certain  range  of  kink  angle  and  length.  The  kink  length  and  the  angle  at  which  the 
kink  closes  decreases  with  an  increase  in  stiffness  parallel  to  the  direction  of  the  main  crack. 
For  stiffness  ratio  between  1  and  10  crack  closure  is  observed  only  for  kink  angles  above  60 
degrees.  If  the  length  and  angle  of  the  kink  is  large  enough  the  kink  will  not  close  at  all  for 
uniaxial  loading. 

Figure  31  shows  the  behavior  of  mode  I  stress  intensity  factor  when  the  stiffer 
direction  is  parallel  to  the  main  crack.  It  can  be  observed,  that  for  very  short  kink  lengths  (/ 
approaching  0),  the  maximum  Ki  does  not  occur  at  0  =  0  but  at  kink  angles  other  than  zero. 
This  implies  that  once  the  crack  turns  to  form  a  kink  it  requires  a  relatively  lower  load  to 
extend  the  kink  further.  But  as  the  kink  length  increases  Kj  starts  decreasing,  which  means 
a  higher  load  is  now  required  to  extend  the  kink.  Gao,  et  al.,  [56]  reported  that  for  infinitesi¬ 
mally  small  kink  lengths,  for  stiffness  ratio  more  than  4  (stiffer  axis  along  the  main  crack 
direction),  the  maximum  Kj  did  not  necessarily  occur  at  0  =  0.  The  present  result  confirms 
that  observation  and  also  shows  that  for  longer  kink  lengths  maximum  Kj  always  occurs  at 
0  =  0  no  matter  what  the  stiffness  ratio  is.  Gao,  et  al.,  [56]  also  observed,  for  an  infinitesimal 
small  kink,  that  when  the  stiffer  axis  was  along  the  main  crack  direction  and  the  stiffness 
ratio  was  four  the  predicted  kink  growth  direction  based  on  the  stress  intensity  factors  was 
physically  unrealistic.  Due  to  this  reason  Gao,  et  al.,  [56]  concluded  that  the  stress  intensity 
based  fracture  criteria  may  not  be  accurate  for  infinitesimally  small  kink  lengths.  For  longer 
kink  lengths,  from  Figures  29, 30  and  3 1 ,  it  can  be  observed  that  maximum  Ki  occurs  always 
at  0  =  0  irrespective  of  fibre  orientation  and  stiffness  ratio.  But  it  has  been  experimentally 
observed  for  polymeric  composites  (like  graphite-epoxy)  that  the  crack  growth  direction  is 
sensitive  to  the  fibre  orientation  direction  since  the  strength  of  the  composite  also  varies  ac- 


cording  to  the  orientation  of  the  fibre.  Hence  it  can  be  concluded  that  a  fracture  criteriabased 
on  Kj  alone  might  not  be  accurate  for  predicting  fracture  in  an  orthotropic  body. 

Figures  32, 33  and  34  shows  the  variation  of  normalized  mode  II  stress  intensity  factor 
Kn  as  a  function  of  kink  length  and  angle  for  different  stiffness  ratio.  Comparing  these  plots 
with  Figure  22  it  can  be  inferred  that  for  longer  kink  lengths  the  behavior  of  Kn  is  similar 
to  that  observed  for  the  isotropic  case  whereas  for  shorter  kink  lengths  Kn  is  affected  by  the 
variation  in  stiffness  ratio  and  hence  behaves  differently. 

Figures  35, 36  and  37  shows  the  variation  of  normalized  T-stress  as  a  function  of  kink 
length  and  angle  for  different  stiffness  ratio.  A  Tcrit/o  value  of  0.5  (which  is  an  assumption) 
is  also  shown  on  these  figures.  The  Tent  for  an  orthotropic  material  has  to  be  experimentally 
determined  as  in  the  case  of  an  isotropic  material  like  PMMA.  Since  experimental  Tent  val¬ 
ues,  for  any  orthotropic  material,  are  not  available  a  value  of  0.5  ksi  is  assumed  for  all  fibre 
orientation  with  respect  to  the  crack.  It  should  be  noted  that  Tcrit  varies  with  fibre  orientation, 
which  is  not  incorporated  into  the  current  assumption  of  T^it  values.  However,  experiments 
can  be  performed  on  specific  orthotropic  materials  to  determine  these  values  which  can  be 
incorporated  into  the  model.  From  these  plots  it  can  be  observed  that  for  short  kink  lengths 
the  T-stress  is  greater  than  Tcrit  for  all  kink  angles  and  stiffness  ratios.  But  for  longer  kink 
lengths  the  T-stress  can  be  either  above  or  below  Tcrit  depending  on  the  kink  angle,  and  the 
material  stiffness  ratio.  Figure  35  shows  the  variation  of  T-stress  for  the  case  when  the  stiff- 
er  direction  is  perpendicular  to  the  direction  of  the  main  crack.  Using  the  modified  T-stress 
based  fracture  criterion,  if  a  kink  is  formed  at  a  direction  to  the  main  crack  and  if  the  kink 
length  is  small  then  it  is  likely  that  the  kink  will  not  run  in  its  initial  direction  but  will  turn 
away  from  that  direction  since  the  T-stress  is  greater  than  Tcrit-  If  the  kink  length  is  long 
enough  (more  than  30%  of  the  main  crack  length)  then  the  kink  can  extend  along  its  initial 
direction  up  to  an  kink  angle  of  80  degrees. 


Figure  29.  Variation  of  Normalized  Ki  with  Normalized  Kink  Length  /  (k  =  0,  E1/E2  =10, 
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Figure  36  shows  the  variation  of  normalized  T-stress  for  the  case  when  the  stiffness 
ratio  is  1.1.  The  behavior  of  the  T-stress  is  similar  to  that  observed  in  the  isotropic  case. 
Figure  37  shows  the  variation  of  normalized  T-stress  when  the  stiffer  axis  is  parallel  to  the 
main  crack  axis.  From  this  plot  it  can  be  inferred  that  for  short  kink  lengths,  since  the  T- 
stress  is  greater  than  Tcru,  the  kink  does  not  run  along  the  initial  kink  growth  direction.  For 
longer  kink  lengths  (more  than  20%  of  the  main  crack  length)  and  for  kink  angles  below  30 
degrees  the  T-stress  is  less  than  Xcrit*  This  means  that  for  sufficiently  longer  kink  lengths 
the  kink  can  run  in  its  initial  kinking  direction  up  to  an  angle  of  30  degrees.  For  angles  above 
30  degrees,  since  the  T-stress  is  greater  than  Tcrit,  the  kink  turns  from  its  initial  direction. 

Figures  38,  39  and  40  shows  the  variation  of  normalized  Ki,  Kn  and  T-stress  respec¬ 
tively  as  a  function  of  material  angle  for  a  kink  length  which  is  0.2%  of  the  main  crack  length 
and  a  stiffness  ratio  of  10.  From  Figures  38  and  39  it  can  be  inferred  that  when  the  kink  angle 
is  zero,  i.e.  when  the  kink  coincides  with  the  axis  of  the  main  crack  Kj  and  Kn  are  insensitive 
to  the  variation  in  material  angle.  For  all  other  kink  angles  there  is  some  variation  in  Kj  and 
Kn-  Figure  38  shows,  as  stated  earlier,  that  it  is  possible  for  Kj  to  be  maximum  at  kink  angles 
other  than  zero.  Figure  39  shows  that  Kn  can  also  be  zero  at  kink  angles  other  than  zero. 
Figure  40  shows  some  interesting  variations  in  T-stress  as  compared  to  Kj  and  Kn-  The  T- 
stress  is  sensitive  to  variations  in  material  angle  for  a  kink  angle  of  zero  degrees,  i.e.  when 
the  axis  of  the  kink  coincides  with  the  axis  of  the  main  crack.  For  all  other  kink  angles  the 
T-stress  shows  large  fluctuations  as  the  material  angle  is  varied.  Up  to  a  kink  angle  of  60 
degrees  the  maximum  T-stress  below  Tcrit  occurs  when  the  kink  almost  parallels  the  stiffer 
axis  direction.  For  example  the  maximum  T-stress  below  Tcrit  for  a  kink  angle  of  20  degrees 
occurs  when  the  material  angle  is  close  to  25  degrees.  For  kink  angles  more  than  60  degrees 
the  T-stress  remains  above  Tcrit  for  all  material  angles. 
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Albritton  and  Goree  [5]  performed  experiments  on  standard  compact  tension,  single 
edge  notched  tension  and  center  crack  tension  silicon  carbide  reinforced  aluminum  coupons. 
These  test  coupon  had  fibres  running  parallel  to  or  perpendicular  to  the  load  direction.  While 
crack  turning  was  observed  in  CT  test  coupons  where  the  fibres  ran  parallel  to  the  load  it  was 
not  observed  for  any  of  the  coupons  for  which  the  fibres  ran  perpendicular  to  the  load.  From 
Figures  35  and  37  it  can  be  seen  that  for  a  kink  angle  of  zero  degree  the  T-stress  has  a  maxi¬ 
mum  value  below  Tcrit  when  the  stiffer  axis  is  parallel  to  the  crack  (i.e.  perpendicular  to  the 
load).  Hence  the  crack  in  the  CT  test  coupon,  where  the  stiffer  fibres  are  perpendicular  to 
the  load,  runs  stably  when  compared  to  the  case  when  the  fibres  are  parallel  to  the  load.  Al¬ 
britton  and  Goree  [5]  also  noted  that  when  the  fibre  volume  fraction  increased  for  the  case 
where  the  fibres  ran  parallel  to  the  applied  load,  the  angle  at  which  the  crack  turned  in  the 
finite  coupon  increased.  This  dependence  of  crack  turning  angle  on  the  fibre  volume  fraction 
can  be  attributed  to  the  T-stress  as  explained  in  the  following  lines.  Increasing  fibre  volume 
fraction  primarily  results  in  increase  in  stiffness  along  the  fibre  direction.  From  Figures  35 
and  36  it  can  be  seen  that  as  the  stiffness  increases  in  a  direction  parallel  to  the  load  the  maxi¬ 
mum  kink  angle  at  which  the  T-stress  is  below  Tcrit  increases,  i.e.  whereas  for  E1/E2  =1.1 
the  T-stress  is  below  Tcrit  up  to  a  kink  angle  of  60  degrees,  for  E1/E2  =  10  (stiffer  axis  parallel 
to  the  load  direction)  the  T-stress  is  greater  than  Tcrit  up  to  a  kink  angle  of  80  degrees.  There¬ 
fore,  for  a  sufficiently  long  kink,  the  angle  at  which  the  kink  can  turn  and  run  increases  with 
the  increase  in  the  volume  fraction. 

It  should  also  be  noted  that  for  an  orthotropic  material  there  exists  orthotropy  in 
strength  also.  Therefore  the  crack  can  turn  run  in  a  direction  where  the  material  is  weak 
which  is  usually  parallel  to  the  stiffer  axis.  By  calculating  the  Tcrit  value  experimentally  the 
influence  of  strength  and  stiffness  of  a  material  are  also  incorporated  in  to  the  fracture  model 
through  the  Tcrit-  Hence  by  comparing  the  theoretically  determined  T-stress  values  with  Tcrit 
values  the  influence  of  strength  and  stiffness  are  also  taken  into  account. 
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The  current  model  is  a  macroscopic  model  where  phenomenon  such  as  fibre-matrix 
debonding,  fibre  or  matrix  cracking  cannot  be  accounted.  But  these  can  be  analyzed  by  using 
a  microscopic  model.  Therefore  this  model  should  be  used  with  caution  when  applied  to 
such  materials  where  the  above  mechanisms  dominate. 


CHAPTER  Vn 


CONCLUSIONS 

The  branched  crack  problem  in  an  isotropic  and  anisotropic  material  is  solved  using 
the  method  of  dislocations.  Stress  singularities  at  the  comer  where  the  branched  crack  and 
the  main  crack  meet  is  investigated  and  the  singularity  is  shown  to  be  given  by  Williams’ 
[19]  characteristic  equation  for  the  isotropic  case.  An  accurate  method  is  developed  to  evalu¬ 
ate  the  T-stress  at  the  tip  of  the  branched  crack.  It  is  shown  that  the  T-stress  at  the  branched 
crack  tip  is  insensitive  to  the  order  of  the  singularity  at  the  comer  where  the  kink  and  main 
crack  meet,  for  both  isotropic  and  anisotropic  materials.  It  is  also  verified  that  the  mode  I 
(Ki)  and  II  (Kn)  stress  intensity  factors  at  the  kink  tip  are  insensitive  to  the  stress  singularity 
at  the  reentrant  wedge  comer  of  the  branched  crack  as  weU. 

The  values  of  Kj,  Kn  and  the  T-stress  evaluated  in  front  of  the  kink  tip  are  studied  as 
a  function  of  kink  to  main  crack  length  ratio,  kink  angle  and  the  loading  parameters  for  the 
isotropic  case;  for  the  anisotropic  case  the  influence  of  relative  stiffness  property  is  also  in¬ 
cluded.  Based  on  the  experiments  performed  by  Richardson  [94],  Cotterell  and  Rice’s  [7] 
T-stress  criteria  is  modified  to  include  an  experimentally  determined  value  of  Tent  above 
which  a  crack  deviates  from  its  plane  of  growth.  It  is  proposed  that  a  crack  will  continue 
in  its  plane  of  growth  if  the  T-stress  is  less  than  Tent,  or  deviate  from  its  plane  if  the  T-stress 
is  more  than  Tcrit-  For  PMMA  (which  is  isotropic),  from  Richardson  [94],  a  Tcrit  value  of 
0.5  ksi  was  experimentally  determined. 

The  modified  T-stress  based  fracture  criteria  is  applied  to  the  branched  crack,  in  both 
isotropic  and  orthotropic  material,  and  the  behavior  of  the  branched  crack  is  discussed.  The 
following  conclusions  are  arrived  at  from  the  above  study.  Considering  an  isotropic  material 
first,  it  is  observed  that  if  a  crack  subjected  to  either  uniaxial  or  biaxial  loading  deviates  from 
its  path  to  form  a  kink,  the  length  and  angle  of  the  kink  with  respect  to  the  main  crack  deter- 
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mines  how  the  kink  behaves.  For  a  short  initial  kink  length  (/  tending  to  0)  and  for  all  kink 
angles  (0<6<90),  under  uniaxial  loading,  the  kink  aligns  itself  with  the  main  crack,  but  under 
biaxial  loading  it  deviates  away  from  the  main  crack.  Therefore,  for  a  biaxial  loading  state, 
the  kink  direction  is  sensitive  to  the  applied  transverse  load  Oxx°°‘  For  a  longer  kink  length, 
the  kink  can  continue  growing  in  the  initial  direction.  For  example  consider  an  isotropic  ma¬ 
terial  like  PMMA  (Tcrit  =  0.5  ksi) ;  for  a  kink  angle  of  less  than  60  degrees  and  a  kink  length 
which  is  more  than  40%  of  the  main  crack  length  the  kink  can  extend  in  the  initial  direction 
it  started. 

Unlike  the  isotropic  case,  Tent  for  an  orthotropic  material  is  sensitive  to  the  material 
stiffness  and  strength  of  the  material,  i.e.  the  Tent  varies  as  a  function  of  orientation  of  the 
fibre  with  respect  to  the  applied  load.  Hence  a  range  of  Tent  values  has  to  be  experimentally 
determined  in  order  to  apply  this  modified  fracture  model  to  the  orthotropic  material.  In  this 
study,  for  an  orthotropic  material  (due  to  the  lack  of  experimental  Tcrit  values)  a  Tcrit  value 
of  0.5  ksi  is  assumed  for  all  fibre  orientation.  Hence  the  following  conclusions  are  arrived 
at  from  the  above  assumptions.  The  behavior  of  short  kinks  ( /  tending  to  0)  is  similar  to  the 
isotropic  case  for  both  uniaxial  and  biaxial  loading  conditions.  Hence,  it  can  be  inferred  that, 
the  relative  stiffness  does  not  influence  the  kink  behavior  when  the  kink  is  small.  On  the  con¬ 
trary  for  a  longer  kink  length,  under  uniaxial  or  biaxial  loading,  the  behavior  of  the  kink  is 
influenced  by  the  relative  stiffness.  For  an  uniaxial  load,  higher  the  stiffness,  more  is  the 
tendency  for  the  kink  to  run  parallel  to  the  stiffer  direction. 

Hence,  given  a  material,  the  relative  dimensions  of  the  crack  and  the  kink,  and  the  Tcrit 
value  (which  can  be  determined  as  in  Richardson  [94]),  the  present  study  can  be  used  to  pre¬ 
dict  how  the  kink  will  behave  when  subjected  to  different  loading  conditions  based  on  the 
modified  T-stress  fracture  model. 


APPENDICIES 


86 


APPENDIX  A 
ISOTROPIC  MATERIAL 
Analytical  Formulation  for  a  Dislocation  Function 
Consider  an  infinite  plane,  as  shown  in  Figure  A.  1 ,  with  a  dislocation  at  point  Zq  with 
respect  to  the  principal  axes  (X-Y).  Fix  another  co-ordinate  system  (X’-Y’)  with  Zq  as  ori¬ 
gin.  The  most  general  form  of  the  dislocation  stress  function,  which  is  multi  valued,  can  be 
written,  following  Muskhelishvili  [15],  as  follows. 


Y  Y’ 


Figure  A.l.  Geometry  of  Dislocation  Function 


4>d(z)  =  Yiln(z')- 

(A.1) 

(A.2) 
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where,  71  and  72  are  complex  constants,  z  and  z’  are  complex,  z  is  written  either  as  te‘®^  or 
xe‘®^.  (‘t’  and  ‘x’  are  used  interchangeably  in  the  development  of  equations  below).  The 
branch  cut  for  the  multi  valued  function  can  be  taken  to  be  along  the  X’  axis.  The  displace¬ 
ment  can  be  be  expressed  in  terms  of  the  stress  functions  as, 

u  +  iv  =  ^[x(j)D(z')  -  z'<I>d(z')  -  x1)d(z')]’  (A.3) 

where,  E= Youngs  modulus.  For  one  circuit  (Lj)  around  the  dislocation,  from  the  side 
to  the  '+’  side  in  anti-clockwise  direction,  the  net  change  in  the  displacement  is  given  as, 

u  +  iv  =  •^[>t(|)D(z')  -  zO^Cz')  -  ^d(z')]  ’  (A-4) 

which,  after  substituting  for  the  stress  functions,  can  be  written  as, 

[u'^  -  u~]  -f  i[v'^  -  v“]  =  ^[X7i  -1-  72]-  (A.5) 

Let, 

[b]  =  [u+  -  U-]  4- i[v+  -  V-],  (A.6) 

where,  [b]  is  the  Burgers  vector.  Therefore, 

f  [X7i  +72]  =  [b]-  (A.7) 

The  sum  of  external  forces  acting  on  a  contour  (Li)  around  the  dislocation  should  be 
zero.  This  is  written  in  terms  of  the  stress  functions  as, 

^  Forces  =  -  i[4)D(z')  -  z'^^(z')  +  =  0’ 

which  implies  that. 


Yi  =  Y2-  (A.9) 

Now  the  constants  71  and  72  are  uniquely  determined  and  the  stress  functions  can  be  ex¬ 
pressed  as. 


“  :ti(x  +  1)^”^^'^’ 


(A.  10) 
(A.  11) 
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Changing  the  coordinate  system  from  X’-Y’  to  X-Y  coordinate  system,  the  stress 
functions  for  the  dislocation  can  be  written  as, 


^  ^  [b]E  ,  , 


^d(z)  =  -  r- 


[b]E 


■ln(z  -  Zo)  -  Zo 


[b]E 


:iTi(x  +  °jti(x  +  l)(z  —  Zo)’ 

by  using  the  expressions  given  in  Muskhelishvili  [15].  Let, 

0(t  )  =  WE 

where,  to  is  some  arbitrary  variable.  Therefore, 

(})d(z)  =  0(to)ln(z  -  Zo)’ 

^  ©(to) 


o^dCz)  =  0(to)ln(z  -  Zo)  -  Zo 


’(z-zo)’ 
(|)'d(z)  =  4)d(z)  =  0(to)/(z  -  Zo)> 


(A.12) 

(A.13) 

(A.14) 

(A.15) 

(A.16) 

(A.17) 

(A.18) 


tl^'D(z)  =  'Pd(z)  =  ©(to)/(z  -  Zo)  +  0(to)zo/(z  -  Zo)A 
An  infinitesimal  burgers  vector  can  be  mathematically  represented  as  [bjdto,  where 


[b]  =  -  u  ]  +  i[v'^  -  V  ]• 


(A.19) 


This  burgers  vector  can  be  integrated  over  a  range  to  represent  a  crack. 


Formulation  of  Stress  Functions  for  the  Main  Crack  Which  Opens  Up  Due  to  the  Pres¬ 
ence  of  a  Dislocation 

The  stress  fields  on  the  main  crack  due  to  a  dislocation  away  from  the  crack  are  next 
evaluated  by  using  the  dislocation  stress  stress  functions.  The  normal  and  shear  stresses  act¬ 
ing  on  the  surface  of  the  main  crack  due  to  a  dislocation  are  given  as. 


o„(x,  0)  -  io,,(x.  0)  =  ^  (A,20) 

where,  -2c  <  x  <  0.  The  stresses  acting  on  the  main  crack  surface  is  now  known.  Hence 
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the  boundary  value  problem  of  a  pressurized  crack,  shown  in  Figure  A.2,  can  be  solved  using 
techniques  given  in  Muskhelishvili  [15]. 


Y 


Figure  A.2.  Geometry  of  Pressurized  Main  Crack 

Following  Muskhelishvili  [15]  the  normal  and  shear  stresses  due  to  the  main  crack 
alone  can  be  written  as, 

Oyy(z)  -  iaxy(z)  =  Om(z)  +  d)M(z)  +  z«)'m(z)  +  (A.21) 

Let, 

Qm(z)  =  (A-22) 

Then, 

c^yy(^)  ~  icrxy(z)  =  Oj^(z)  +  Qj^(z)  +  (z  —  (z)-  (A.23) 

Let  +  refer  to  the  upper  side  of  the  main  crack  and  -  refer  to  the  lower  side.  Using 
equation  A.23  the  stresses  on  the  upper  and  lower  side  of  the  main  crack  can  be  written  as. 
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a+y  -  iOx+y  =  0+(x)  +  Quix),  (A.24) 

(Jyy  ~  i(?xy  ~  ^m(x)  +  Qj^(x)-  (A. 25) 

Adding  equations  A.24  and  A.25, 

^yy  ^yy]  ~  ^[t^xy  tJxy]  =  [<I*j^(x)  +  Qjyj(x)]  +  [Ojyf(x)  +  Qf^(x)]  •  (A. 26) 
Similarly  subtracting  equation  A.24  from  A.25, 

2  ^  ^yy  ~  ”  2  ~  ~  “■  ^  ■” 

Using  theorems  given  in  Muskhelishvili  [15]  A.27  can  be  written  as, 

[OmCz)  -  Qm(z)]  =  0-  (A.28) 

Using  A.26,  A.20  and  A.28, 


^m(z)  ~ 


1  f  X  +  (t)p(t,Zo)dt  c 
2jriX(z)  J  t  -  z 

L 


(A.29) 


where,  Om(z)  is  a  sectionally  holomorphic  function  which  is  holomorphic  everywhere  ex¬ 
cept  on  L  (0  to  -2c).  C  is  a  constant  which  has  to  be  determined  by  using  the  condition  of 
single  valuedness  of  displacement.  Also  note  that  x  has  been  replaced  by  t.  The  above  func¬ 
tion  Om(z)  can  be  evaluated  by  using  contour  integrals.  As  an  example,  one  term  in  p(t,Zo) 
will  be  evaluated  below. 

Consider  the  first  term  in  the  above  integral  (A.29)  which  is  repeated  below  with  all 
the  terms. 


1 


2KiX(z) 


fxi 

t  - 
J 


(t) 

z 


0(to) 


0(to)  Q(lo)(Zo  Zq) 


t  -  Zn 


(t  -  Zo)2 


dt  -h 


X(z) 


(A.30) 


F(z,  Zo) 


1 


27i:iX(z) 


fx+(t) 

©(to)  1 

1  t  —  z 

t  -  Zo 

dt) 


(A.31) 


L 

where  X(z)  is  a  sectionally  holomorphic  function  and  its  branch  cut  is  defined  as  shown  in 
Figure  A.3. 
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Figure  A. 3.  Definition  of  the  Function  X(z) 

X(z)  =  yz(z  +  2c)-  (A.32) 

X(z)  =  (A.33) 

The  branch  cut  for  is  taken  such  that,  0  <  0a  <  27t  and  similarly  the  branch  cut  for 
z+2c  is  taken  such  that  0  <  Bj  <  27t.  Therefore, 

X+(t)  =  yt(2c  +  t)>  (A.34) 

X-(t)  =  -  yt(2c  +  t)-  (A.35) 

Consider  a  contour  around  the  line  of  integration  as  shown  in  Figure  A.4  in  which  the 
contour  F  divides  the  entire  infinite  region  into  two  parts.  S’*"  and  S“. 


Figure  A.4.  Contour  F 

^  is  complex  and  lies  on  the  contour  F.  The  line  integral,  equation  (A.31)  can  be  writ¬ 
ten  in  terms  of  a  contour  integral  as  shown  below. 
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i 

■  0(to)  ■ 

d^  =  ( 

c 

X(|) 

■  0(to)  ■ 

d^+- 

L- 

rx+(t) 

r  0(to)  1 

J  ?  -z 

r 

3  ~ 

^  1  -  z 

1 

3  ~ 

to  -  z 

1- 

t  -  Zo 

+ 


'  0(to)  ‘ 

d£+ 

0(to)  1 

J  ^  -z 

-  Zo_ 

+  to  -  Z 

t  -  Zo 

dt. 


(A.36) 


Cj  L- 

It  can  be  shown  that  the  contour  integrals  around  Ci  and  C2  vanish  as  the  contours  are 
shrunk  to  a  point.  Hence  this  leaves  the  line  integrals.  Since  X+(t)  and  X“(t)  are  related,  the 
contour  integral  around  T  can  be  written  as, 


X(|) 

'  0(to)  ' 

d|  =  - 

J 

fx+cto) 

■  0(to)  1 

dt  - 

- 

fx-(to) 

■0(to)l 

)  . 

5  —  Z 

I  -  Zo_ 

to  -  z 

t  -  Zo 

to  -  z 

t  -  Zo 

r  L+  L- 


(A.37) 


r 


Q(to)l  ,,.  ^  ^  f  X  +  (to) 


e(to) 

t  -  Zo 


dt- 


(A.38) 


The  contour  integrals  can  be  evaluated  using  the  theorems  given  in  Muskhelishvili 
[15],  chapter  12.  Let, 


'  0(to)  1 

J  ^ 

0 

SI 

1 

JXP 

r 

where  f(|)  is  defined  as. 

(A.39) 


(A.40) 

f(^)  is  a  sectionally  holomorphic  function  in  S+ ,  but  holomorphic  everywhere  in  S“ 
except  at  Zq  where  it  has  a  simple  pole.  Zq  is  assumed  to  lie  outside  the  contour.  Using  the 
theorems  given  for  contour  integrals  in  Muskhelishvili  [15],  chapter  12,  the  following  can 
be  written. 


m)  =  x(i) 


e(to) 
i  -  Zo 
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1  I  f(^) 


2jti  T  t  -  z 


=  f(z)  -  Gi(z)  -  Goo(z), 


(A.41) 


where,  Gi(z)  is  the  principal  part  of  the  function  at  Zq  and  Goo(z)  is  the  principal  part  at  infin¬ 
ity.  Note  the  change  in  sign  in  equation  A.41  when  compared  to  the  theorem  given  in  Musk- 
helishvili  [15],  chapter  12,  since  the  direction  of  integration  considered  for  this  study  is  re¬ 


versed.  Gi(z)  is  evaluated  as  follows. 

Ai 

Gi(z)  = 

where  A]  is  the  residue  of  the  function  f©  at  Zq. 

,X(^)0(to) 

A  =  lim(i  -  Zo)-^ - 

‘  5-Zo  (?  ~  Zo) 

A,  =  X(zo)0(to)- 


(A.42) 


(A.43) 


(A.44) 


To  evaluate  Goo(z)  the  function  f(|)  is  expanded  in  terms  of  1/|  and  evaluated  as  ^  ap¬ 


proaches  oo. 


Goo(z)  =  lim  1  +  ^  1  -  ®(to) 


Goo(z)  =  0(to)- 


(A.45) 


(A.46) 


Therefore, 


-  0(.„> 


2n;iT  ^  -  z 


(A.47) 


Using  similar  procedures,  all  the  line  integrals  given  in  equation  A.30  can  be  evaluated 


to  give. 


<E)m(z)  =  -  0(to)F*(z,  Zo)  +  0(to)F  (z,Zo)  +  0(to)(zo  -  Zo)G  (z,Zo))  + 


(A.48) 


where. 


P*.  s  _  i  _1 _ X(Zo) _ 

t  {z,Zo)  2Z-Zo  X(Z)(z  -  Zo)  X(Z) 


(A.49) 
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G*(z,  Zo)  =  ^o)’ 


(A.50) 


and  C  is  still  unknown  and  can  be  evaluated  from  the  condition  that  the  displacements  are 
single  valued.  The  condition  for  single  valuedness  of  displacement  is  given  as, 


x|[Od(z)  +  ^M(z)]dz  -  |[Qd(z)  +  ^M(z)]dz  =  0- 


(A.51) 

I  ''  iVlN 

r  r 

Contracting  these  contours  along  the  line  of  discontinuity  0  to  -2c  and  0  to  / 

(Figure  A.5),  the  contour  integrals  can  be  written  (using  the  relation  given  in  A.28)  as, 


(x+  1) 


[<Dj:i(t*)  -  <DM(t*)]dt*  +  X 


-  OD(t*)]dt* 


Lj  +  L2 


Li  + 


+  I  [Q"(t*)  -  Q;(t*)]dt*  =  0, 


(A.52) 


Lj+Lj 


Figure  A.5.  Contour  F  (L1+L2) 


where  the  +  and  -  are  the  limiting  values  of  the  sectionally  holomorphic  function  on  the  line 
of  discontinuity.  Od  and  Qd  are  defined  as, 

^,(z)  =  (A.53) 


Qd(z)  =  Od(z)  +  z^>'d(z)  +  Wd(z). 


therefore. 


^d(^)  ~ 


©(to) 
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(A.54) 

(A.55) 


(z  -  Zo) 

Using  +  and  -  limiting  values  of  the  function  in  the  first  term  in  equation  A.52,  the 
following  is  obtained. 

0 

(x+1)  [^^(t*)  -  Oj;,(t*)]dt*  =  (x  +  1) 

L|  t  =  —  2c 


X(Zo)©(to)  X(Zo)0(to) 

X  +  (t)(t  -  Zo)  X  +  (t)(t  -  Zo) 


2©(to)  ©(to)(Zo  +  c)(Zo  -  Zo)  X(Zo)©(to)(Zo  -  Zq)  2C 
X  +  (t)  X(Zo)(z  -  Zo)X  +  (t)  X  +  (t)(t-Zo)2  x  +  (t) 


dt  (A.56) 


where,  t*  =  t.  The  above  integrals  can  be  evaluated  using  contour  integrals  (see  Appendix 
D,  equations  D.l  -  D.13)  and  are  given  as, 

0 

(A.57) 


dt 


Jii 


X  +  (t)(t  -  Zo)  X(Zo) 


t=  -2c 
0 

1 

t=  -2c 


dt  _  _ 
X  +  (t) 


=  —  jti 


(A.58) 


Substituting  the  above  integrals  in  equation  (A.56)  the  following  is  obtained. 
(x  +  1) 


-  ^)M(t*)]dt*  =  (X  +  1)(-  4xiC  -  2m©(to))-  (A.59) 


The  only  terms  from  Om(z)  which  contributes  on  L2  are  the  following  integrals. 


/ 

-4 


©(to) 


-dt^  +  ^ 


1  X(Zo)©(to) 


(t*-Zo)  2  X  +  (t*)(t*-Zo) 


dt 


where  t*  e  L2.  Their  +  and  -  limiting  values,  using  Plemelj  formula,  are  given  below. 
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X(zo)0(to) 
X(t*)(t*  -  Zo) 


dt*  =  X(zo)0(to) 


/ 

I 


^  --  -dt* . 


X(t*)(t*  -  Zo) 


0 

0 

+  , 

1 

r 

/ 

r 

=  Jti0(to)  +  X(zo)0(to) 

0 

0 

1  -dt*. 


X(t*)(t*  -  Zo) 


X(Zo)0(to) 


and 


1 


X(t*)(t*  -  Zo) 


-dt 


=  -  3ti0(to)  -I-  X(zo)0(to) 


1 


X(t*)(t*  -  Zo) 


■dt* 


-dt*  =  0(to) 


(t  *  -  Zo) 


dt* 


(t*  -  Zo) 


0(to) 


0(to) 


I 

L, 

I 


dt* 


(t*  -  Zo) 


dt* 


=  ra0(to)  -f  bf% 


(t  -  Zo) 


=  -  jri0(to)  -K  bf  , 


where  bf*  is  a  bounded  function.  Therefore, 


2 


(X  -h  i)i  [a>j:i(t*)dz  -  OM(t*)]dt*  =  0 


(A.60) 


Consider  the  next  integral 


xj  [O^Ct*)  -  Opi 


,(t  )]dt  -f 


[Q^O*)  -  QD(t*)]dt* 


L|  +  L2 


+L2 


These  integrals  do  not  have  different  limiting  values  on  Li  but  only  on  L2.  Therefore  using 
Plemelj  formula  again  the  limiting  values  of  the  above  integrals  are  given  below. 
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<I>D(z)dz  =  0(to) 


(Z  -  Zo) 


,  X  =  Jti0(to)  +  bf% 


®(to)  .  =  -  3ti0(to)  +  bf, 


where  bf*  is  a  bounded  function.  Therefore, 


K  )dz  -  ^>j3(t  )]dt  =  x2jti0(to)- 


Similarly, 


=  2jti0(to)- 


Therefore, 


X  I  [OpCtOdz  -  <DD(tl]dt^  +  I  [QD(t*)  -  ^D(t*)]dt* 


=  (x  +  l)23Ti0(to)- 

Adding  equations  A.59,  A.60  and  A.61  and  equating  to  zero, 
(x  +  1)(-  2jtiC  -  4jri0(to)  +  2jii0(to))  =  0- 
C  =  -  0(to)- 


(A.61) 


(A.62) 

(A.63) 


The  constant  C  can  be  evaluated  using  a  different  method  which  is  the  same  as  that  used 
by  Lo  [24] .  Consider  equation  A.5 1  in  which  the  contour  is  assumed  to  enclose  the  branched 
or  kinked  crack  also.  Therefore, 
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H  (()[Od(z)  +  <I>M(z)]dz  - 


^[Qd(z)  + 


z  =  0- 


(A.64) 


where  the  contour  T  extends  to  oo.  Expanding  the  above  integrals  in  terms  of  1/z  and  1/z, 


+  ^'D(z)]dz  =  Jim  +  0(to)dto]  +  o^-^jjdz. 

r  r  ^ 

(|)[Qd(z)  +  Q;^(z)]dz  =  =  lim 
J  Z-*oo 


r 

If, 


=  [C  +  0(to)dto]  +  o(-^ 
^  Vz^ 


(A.65) 


(A.66) 


C  =  -  0(to),  (A.67) 

then  the  contour  integrals,  in  both  of  the  above  equations,  vanishes  as  z  tends  to  oo.  There¬ 
fore  the  functions  Om(z),<1>d(z),  'Pm(z)  and  ¥d(z)  are  single  valued  in  the  entire  plane,  in¬ 
cluding  infinity,  except  along  the  main  and  kinked  crack  respectively. 

Using  equation  (A.67)  in  (A.48), 

^m(z)  =  -  j0(to)F(z,Zo)  +  0(to)F(z,Zo)  +  0(t;)(zo  -  Zo)G(z,Zo)l-  (A.68) 

where. 


_J _ X(Zo) 

z  -  Zo  X(Z)(z  -  zo) 

G(z.  Zo)  =  ■^FCz,  Zo)» 
and 

^m(z)  =  +  z^^  -b 

therefore, 

Vm(z)  =  ^m(z)  -  Om(z)  -  zd)M'(z). 


F(z,  Zo)  =  i 


(A.69) 

(A.70) 


(A.71) 

(A.72) 


Formulation  of  the  Solution  for  a  Infinite  Plane  with  a  Crack. 


The  next  step  is  to  derive  the  equations  for  an  infinite  plate  with  a  crack  which  is  sub¬ 
jected  to  remote  infinite  biaxial  stresses.  The  steps  for  solving  this  problem  follows  directly 
from  Muskhelishvili  [15],  chapter.  19  (pg.  506  -  508). 


Pj(z)  =  Clo^  +  Clp 


r  - 
’-LO  -  2 


(A.73) 

(A.74) 

(A.75) 


Cli  is  a  constant  which  has  to  be  determined  from  single  valuedness  of  displacement 

condition.  Single  valuedness  of  displacement  condition  is  expressed  as, 

0 


Pi(t) 
X  +  (t) 


dt  =  0- 


-2c 

0 

1 

-2c 


+  ^Ll 

X  +  (t) 


dt  =  0- 


(A.76) 


(A.77) 


The  above  integral  can  be  evaluated  using  contour  integrals  (see  Appendix  D,  D.l 
D.25)  and  the  following,  Cli=Ci^c  is  obtained.  Therefore, 

CO" 


Cli  = 


VWyy 


<1>l(z)  = 


a“  (z  -f  c)  (a“ -a' 


Ql(z) 


Therefore, 


'Pl(z)  = 


J 

2X(z) 

+  — 

Ql(z)  -  OJlCz)  -  zOL(z) 

Oy“((Z  +  C) 

(Ox“x-a;y) 

2X(z) 

4 

ZC^Gyy 

(Ox“x-Oy”) 

(A.78) 

(A.79) 

(A.80) 

(A.81) 


2[X(z)]3/2 


(A.82) 
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Let, 


=  koy”. 

(A.83) 

2X(z)  4 

(A.84) 

,,,  ""yO'  -  1) 

^  2[X(z)]3/2  2 

(A.85) 

Formulation  of  the  Solution  for  the  Kinked  Crack. 

The  final  step  is  to  satisfy  the  boundary  condition  on  the  kink.  For  this,  the  Plemelj 
formulae  are  used.  The  dislocation  stress  functions  are  the  only  ones  which  give  different 
limiting  values  on  the  kink  surfaces.  The  stresses  due  to  the  dislocation  are  given  by, 

(z)  +  iag(z)  =  Od(z)  +  +  Wd(z)].  (A.86) 

The  *+’  and  limiting  values  of  the  stresses  on  the  kink  are  given  by. 


o^0(t)  +  ict^(t)  —  T  JTi©(t)  + 


©(to)e 


-ie 


(t  -  to) 


■dto  ±  Jti©(t)  + 


i 

J  (T^*" 


±  Jti©(t)  + 


i 

{  0(to)e-i9 

J  (to  -  t)  ‘ 


"dto  ^  7ri0(t)  + 


©(to)ei9 
(t  -  to) 


dto 


(A.87) 


tJ00(t)  +  ioS(t) 


re'' 


“  -  -._[®(to)e“*®  , 

-  7-5^*° 


The  above  integral  is  interpreted  as  a  Principal  Value  integral.  The  other  stress  func¬ 
tions  are  regular  on  the  kink.  Hence  the  stresses  on  the  kink  can  be  written  as, 

i  ^ 

r  t 

aee(z)  +  ia,0( 


/  I 

9(2)  =  +  j  M(t,  to;  a(to))dto  +  cr^(z)  +  io^^Cz)  =  0- 


(A.  88) 


Where, 
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M(t,to,a(to))  =  <Dm(z)  +  +  e2i9[(z  -  z)<Dm'(z)  -  (Dm(z)  +  ^(z)],  (A.89) 

Oe”e(z)  +  ia^(z)  =  «)l(z)  +  +  e2i0[z<DL'(z)  +  Wl(z)].  (A.90) 


APPENDIX  B 


EVALUATION  OF  THE  SINGULAR  INTEGRALS 
First  the  formula  used  for  evaluating  the  singular  integral  is  given.  The  derivation  can 
be  found  either  in  Muskhelishvili  [83],  chapter  4  or  Kaya  [87].  Then  the  application  of  the 
formula  to  some  typical  singular  integrals  is  shown. 


Formula 


0(to) 
to  -  Z 


where, 


0(to) 


P*(to) 

(to)“(/  -  to)P 


z  =  te'®%  Zo  =  toC^®  >  (B.3) 

and  P*(to)  is  a  bounded  function  in  the  interval  0  -  /  and  is  not  zero  at  the  end  points  (0  and 
/),  and  a,  P  <  1.  z  is  any  point  in  the  infinite  complex  plane  but  not  on  0  -  /  (Figure  1,  page 
2).  The  singular  behavior  of  the  above  integral  after  normalizing  the  integral  with  respect 
to  the  kink  length  /  as  z  approaches  0  is  given  as, 

f  f^k  =  -  ,i„X  +  (B.4) 

t^o  J  lo  -  zj  ®  2sm(:n;a)  t^ot 

where,  ©*(t)  tends  to  a  definite  limit  as  t  tends  to  zero  and  P(0)  =  P*(0)//“+P.  The  other  inte¬ 
grals  can  be  evaluated  in  a  similar  manner. 
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Some  Typical  Singular  Integrals 

For  all  the  integrals  given  below  the  length  parameter  is  normalized  with  respect  to 
the  kink  length  /  and  P(0)  = 


I,  =  limi  1  dto. 


‘1  -  ;“o2  (z  -  zo) 


t-o2  J  ta(l  -  to 


(tp)  dtp 

-  to)P  (z  -  Zp)’ 


t-o2  I  t»(l  -  t„)M  (te®.  -  t„e») 


limi  f  ^2^ 
t-o2  J  t»(l  -  to 


_  dtp 

)P  (tei(ea-e)  _  t^) 


which  can  be  written  as, 


1 

limi  _ 

t-o2  J  [ta 


v(tp)  dtp 

t«(l  -  to)P  (tp  -  Zi) 


where,  zi  =  te‘^®®  and  7(tp)  =  -P(to)e-*®.  Using  the  formula  given  above  (B .4)  for  singular 


integrals. 


-  3Tv(0)ei’^«  1  , 

— T — — — T  lim— 77  +  bii» 

—  2sin(3Ta)  t— qz?  ‘ 


where  bfi  is  bounded  function.  Substituting  for  7(0)  and  simplifying, 

JtP(0)e“4a(ea-Jt)]e-i[e(i-a)]  .  I 

— - - — — r -  lim-77  +  bfi- 

2  sin(jra)  t-*o  t  ‘ 


Consider  integral  I3. 


T  -  11  X(Zp)0(tp) 

I3  lini_  I  I  „  \  dto 


;ro2  X(z)(z  -  Zo) 


0 
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1  f  P(to)ytoe‘6(toe‘9  +  2c) 

=  1™:^  - ============== — ^ - — 

t^o  2  J  (tQ)“(l  —  to)P  yte‘®a(te‘0a  +  2c)  (te'®^  —  toC’®) 


which  can  be  written  as, 


=  lim-sS-  f  (3(to)y(toeie  +  2c) 

t^o2y^  J  (to)““^(l  —  to)l^(te‘®‘>  —  toC'®) 


=  lim-^  (  P(to)y(toe‘®  +  2c) 

t™2y2ct  J  (to)““2(l  -  to)P(tei®»  -  toei®) 


which  can  be  written  as. 


=  lim — i —  - - 

t-02y^J  _  yP  (to  -  z,) 


where,  zi  =  te'(®^‘®^  and 


v(to)  =  -  P(to)y(toe‘®  +  2c)e-i(^). 

Again,  using  the  formulae  given  above  (B.4)  for  singular  integrals  and  by  simplifying  the 

resulting  expressions,  the  following  is  obtained. 

_  JtP(0)ie-ti«(®a-^)]e-il®0  -«)]  ,•  J_  (B  6) 

2cos(jta) 

where  bf4  is  a  bounded  function.  Integrals  with  l/(z  -  Zq)^  terms  can  be  treated  as  follows. 


limljCz)  =  limi  f  Q(to)(Zo_  Zo)^^^^ 
,_»0  t-<«o2  I  (z  -  ZoY 


lirnIsCz)  =  - 


d  1  I  0(to)(zo  -  Zo) 


dz2  I  (z  —  Zq) 


The  singular  part  of  the  integral  can  now  be  evaluated  using  the  formula  for  singular 


integrals  and  is  given  below. 
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_  Ji:P(0)i(l  -  a)e-iW0a-Jt)]e-i[e(a-2)]  ^ 

- ^ +  bL» 

sin(jta)  t^o  t  ^ 

where  bfs  is  a  bounded  function.  All  the  other  integrals  can  be  evaluated  in  a  similar  manner 


as  given  above. 
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APPENDIX  C 

ANISOTROPIC  MATERIAL 

Formulation  of  Governing  Equations  for  an  Anisotropic  Elastic  Body 
The  Hooke’s  law  for  an  anisotropic  material  possessing  one  plane  of  elastic  symmetry  which 
is  perpendicular  to  the  z  axis  is  given  as, 

^xx  “  Cil^xx  d"  ^I2^yy  Cl6^xy  "h  Cis^^zz’  (0.1) 

(C.2) 

(0.3) 


Eyy  —  C2]Oxx  "I"  ^22^yy  ^26”^xy  "I"  C23OZZ’ 


Ezz  ~  CiS^xx  "i"  ^27pyy  ^36"^xy  "I"  C33OZZ’ 


2Exy  ~  CigCTxx  "h  ^2(Pyy  ^66'^xy’ 


(0.4) 


where,  Cij’s  are  the  material  constants.  If  plane  stress  conditions  are  assumed,  where  Ozz  = 
0,  the  following  equations  are  obtained. 

Exx  ~  CjjOx  +  ^I2^y  ^16'^xy’  (0.5) 

(0.6) 

(0.7) 


Eyy  ~  ^12^x  “I"  ^21^y  ^26'*'xy’ 


Ezz  ~  CJ3OXX  “f  E23Cfyy  +  C3gTxy’ 


2Exy  —  CjgOx  +  C25CJy  +  CggTxy- 


(0.8) 


If  plane  strain  conditions  are  assumed  then  £zz=0.  Solving  for  Ozz  using  equation  C.3  the 
following  relationships  are  obtained. 


Exx  &l]Oxx  "f  ^I2^yy  ^16^xy’ 


Eyy  —  aj2CJxx  "b  ^22®yy  ^26^xy’ 


2Exy  ~  “b  ^26^yy  ^66^xy' 


The  Cij’s  are  related  to  the  aij’s  by  the  following  relationships. 


ail  = 


C11C33  0^2 


"33 


ConC 


2T-33 


—  (^2 
23 


^22 


"33 


(0.9) 

(C.10) 

(0.11) 

(0.12) 

(C.13) 


ai2  - 


_  ^12^33  ^13^23 


"33 


(0.14) 
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„  _  ^26^33  ^23^36 

(C.15) 

^16^33  “  ^^13^^36 
^16  ’ 

(C.16) 

C66C33  -  cj^ 

““  c55  ' 

(C.17) 

1  The  equilibrium  equations  for  a  2-D  plane  problem  where  the  body  forces  are  ne- 

glected  are  given  by, 

dOxx  ,  f. 

-  "T  .  —  U? 

dx  dy 

(C.18) 

dOyy  dXxy 

dy  dx 

(C.19) 

The  compatibility  equation  is  given  by, 

ay2  dx^  dxdy  ' 

(C.20) 

Let, 

a2U(x,  y) 

"  ay2  ' 

(C.21) 

52U(x,  y) 
ax2  ■ 

(C.22) 

.  _  52U(x,y) 

dxay 

(C.23) 

These  equations  satisfy  the  equilibrium  equations.  Substituting  these  in  equations  C.5,  C.6 
and  C.8  and  then  substituting  them  back  in  the  compatibility  equation  C.20  the  following 

homogeneous  partial  differential  equation  is  obtained. 

a^u  0  a^u  54u  ^  34u  _  a4u 

"“ax'*  ^"“ax’ay  + ‘^"“^"«>ax2ay2  ^""^axay** 

=  0 

(C.24) 

The  solution  for  the  above  partial  differential  equation  can  be  written  as. 

U  =  F(x  +  sy) 

(C.25) 

Substituting  the  above  solution  in  equation  (C.24)  the  following  equation  is  obtained. 
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Ciis4  -  2Ci6s3  +  (2Ci2  +  C66)s2  -  2C26S  +  C22  =  0-  (C.26) 

Based  on  energy  principle  Lekhnitskii  has  shown  that  the  roots  of  the  above  equation 
are  always  complex.  Let  si,  S2,  S3  and  S4  be  the  four  complex  roots  (S3  =  si,  S4  =  $2)  of  the 
above  equation.  Then  the  complete  solution  of  partial  differential  equation  can  be  written 
as, 

U  =  Fi(x  +  Sjy)  +  F2(x  +  S2y)  +  F3(x  +  s^y)  +  F4(x  +  S4y)  (C.27) 


The  following  conditions  hold  for  the  complex  roots  Si,  i=1..4.  Let, 


Sj  =  ttj  +  iPi 

(C.28) 

§2  =  «2  +  iPz 

(C.29) 

where,  ai,  a2.  Pi,  P2  ^re  all  real. 

Pi  >  0,  P2  >  0  and  Pi  P2 

(C.30) 

Let, 

Zj  =  X  +  Sjy,  Z2  =  X  +  S2y 

(C.31) 

then. 

U  =  Fi(Zi)  +  F2(Z2)  +  Fi(Zi)  +  F2(Z2> 

(C.32) 

Let, 

dz', 

(C.33) 

1]  = 

(C.34) 

Using  C.32  to  C.34  in  C.21  to  C.23  the  stresses  can  be  written  as. 

2Rc[s^  (j)  (Zj^)  +  82^  (^2)3* 

(C.35) 

Oyy  =  2Re[(})'(Zl)  +  llj'(Z2)]- 

Txy  =  -  2Re[Si  (l)'(zi)  +  S2  ilj'(z2)]-  (C.37) 
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Substituting  the  above  relations  in  the  Hooke’s  law  (C.9-C.1 1)  and  after  integration  and  ne¬ 
glecting  the  rigid  body  displacements  and  rotations  the  following  displacement  equations  are 
obtained. 

u(x,y)  =  2Re[pj(j)(zi)  +  (C.38) 

v(x,y)  =  2Re[qi(j)(zi)  -f  q2^(z2)]-  (C.39) 

where, 

Pi  “  ^11^1  +  Ci2  -  CjgSi  • 

Pi  ~  ^11^2  *^12  “  ^16^2  ■ 

Ci2s2  -f  C22  -  C26S1 

= - s^ - 

_  ^^12^2  ^22  ~  ^26^2 


Analytical  Formulation  For  Dislocation  Function 
The  analytical  formulation  for  the  dislocation  function  is  given  below.  Consider  an 
infinite  plane,  as  shown  in  Figure  C.l,  with  a  dislocation  at  point  Zqi  (zoi  =  Xoi  +Siyoi,  i  = 
1  ..4)  with  respect  to  the  principal  axis  (X-Y).  The  most  general  form  of  the  dislocation  stress 
function,  which  is  multi  valued,  can  be  written  as  follows. 

(1)d(zi)  =  Aln(zj  -  Zoi>  (C.40) 

tji£)(z2)  =  Bln(z2  -  Zo2)*  (C.41) 

In  the  above  equations  (C.40,  C.41)  zi,  Z2,  Zqi,  Zo2,  A  and  B  are  complex.  The 
constants  A  and  B  can  be  determined  by  using  the  conditions  that  there  is  a  discontinuity  in 
displacements  for  one  revolution  around  a  contour  Li  which  goes  around  the  dislocation  and 
by  satisfying  the  condition  that  the  resultant  force  vector  around  a  contour  is  zero.  Let  the 
jump  in  discontinuity  in  displacement  be  represented  by  Bx  and  By  which  are  the  X  and  Y 
components  of  the  Burgers  vector.  The  ‘+’  and  superscript  represents  the  upper  and  lower 


side  of  the  dislocation  or  the  discontinuity  in  displacement.  The  resultant  force  vector  is  giv¬ 
en  by, 


X  +  iY  =  -  i[(l  -  is,)(|)D(Zi)  +  (1  -  iSj)(j)D(z,)  +  (1  +  is2)tiJD(z2)  +  (1  + 

(C.42) 

Using  equations  C.40,  C.41  in  C.42  the  following  equations  are  obtained. 

(1  +  iSi)A  -  (1  +  isi  )A  +  (1  -h  is2)B  +  (1  +  i%)B  =  0-  (C.43) 

(1  -  is,  )A  -  (1  -  iSi)A  +  (1  -  is2)B  -  (1  -  iS2)B  =  0-  (C.44) 


(p,  -  iq,  )A  -  (p,  -  iqi)A  +  (p2  -  iq2)B  -  (P2  -  iq2)B  =  — 


Bx  -  iBy 
2jT;i 


(C.45) 


__  _  Bx  "b  iBy  A^\ 

-  (p,  +  iq,  )A  +  (Pi  +  iqi)A  -  (p2  +  iq2)B  +  (P2  +  iq2)B  =  — — •  (C-46) 

The  above  four  equations  are  solved  using  Maple  [88]  and  the  following  expressions 
are  obtained  for  A  and  B. 

S](BxC22  —  S2S3S4CjjBy) 


A  = 


B  = 


(Sj  —  S2)(Sj  S3)(s,  S4)C22Cij 

S2(BxC22  ~  S,S3S4C,,By) 

(Sz  S,)(S2  ~  S3)(S2  ”  S4)C22Cjj 

If  Ci6  =  C26  =0  then, 

"22 


(C.47) 

(C.48) 


^1^2  —  ^1^2  ~ 


(C.49) 


Using  the  above  equation  in  the  expressions  for  A  and  B,  the  dislocation  density  function 
is  given  as, 

SjBx  -  By  ^  ^ 


(s,  -  SzKsj  -  SzXs,  -  84)0,1 
^0(^2)  -  _  s^)ci,  ^”^^2  2:02) 


(C.51) 


Therefore, 
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Od(Zi)  =  (j)'D(Z]) 


SjBx  By 


(Si  -  SjXSj  -  SjXSi  -  84X11  (Zj  -  Zqi) 


~  X  0(^2)  ■“ 


($2  -  Si)(S2  -  S3)(S2  -  84)011  (Z2  “  Z02) 


(C.52) 


(C.53) 


An  infinitesimal  Burgers  vector  can  be  mathematically  represented  as  [Bx  ]dto  and 


[Byjdto,  where 


(C.54) 


(C.55) 


This  Burgers  vector  can  be  integrated  over  a  range  to  represent  a  crack  in  that  range. 


Formulation  of  Stress  Functions  for  the  Main  Crack  Which  Opens  Up  Due  to  the  Pres¬ 
ence  of  a  Dislocation 

The  stress  fields  on  the  main  crack  due  to  this  dislocation  are  next  evaluated  by  using 
the  dislocation  stress  functions  (C.52,  C.53).  The  negative  of  this  stress  field  is  applied  on 
the  main  crack  and  the  stress  field  due  to  the  interaction  between  the  main  crack  and  the  dis¬ 
location  is  derived.  The  development  of  these  stress  functions  follows  the  method  formu¬ 
lated  by  Krenk  [52].  Following  the  formulation  of  Krenk  [52]  four  new  complex  stress  func¬ 


tions  are  defined  in  terms  of  the  known  stress  functions  as  follows. 

20i(zi)  =  OCzi)  +  0(Zi> 

202(Z2)  =  'P(Z2)  +  W(Z2)- 

2Qi(Zi)  =  W(zi)  -  W(zi> 

2^2(22)  =  W(Z2)  -  W(Z2)- 


(C.56) 

(C.57) 

(C.58) 

(C.59) 


The  normal  and  shear  stresses  can  be  defined  using  the  newly  introduced  stress  func¬ 


tions  as. 


Oyy  {0j(Zj)  "h  Q2(Zj)  -f"  02(Z2)  +  ^2(^2^ 
+  02(Zj)  ~  Q2(Zj)  -|-  02(Z2)  ~  Q2(Z2)} 


(C.60) 
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Oxy  ~  {Si0i(Zj)  +  S]Qj(Z])  +  $202(^2)  S2^2(^2) 

+  $30, (Zj)  -  S3Qj(Zj)  +  S402(Z2)  -  8402(^2)}  (C.61) 

Let  the  stresses  on  the  crack  surface  be  represented  by  the  following  functions, 


Figure  C.2.  Geometry  of  Pressurized  Main  Crack 
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Let  the  main  crack  extend  from  -c  to  +c  along  Y=0  axis.  Then, 

Oyy  =  {[0J  +  Qj  +  ©2  4  ^2]^  [®1  ~  ^1  ®2  ~  ^2^  ' 

Oyy  =  {[0J  “I"  £^2  "I"  ©2  4  ^2^  f®l  ~  ®2  ~  ^2^^'  (C.67) 

Oxy  ~  ~  {[Sj©]  “I"  SjQi  +  S2©2  "t"  ^2^2]^ 

—  [Sj©!  —  SjQj  +  S4©2  —  S4Q2]~}  (C.68) 

CJ^y  ”  "h  "H  S2€)2  4  ^2^^2^ 

—  [S3©i  -  S3Q1  +  S4©2  -  84^2]^} 

Therefore, 

f(t)  =  {[©1  +  ©2]+  +  [©1  +  ©2]">-  (C.70) 

“  h(t)  =  {[©jCXj  +  Q|i(3j^  +  ©2^2 

+  [©jCXj  +  Qjipi  +  ©20t2  "t”  ^2^p2^  (C.71) 

g(t)  =  -  {[Qj  +  Q2]+  -  [Q,  +  Q2]">- 

—  k(t)  =  {[Q|Otj  +  ©|iPj  +  £^2^2  ®2^p2^ 

4-  [QjOCj  +  ©jip^  +  Q2CX2  4  ®2^p2^^}  (C.73) 

where,  -c  <  t  <  c. 

Using  Muskhelishvili  [15],  the  most  general  solution  for  above  problem  can  be  written 


as, 


©j(z)  +  ©2(z) 


1 


I: 


f(t)  ^  P(z) 


2raX(z)  I  X  +  (t)(t  -  z)  X(z) 

h(t) 


(C.74) 


aj©,(z)  +  a2©2(z)  +  iPiQjCz)  +  ip2^2(z)  =  “  2jtiX(z) 


dt  + 


Q(z) 


Qj(z)  +  ^2(2)  ~ 


iX(z)  J 


g(t) 

2jtiX(z)  |^(t  -  z) 


^X  +  (t)(t-z)'^  X(z) 

(C.75) 

dt,  (C.76) 


aiQi(z)  +  a2Q2(z)  +  i|3i0i(z)  +  =  2^^  j  X+m- 

where, 


(C.77) 


n 

P(z)  =  ^  Anz"-  (C.78) 

n  =  0 
n 

Q(z)  =  ^  Bnz".  (C.79) 

n  =  0 

where  n  is  the  number  of  cuts  (i.e.  cracks).  Now  from  equation  C.56, 


20i(zi)  =  d)(z,)  +  a>(zi> 

(C.80) 

20i(zi)  =  d>(z,)  +  $(zj> 

(C.81) 

Therefore, 

0j(zi)  =  0i(zi). 

(C.82) 

Similarly, 

(C.83) 

Q|(z^)  =  ~  Qj(z^)- 

(C.84) 

Q2(Zi)  =  ^2(^1)* 

(C.85) 

Hence, 

+  ©2(2^)  —  0j(z^)  +  02(Zj)* 

(C.86) 

^1(^1)  ^2(^2)  “  ”  ^2(^2)* 

(C.87) 

From  equations  C.74  and  C.75, 

P(z)  =  P(z). 

(C.88) 

Q(z)  =  Q(z)- 

(C.89) 

Therefore  the  constants  in  P(z)  and  Q(z)  are  all  real. 


In  the  present  case  n  (the  number  of  cracks)  is  1  therefore,  from  equations  C.78  and 
C.79  there  are  4  undetermined  constants.  By  enforcing  the  condition  that  the  stresses  at  infin- 


ity  vanish,  Ai  and  B  i  vanishes.  Aq  and  Bq  have  to  be  determined  from  the  single  valuedness 
of  displacement  conditions.  Since  o'^yy=o“yy  and  o'^xy=crxy^  therefore  g(t)  =  k(t)  =  0.  The 
stress  functions  can  now  be  written  as, 


(S2  -  Si)^m(Zi)  = 


(Si  -  S2)^m(Z2)  = 


1 


2jtiX(Zj) 


X  +  (t)(s2f(t)  +  h(t))  S2A0  +  Bo 
—at  + 


(t  -  Z]) 


23TiX(Z2) 

where  -c  <  L  <  c  and 


X  +  (t)(Sif(t)  +  h(t))  .1^0 

-  — dt  + 


X(Zi) 
S,Ar,  +  B 


0 


(t  -  Z2) 


X(Z2) 


(C.90) 

(C.91) 


_  (Si  -  S2)fi  (S3  -  S2)fi  (S4  -  S2)f2 


ms,  +  h(t)  =  ^ 


03 


^04 


f(t)Si  +  h(t)  = 


(S2  -  Si)f2  ^  (S3  -  Si)fj  (S4  -  Si)f2 


t  —  Z 


+ 


02 


t  —  Z, 


+ 


03 


t  —  Z 


(C.92) 

(C.93) 


04 


where. 


_  1 


S]Bx  -  Bv 


fo  = 


2:ti(Sj  -  S2)(si  -  S3)(s,  -  34)011 
1  S2BX  —  By 


(C.94) 


(C.95) 


2m  (S2  -  Si)(S2  -  S3)(S2  -  84)011 

X(z)  is  a  sectionally  holomorphic  function  and  its  branch  cut  is  defined  as  shown  in 
Figure  C.3. 


z=te'®a 

z+c=rie‘®i 

z-c=r2e‘®2 

- ►  X 


Figure  C.3.  Geometry  for  the  Function  X(z) 


From  Figure  C.3, 


X(z)  = 
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(C.96) 

(C.97) 


X(z)  =  yrir2e'(02+®>)- 

The  branch  cut  for  z+c  is  taken  such  that,  0  <  0i  <  2k  and  similarly  the  branch  cut  for 
z-c  is  taken  such  that  0  <  02  <  2k.  Therefore, 

X  +  (t)  =  iv/c^  -  t2,  (C.98) 

X-(t)  =  -iyc2-t2.  (C.99) 

Using  the  methods  described  in  Appendix  A  (A.29  -  A.50)  equations  C.90  and  C.91 
can  be  evaluated  and  are  given  below. 

(H  -  Si)^m(z)  = 


4 


+ 


+ 


Similarly, 


4 


+  i 


'(Si  -  S2)fi 

X(Zoi)(Si  -  S2)fi 

(S,  -  S2)fi'| 

Z  -  Zqi 

X(z)(z  -  Zqi) 

X(z)  J 

(S3  -  82)?! 

X(Zo3)(S3  -  82)?! 

(S3  -  S2)f j 

Z  —  Zq3 

X(z)(z  -  Z03) 

X(z) 

J 

(S4  -  S2)f2 

X(Zo4)(S4  -  S2)f2 

(S4  -  82X2 

z  —  Zg4 

X(z)(z  -  Z04) 

X(z) 

^0  "b  Fq 

X(z) 

f(S2  -  Si)f2 

X(Zo2)(S2  -  Si)f2 

(S2  -  Si)f2l 

[  Z  ~  ^02 

x(z)(z  -  Z02) 

X(z)  J 

(S3  -  Sj)fi 

X(Zo3)(S3  -  Si)fi 

(S3  -  Si)fj 

Z  ~  ^03 

X(z)(z  -  Z03) 

X(z) 

> 

(S4  —  Sj)f2 

X(Zo4)(S4  -  Si)f2 

(S4  -  Si)f2 

Z  -  Zo4 

X(z)(z  -  Z04) 

X(Z) 

(C.lOO) 


+ 


2 

SiAo  +  Bo 

X(z) 


(C.lOl) 
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The  constants  Aq  and  Bq  are  determined  from  the  condition  that  the  displacements  are 
single  valued.  This  condition  can  be  expressed  as, 

j  2Re[(p,OD(Zl)  +  Pl<I'M(Zl))dZi  +  (P2^d(Z2)  +  P2^M(Z2))dZ2]  =  0-  (C.102) 

|2Re[(qiOD(Zi)  +  qi<DM(Zi))dZi  +  (q2^D(z2)  +  q2'I'M(z2))dz2]  =  0-  (C.103) 

where  the  contour  is  drawn  around  the  two  cuts  as  shown  in  Appendix  A  (Figure  A.5,  page 
95).  Shrinking  the  contours  around  the  cuts,  the  following  integrals  are  obtained. 


J  2Re[((pi<D  +  (t)  +  PiO+(t)  -  Pi<D5(t)  -  Pi<DM(t)) 

L,+L2 

+  (P2'PD(t)  +  P2^m(0  -  P2^D(t)  -  P2'I'M(t))]dt  =  0-  (C.104) 

I- 

J 

Lj  +  L2 

+  (q2W  +  (t)  +  q2W  +  (t)  -  qjW^Ct)  -  q2WM(t))]dt  =  0-  (C.105) 

Now  consider  equation  C.104.  Like  the  isotropic  case,  Om,  Wm  gives  different  limit¬ 
ing  values  on  Lj  only  and  Op  gives  different  limiting  values  on  L2  only.  Therefore,  using 

the  Plemelj  formulae, 


j[Pi^D(0  -  Pi<I>D(t)]dt  =  pi2rafj- 
Li  Li 


(C.106) 


(Si  -  S2)fi 
X+(t) 


X(Zo3)(S3  S2)fi  _  (S3  S2)fj  X(Zq4)(S3  —  S2)f2  (S3  S2)f2 

X  +  (t)(z-Z03)  X  +  (t)  X  +  (t)(z-Z04)  X+(t) 

S2A0  +  Bpj  Pjdt 
X  +  (t)  J(S2-Si) 


(Cl  07) 
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The  above  integrals  can  be  evaluated  using  contour  integrals  which  results  in  the  following 
expression. 


|[pi^mO)  -  piOM(t)]dt 
Li 


23tipj 
(§2  -  Sj) 


fl(Sj  -  S2)  +  fi(S3  -  S2)  +  f2(S4  -  S2)  -  (S2A0  + 

Similarly, 

r 

[P2^dW  ~  P2^D(t)]dt  =  P22jtif2- 


(C.108) 


(Cl  09) 


[P2^M(t)  “  P2^M(t)]dt  =  jf2(S2  -  Sj)  +  fi(S3  -  Sj) 


+  f2(S4  -  Si)  -  (SiAo  +  Bo)l^^" 


s?) 


(C.llO) 


Using  the  above  integrals  in  equation  C.104  the  following  equation  is  obtained. 

2Re(fi(S3  -  S2)Pi  +  f2(S4  -  S2)pi  -  fi(S3  -  Si)p2 

-  f2(S4  -  Si)p2  -  (S2A0  +  Bo)Pi  +  (SiAo  +  Bo)P2}  =  0  (C.lll) 

This  can  be  written  as, 

2Re<(fi  +  f2)(P2Si  -  p,S2)  +  (f,S3  +  f2S4)(Pi  -  Pj) 


+  Ao(p2Si  -  P1S2)  +  Bo(P2  -  Pi)>  =  0 


By  observation. 


Ao  =  -  (fi  +  f2) 


Bo  =  (Ssfj  +  S4f2) 


Since  Aq  and  Bq  are  real, 

Ao  =  -i(fi  +  f2  +  fi  +f2)  (C-n2) 

Bo  =  2^flSi  +  ^2^2 


(C.113) 
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Using  these  in  C.lOO  and  C.lOl  the  following  are  obtained. 


(S2  -  s,)<Dm(z)  =  ^  z-4 


2  z  -  Z03 


+  i 

21  Z-Z04 


X(Zoi)(Si 

-  S2)fl 

(s 

-  S2)fl' 

X(z)(z 

-  Zqi) 

X(Z) 

X(Zo3)(S3 

-  S2)fi 

(S3 

-  S2)fi 

X(z)(z  - 

"  Z;o3) 

X(Z) 

J 

X(Zo4)(S4 

~  S2)f2 

(S4 

-  82)^2' 

X(Z)(Z  - 

-  Z04) 

X(Z) 

> 

'3  “  S2)  + 

(2(84  - 

S2) 

2X(z) 


Similarly, 


(Cl  14) 


rc  —  C  Mir  _  X(Zq2)(S2  Sj)f2  (S2  Sj)f2 

^  ^  2 1  z  -  Z02  X(z)(z  -  Z02)  X(Z) 


+  i.  (S3  -  Si)fi  X(Zo3)(S3  -  Si)fi  (S3  -  Si)f^ 
2  z  -  Zo3  X(z)(z  -  Z03)  X(z) 


1  (S4  -  Si)f2  _  X(Zo4)(S4  -  Si)f2  _  (S4  -  Si)f^ 

2  I  z  -  Zo4  X(z)(z  -  Z04)  X(z) 

f2(S2  -  S|)  +  fi(S3  -  Sj)  +  f2(S4  -  $1) 

2X(z) 

Simplifying  the  above  expressions  the  following  are  obtained. 


(C.115) 


SiBx  Bv 


0^(z)=  j-j _ 1-  s^Bx-By  1  X(z,i) 

Cu4jti(s2  -  Sj)  I  (Sj  -  S3)(Sj  -  S4)  Z  -  Zoi  I  ^  X(z) 


S3BX  By 
(S3  —  S4)(S3  —  Sj) 


X(2o3) 

X(z) 


S4BX  By  J 

(§4  -  S3)(S4  -  Sj)  Z  -  Zo4 


X(Zo4) 


q.,„(z)  =  11  1  I  1  I  _ 

^  C]  1  4jri  (Sj  -  S2)  (S2  -  S3)(S2  -  S4)  Z  -  Z02  1  X(z) 


(C.116) 

X(Zo2)\ 
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SjBx  By 


1 


(S3  -  S4)(S3  -  S2)  Z  ^03 


By 


1 


1  - 


1  - 


X(Zo3)^ 

X(Z)  ^ 

X(Z,4)' 


(Cl  17) 


(S4  -  S3)(S4  -  S2)  2  ^04  \  ^  X(Z) 

The  equations  given  above  are  the  same  as  that  obtained  by  Obata  et.al  [54]  using  a  different 
procedure. 


Formulation  of  Stress  Functions  for  the  Main  Crack  Which  Opens  Up  Due  to  Loads  at 

Infinity 


The  stress  functions  for  the  problem  of  the  main  crack  with  remote  loadings  (C>l,  ^l) 
can  be  obtained  using  the  same  approach  as  discussed  in  the  previous  section  (see  also  Krerik 
[52])  and  are  given  below. 


Q(z)  +  S2P(z)  S2r'  +  A 

(S2-Si)X(z)  (S2-Si)’ 

(C.118) 

Q(")  +  SlP(z)  SlP'  +  A 

(Si  -  S2)X(Z)  (S,  -  S2)’ 

(C.119) 

where,  T  and  A  are  pure  imaginary  constants  which  have  to  be  determined  from  remote 

boundary  conditions  and 

n 

P(z)  =  X  AnZ" 

n  =  0 

(Cl  20) 

n 

Q(z)  =  y  B„z", 

(C121) 

n  =  0 

where  n  is  the  number  of  cuts  or  crack  and  An  and  Bn  are  real.  In  the  above  case  n  =  1  there- 

fore. 

P(z)  =  Ao  +  AjZ^ 

(Cl  22) 

Q(z)  =  Bo  +  BjZ- 

(Cl  23) 

Using  equation  C.  1 1 8  and  C.  1 1 9  in  the  expression  for  stresses  (C.35  to  C.37)  at  infinity 
(i.e.  z  tends  to  00)  the  following  relations  can  be  obtained. 
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rr” 

Ai 

(C.124) 

(^xv 

(C.125) 

2Re[(Si  +  S2)(A  -  Bj)  +  s,S2(r'  -  Aj)}  =  o~- 

(C.126) 

r'  =  ic, 

(C.127) 

A  =  iD- 

(C.128) 

Therefore, 

^  00  Q  q:  00 

C[0|P2  +  PiOj  +  D[P,  +  pj  =  ^  +  ^(a,  +  a,)  +  -f^Ca.aj  -  p.p^),  (C.129) 
where  C  or  D  have  to  be  chosen  arbitrarily  from  which  the  other  constant  is  determined.  Now 
Ao  and  Bo  have  to  determined  from  the  single  valuedness  of  displacement  condition  which 
can  be  expressed  as, 

j  ZReKp.l’^d)  -  Pitf  (t))  +  -  P2l'C(t)))dt 

L, 

where  Lj  e  {-a,  a}.  After  substituting  the  limiting  values  of  the  stress  functions  and  using 
the  following  integral  formulas  (See  Appendix  D,  D.26  -  D.38), 

f  dt  _ 

J 

L, 

1-^  =  0 

J  xut) 

L, 

Ao  =  0  and  Bo  =  0  is  obtained.  Assuming  C  =0,  the  stress  functions  for  remote  loading  can 
be  expressed  as. 


3)l(z)  =  2  Re 


+  SoO. 


2<Jyy  2 


2(S2  -  s,)  X(z) 


+  i 


.  +  ((Xja2  +  (cti  +  <^2)0^ 


2(S2  —  Sj)(pj  +  P2) 


(C.130) 
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Wl(z)  =  2Re| 


+  SiOy“  2  ^  .Ox”  +  («ia2  -  PlP2)qy”y  +  («i  +  «2)0'x°°y 


2(Si  -  S2)  X(z) 


2(Si  -  S2)(Pi  +  h) 


Since  in  the  present  study  a°°xy  =  0, 


<1)l(z)  =  2Re| 
^l(z)  =  2Re( 


s^o 


z'-’yy 


z  .Ox”  +  (aia2  -  PiP2)o 

^  +  1 


00  1 

yy 


2(s2  -  Si)X(z)  2(s2  -  Si)((3j  +  (32) 


^1*^^  Z  _j_  (tti0t2  PiP2)^: 


00  ' 

yy 


2(si  -  S2)X(z)  2(si  -  S2)(|3i  +  |32) 


(C.131) 

(C.132) 

(C.133) 


Formulation  of  the  Singular  Integral  Equation 
Bx(to)  and  By  (to)  can  be  treated  as  Green’s  function  for  a  point  dislocation  and  can  be 
integrated  over  a  given  range  to  represent  a  crack.  Then  the  only  boundary  condition  that 
remains  to  be  satisfied  is  the  stress  free  condition  on  the  branched  crack.  This  is  expressed 
as, 

O00(r,  0)  =  0, 


0re(r,0)  =  0- 

where,  0  <  r  <  /  and  0  is  the  kink  angle  and  I  is  the  kink  length.  Using  the  stress  transformation 
relations  given  below, 

O00  =  Oxx  cos^  0  +  Oyy  sin^  0  +  2axy  cos  0  sin  0i 

~  (%y  “  Oxx)  sin  0  cos  0  +  2Oxy(cos^0  —  sin^0)’ 

and  substituting  for  Oxx,  Oyy  and  Oxy  the  following  relations  are  obtained. 

[2  2" 

^>i(zi)(sj  sin0  +  COS0)  +  'Rj(z2)(s2  sin  0  +  cos0) 

a*g  =  2Re|'<l>j(zi)^(l  -  Sj)sin0cos0  -  Sicos20j 


+  ^j(z2)((l  —  S2)sin0cos0  —  S2cos20j 


®0e  “  ^60  *^00 
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where  i  =  D,  M  and  L.  Substituting  for  the  stress  functions  and  simplifying,  the  following 
pair  of  singular  integral  equations  are  obtained. 


i 

^69  ““j"' 


MiiBx(to)  +  Mi2By(to) 


2  2 

+  2 Re  (sjSinO  +  cos0)  ^>]^(zj)  +  (s2sin0  +  cos0) 


+  2Re  (sjsin0  +  cos0)^<I>l(zi)  +  (s2sin0  +  cos  0)^'Fl(z2) 


M2iBx(to)  "f  M22By(to) 


where. 


i 

+  J  2Re  ^(l  —  Sj)sin0cos0  -  Sj  cos20j<I)j^(Zj) 


+  ^(l  -  S2]sin0cos0  —  S2COS20 


+  2Re  1 


Sj)sin0cos0  —  Sj  cos20jOL(z2) 


+  f(l  —  S2)sin0cos0  —  S2COS20 


Zj  =  a  +  t(cos0  +  SiSin0), 


Zqj  =  a  +  to(cos0  +  SiSin0), 


and  i  =  1 ,  2.  Also 


(C.134) 

(C.135) 
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(§2  -  Si)(S2  -  S3)(S2  -  S4) 


(cos0  +  S2sin0) 


(C.136) 


Mn  =  —  2Re 


■(cos0  +  s.  sin0) 


^12  -  -  (S,  -  S2)(S,  -  S3)(S,  - 


(S2  -  Si)(S2  -  S3)(S2  -  S4) 


(cos0  +  S2sin0) 


(C.137) 


^21  2;riCii  ](si  -  S2)(si  -  83X81  -  84) 


1  —  82|8in0CO8  0  —  S|CO8  20 


(co8  0  +  $1  8in0) 


^1  —  s2jsin0co80  —  82COS20 

(82  —  Si)(82  —  S3)(S2  -  84)  (COS0  +  82  8m0) 


L  (C.138) 


M22  2Re  2ji;iCii"  (sj  -  S2XS1  -  83)(8i  -  84) 


^1  —  82jsin0CO80  -  8iCO820 
(co8  0  +  81  sin0) 


1  —  s2|sin0co8  0  —  80  cos  20 


(82  —  Si)(S2  —  S3)(S2  —  84)  (COS0  +  82  sin  0) 


.  (Cl 39) 
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APPENDIX  D 

EVALUATION  OF  SOME  CONTOUR  INTEGRALS 
Given  below  are  some  of  the  contour  integrals  which  occur  during  the  development 
of  the  elasticity  solution  discussed  in  chapter  II  and  III.  Consider  the  integral. 


A  =  f  (D-1) 

X  +  (t) 

J 

L 

where,  by  superscript  +  is  meant  the  limiting  value  of  function  X(t)  on  L  where  L  (-2c  < 
t  <  0)  is  defined  in  Figure  D.l  and 

X  +  (t)  =  yt(t  +  2c)-  (D.2) 


Figure  D.l .  Contour  Assumed  for  Evaluating  Integral  D.l 
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Using  Cauchy’s  integral  theorem  [89], 


'  d?  I  di 

J  X(^)  ^  f  X(?) 

r,  r, 


=  0- 


(D.3) 


Where  the  contour  is  shown  in  Figure  D.l .  Note  that  the  contours  F i  and  r2  are  of  opposite 
directions.  The  direction  assumed  for  r2  is  taken  to  be  positive.  The  outer  contour  r2  can 
be  extended  to  infinity  and  evaluated  and  similarly  the  inner  contour  T i  can  be  shrunk  on-to 
the  line  L  across  which  X(^)  is  discontinuous.  The  contours  on  Li  and  L2  cancel  each  other 
in  the  limit.  Consider  the  integral  on  r2  which  can  be  written  as 


I  -  0(i) 


dl. 


as  ^  tends  to  00.  Let, 

I  =  Re‘®,  d^  =  iRe*®d0. 

The  0(l/^2)  term  vanishes  as  ^  tends  to  00.  Therefore, 

2n 


1 


=1 

0 


d|  =  i  d0, 


=  2jTi- 

Consider  the  integral  on  which  can  be  written  as, 

0 

/■  .c  r 

dt 


T]  -2c 


+ 


d?  .  I  d^ 


+ 


x+(t)  ■  JX(^)  Jx© 

c,  C2 


I  d^ 

‘'x© 

J 


Let, 


I  =  Rc  +  2c,  d^  =  Rc  e“^c,id0 


c,’ 


(D.4) 

(D.5) 

(D.6) 

(D.7) 

(D.8) 

(D.9) 


where  Rci  is  the  radius  of  the  contour  Ci  and  0ci  is  the  angle  from  the  positive  direction  for 
the  radius  Rci.  Then  as  the  radius  Rci  tends  to  zero. 
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Similarly, 


f 

Cl 


f 


Rce'0c,id9c^ 


y(Rc  eie  +  2c)(Rc  e*  +  4c 


=  0(Rc,). 


(D.IO) 


Rc  ei®c2id0p 


Rc/^<^2)(Rc/^c,  +  2c) 


=  o  yRcj- 


(D.ll) 


where  Rc2  and  0c2  have  the  same  meaning  as  described  above.  Therefore  the  integrals  on 
contour  Ci  and  C2  vanish  as  their  radius  vanishes.  Hence 

0 


dt 


X+(t) 


—  —  2jXb 


(DA2) 


( 

I 


dt 


X  +  (t) 


-2c 


=  —  Jtl- 


(D.13) 


-2c 

Consider  the  integral 
0 


I 

-2c 


tdt 


X  +  (t) 


(D.14) 


where  X(t)  is  defined  as  before.  Let  the  contour  be  assumed  as  shown  in  Figure  D.  1 .  There¬ 
fore  as  before. 


ridi  iidi  ^ 

JX(i)'^JX(|)  ’ 

(D.15) 

Fi  r2 

(D.16) 

r2  r, 

(N 

1 

0 

11 

(D.17) 

The  integral  on  Ti  is  given  as. 
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I  =  7  I  tdt  ,  I ,  I 

JX(|)  ^Jx+(t)  Jx®  Jx®' 

T]  —2c  C,  C2 

As  before  the  integral  on  C]  and  C2  vanishes.  Therefore 


(D.18) 


X  +  (t) 


=  JtlC- 


(D.19) 


Consider  the  integral 


_ dt _ 

(t  -  Zo)X+(t)’ 


(D.20) 


where  X(t)  is  defined  as  before.  Let  the  contour  be  assumed  as  shown  in  Figure  D.2.  There¬ 
fore  as  before, 


-  Zo)X(^)  T  (^  -  Zo)X(^)  X(Zo)’ 


(D.21) 


The  right  hand  side  is  obtained  from  the  fact  that  there  exists  a  pole  at  Zqi  and  it  is  negative 
because  the  contour  enclosing  it  is  in  the  negative  direction.  Considering  the  contour  on  r2, 


(D.23) 


The  integral  on  T]  is  given  as, 


(?  -  Zo)X(§) 


V 

1^ 


_dU _ 

Zo)X  +  (t) 


Figure  D.3.  Contour  Assumed  for  Evaluating  Integral  D.26,  D.33 


As  before, 
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f  X(|)  +  f  X®  - 

r.  r, 


Consider  the  integral  on  r2  which  can  be  written  as, 


r 

=  2nh 

The  integral  on  Fj  is  given  as. 


0 


X(?) 


=  2 


tdt 


+ 


l£i+  Ijdl 


x"(t)  Tx(?)-  fm- 


As  before  the  integral  on  Ci  and  C2  vanishes.  Therefore 


C 

I 


dt 


X  +  (() 


=  -  Jtl- 


Consider  the  integral 


C 

I 

-c 


tdt 


X+(t) 


(D.28) 


(D.29) 

(D.30) 


(D.31) 


(D.32) 


(D.33) 


where  X(t)  is  defined  as  before.  Let  the  contour  be  assumed  as  shown  in  Figure  D.3.  There 
fore. 


lim 

|-»00 


=  0. 


(D.34) 

(D.35) 


The  integral  on  Fi  is  given  as. 


(D.36) 


0 
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APPENDIX  E 

NUMERICAL  TECHNIQUES. 

The  numerical  techniques  used  for  the  branched  crack  problem  is  similar  for  both  the 
isotropic  and  anisotropic  case.  Therefore,  the  solution  will  be  shown  for  a  general  case  and 
it  can  be  applied  for  either  the  isotropic  or  the  anisotropic  case.  Consider  the  integral, 
z  Z 

J-  +  I  M[a(to),t,to]dto  =  f(t),  (E.1) 

0  0 


where. 


a(to)  = 


P(to) 

t“(/  -  to)P 


(E.2) 


and  the  first  integral  is  interpreted  in  a  Cauchy’s  principal  value  sense.  The  function  a(to) 
(dislocation  density  function)  is  a  flux  type  function  which  can  be  expressed  in  terms  of  a 
bounded  function  P(to)  as  shown  in  equation  (E.2).  The  denominator  of  equation  E.2  is  also 
known  as  the  fundamental  function  (or  the  weight  function),  and  a  and  P  are  the  stress  singu¬ 
larities  at  the  kink  comer  and  tip  respectively.  M[a(to),  t,  to]  contains  the  terms  which  arise 
due  to  the  interactions  between  the  main  crack  and  the  dislocation  and  f(t)is  due  to  the  re¬ 
motely  applied  load. 

The  above  singular  integral  equation  is  solved  using  three  different  numerical  techniques  as 
given  in:  1,  Erdogan,  et  ah,  [90];  2,  He  and  Hutchinson  [91];  and  3,  Kaya  and  Erdogan  [85]. 
In  the  above  three  methods  the  limits  of  the  integral  equations  are  changed  so  that  the  limits 
are  from  -1  to  1  by  using  the  transformations  given  below. 


to  =  |(lo  +  1),  t  =  |(^  +  1),  dto  =  |dio- 


(E.3) 


135 

Method  I:  Solution  by  Gauss-Chebychev  Integration  Formula  [901 

Following  [90]  the  singular  integral  equation  is  reduced  to  a  system  of  linear  equations 
as  shown  below. 

=  f(ti),  (i  =  l,..n  -  1).  (E.4) 

where, 

Tn(toj)  =  0,toj  =  cos^jr^^2n'  '  ^ 

U„_,(ti)  =  0,  ti  =  cos(f),  (i  =  l,..n  -  1),  (E.6) 

Wj  =  f,  (j  =  l,-n>  (E-7) 

Tn(t)  and  Un_i  (t)  are  the  Chebychev  polynomial  of  the  first  and  second  kind  respectively  and 
Wj  are  the  weights.  This  results  in  (n-1 )  equations  in  n  unknowns.  For  the  complete  solution 
one  more  additional  condition  is  required  which  is  obtained  as  follows. 

In  method  I,  the  singularity  at  both  ends  of  the  branched  crack  or  kink  is  assumed  to  be  1/2. 
Since  the  singularity  at  the  kink  comer,  where  the  kink  joins  with  the  main  crack,  is  actually 
less  than  1/2,  the  value  of  the  bounded  function  P(to)  at  the  kink  comer  is  assumed  to  be  0. 
i.e.  P(-l)  =  0.  With  this  condition  the  above  system  of  linear  equations  can  be  completely 
solved.  This  method  was  used  by  Lo  [24]. 


j=l 


'-oj 


+  M 


Method  11:  Solution  by  Jacobi  Polynomials 
In  this  method  the  correct  singularity  is  assumed  at  the  kink  comer.  The  correct  singu¬ 
larity  is  given  by  the  Williams,  characteristic  equations  [19]  for  the  isotropic  case  and  by 
Bogy’s  [92]  equation  for  the  anisotropic  case.  This  method  is  similar  to  the  manner  in  which 
the  singular  integral  equation  was  solved  by  He  et.al.  [91]  though  in  their  approach  a  simple 
algebraic  polynomial  was  used  instead  of  Jacobi  polynomials.  The  advantage  of  using  Jaco¬ 
bi  polynomials  over  algebraic  polynomials  is  that  the  inherent  properties  of  Jacobi  polyno- 


mials  can  be  used  to  obtain  a  simplified  numerical  scheme.  The  unknown  function  P(to)  is 
expanded  in  terms  of  Jacobi  polynomials  as  shown  below. 


P(to)  =  ^CiP|“>P)(to)- 


(E.8) 


i  =  0 


a  and  p  are  the  stress  singularities  and  Q  is  a  complex  constant.  Substituting  the  above  equa¬ 
tion  in  equation  (E.l)  the  following  system  of  linear  algebraic  equation  is  obtained. 

n  m  p(a.P){'t  .'iW-  " 

't  ^  - i.--l  +  (E.9) 

i  =  0  j=l  '  i  =  0  j=l 

where  toj  is  the  root  of  Jacobi  polynomial  given  by  equation  (E.IO)  and  ti  are  given  by  the 
roots  of  the  Chebychev  polynomials  (equation  E.ll). 

Pm  =  0’  (E-10) 


ti  -  cos(  /  2)’  ^ 


Method  III:  Solution  bv  Double  Series  Jacobi  Polynomials 


The  bounded  function  P(to)  can  be  written  as,  (following  Kaya  et.al  [85]) 

Pij 


(E.ll) 


=  11773- 


,  oj  >  a.  Pi  >  Pj.  (E.12) 


-  t)Pj-Pi(l  +  t)^ 

If  (-l-tti+ai)  or  (-1-Pj+Pi)  <  0  then  the  derivative  of  p(to)  is  not  bounded  at  the  ends  of 
the  integration  limits.  If  that  is  the  case  then  expanding  a(to)  as  shown  below  will  improve 
the  rate  of  convergence  for  the  numerical  scheme. 


a(to)  = 


Pl(to) 


+ 


P2(to) 


-,ai  a,.(E.13) 


(1  -  to)t/2(l  +  to)“i  '  (1  -  to)t/2(l  -b  to)«2’"t  ^  2 

ai  and  a2  are  the  roots  of  the  Williams  characteristic  equation  [19]  for  the  isotropic  case  and 
the  roots  of  equation  given  by  Bogy  [92]  for  the  anisotropic  case.  Pi  (to)  and  P2(to)  are  ex¬ 
panded  in  terms  of  Jacobi  polynomials  as  in  Method  II  and  the  unknown  dislocation  density 
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function  is  evaluated.  The  collocation  points  for  this  method  are  again  given  by  the  roots 
of  the  Chebychev  polynomials  as  shown  below. 


tj  =  cos 


+  nt +  22)’  k-0,....,ni+n2+l- 


(E.14) 


ni  and  n2  represent  the  number  of  terms  in  the  Jacobi  polynomial  series  used  to  represent 
Pi  (to)  and  p2(to)  respectively. 


It  should  be  noted  that  where  as  Kaya,  et  al.,  [85]  used  the  above  method  for  apotential 
type  function  (i.e.  crack  opening  displacement)  this  current  method  is  used  for  flux  type 
function  where  the  unknown  is  the  slope  of  the  crack  opening  displacement.  The  above  three 
methods  are  also  referred  to  in  short  form  as  MI,  Mil  and  Mill  for  method  I,  method  II  and 
method  III  respectively. 
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